Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Bug in LAK if lakes not active in first stress period #40

Open
langevin-usgs opened this issue May 21, 2020 · 0 comments
Open

Bug in LAK if lakes not active in first stress period #40

langevin-usgs opened this issue May 21, 2020 · 0 comments

Comments

@langevin-usgs
Copy link
Contributor

The attached model from Diganta Sarma demonstrates a bug in the Lake package. This is a MODFLOW-NWT model but it looks like the same issue could be a problem in MODFLOW-2005. In the model, lakes are not active in the first stress period so in GWF2LAK7RP, it skips to the line labeled 900. It then tries to compute the conductances for the lake. However, the lake leakances were never read and MODFLOW crashes in subroutine SGWF2LAK7UPW1RPS. It seems to me that the way to fix this might be to add a continue statement labeled 900 after the series of calls in GWF2LAK7RP where the conductances are calculated and remove the 900 label from the statement that has it now.
LakeAqEqSim1.zip

Niswonger sent these fixes.
lakfix.zip

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant