Skip to content

Writing a wel package in a modflow model doesn't work properly #2224

Closed Answered by jlarsen-usgs
YANG-814 asked this question in Q&A
Discussion options

You must be logged in to vote

Have you opened your modflow list file and looked to see if it contains any error messages (e.g., model did not converge, input errors, etc...)

I'd also note that the stress_period_data is usually provided in a dictionary to boundary condition packages, here's an example:

nper = mf.dis.nper
pumping = [
[76 - 1,99 - 1,1 - 1,-10],
[81 - 1,98 - 1,1 - 1,-10],
[85 - 1,102 - 1,1 - 1,-10],
]
wel_spd = {i: pumping for i in range(nper)}
wel = flopy.modflow.ModflowWel(mf, stress_period_data=wel_spd)
mf.write_input()

Your pumping record also looks odd to me. Do you have 85 or more layers in your model?

Replies: 1 comment 10 replies

Comment options

You must be logged in to vote
10 replies
@wpbonelli
Comment options

@YANG-814
Comment options

@wpbonelli
Comment options

@jlarsen-usgs
Comment options

Answer selected by wpbonelli
@YANG-814
Comment options

@YANG-814
Comment options

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Category
Q&A
Labels
None yet
3 participants