-
Notifications
You must be signed in to change notification settings - Fork 116
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
Add BOUNDNAME of terminating sink/source to Particle Track Output #1851
Comments
Hey @cneyens, thanks for this. @aprovost-usgs and I gave it some thought. We see the utility, but there are complications. If a cell has multiple boundary conditions any of which do not have To support this, we'd need to keep face flow constituents distinct. I don't yet have a precise sense for the level of effort involved, but safe to say it's a significant refactor, even if the wider MF6 framework can accommodate it at the moment (still TBD). And the proposed extra column would need a special value indicating when a terminating boundary can't be determined. Or alternatively we'd need to report all candidates. Let us know if you see better ways around this. |
Hi @wpbonelli, that does sound more complicated than I initially thought! We were looking for a way to see which boundary condition terminates a particle. The IZONE functionality gets us there for the most part, so we'll use that. Thanks for the response! |
Is your feature request related to a problem? Please describe.
We've been playing around with the new PRT model in MODFLOW 6. As a feature request, I was wondering if it would be possible to include the BOUNDNAME of the terminating sink/source boundary condition for a given particle in the output file, if a boundname is given.
Describe the solution you'd like
An extra column in the Particle Track Output file containing the BOUNDNAME (if any) of the boundary condition which terminates a particle. While the particle is still active, this column could remain empty or contain a placeholder string.
Describe alternatives you've considered
Using the IZONE functionality. This requires setting up a IZONE array a priori, and also does not distinguish between different boundary conditions within a cell which might terminate the particle trace, e.g. at different cell faces. Lastly, we could always do this in postprocessing by comparing the endpoint locations and the boundary condition locations.
The text was updated successfully, but these errors were encountered: