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
DISV model produced by gridgen not converging when original delr/delc contain decimals #1492
Comments
@dbrakenhoff, thanks for this. I'll take a look. I'm guessing that due to rounding errors the code is creating two vertices very close to one another rather than generating a single vertex. Just a hunch. We also noticed that the MODFLOW 6 DISV package will accept a line as a valid cell, and so we are planning to add a trap for that. |
This comment was marked as off-topic.
This comment was marked as off-topic.
This comment was marked as off-topic.
This comment was marked as off-topic.
Ignore my prior comments, they refer to unrelated issues. Reproduced in #2119 |
Confirmed @langevin-usgs' theory — flopy's gridgen driver should catch these duplicates and remove them |
When the
delr
anddelc
of the original structured grid are not whole numbers, the resulting DISV grid produced by gridgen can cause convergence issues when trying to run the MF6 model. I couldn't find out if there were requirements for the original grid for gridgen, but since it a produces a seemingly valid DISV grid, this seems like a bug?This example showcases the issue, adapted from the example notebook flopy3_gridgen.ipynb. When
delr
/delc
are set to[1.1, 1.2, 1.3, 1.4, 1.6, 1.7, 1.8, 1.9]
the model fails to converge, but[1.25, 1.5, 1.75]
work fine.The text was updated successfully, but these errors were encountered: