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

fix(prt): fix particle reporting #2185

Merged
merged 2 commits into from
Jan 31, 2025
Merged

Conversation

wpbonelli
Copy link
Contributor

@wpbonelli wpbonelli commented Jan 31, 2025

Reporting was incorrect for permanently unreleased particles. We want to report them once whenever they are to be released. #2177 fixed an issue where they were reported every time step, but the solution was incorrect, only reporting them if the release occurred in the first time step.

Also fix reporting for timed-out particles. The new status code 10 introduced by #2177 covers not only the simulation end but the particle reaching its stop time as well.

Update some comments and notes too.

Reporting was incorrect for permanently unreleased and timed-out particles. Also update some development notes.
Copy link
Contributor

@aprovost-usgs aprovost-usgs left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

nice job on this issue. just one last comment (about TODOs)

src/Model/ModelUtilities/TrackFile.f90 Show resolved Hide resolved
@wpbonelli wpbonelli marked this pull request as ready for review January 31, 2025 22:38
@wpbonelli wpbonelli merged commit 32c2e73 into MODFLOW-USGS:develop Jan 31, 2025
6 checks passed
@wpbonelli wpbonelli deleted the prt branch January 31, 2025 22:39
wpbonelli added a commit that referenced this pull request Feb 3, 2025
Use enums for particle status and event codes, following #2188. Also remove an unnecessary/incorrect termination condition left over from #2185 in the base MethodType%check() routine — stationary particles get terminated at the end of the simulation by the tracking loop with timeout status, previously status 5 (no exit face) was getting reported.
Copy link
Contributor

@aprovost-usgs aprovost-usgs left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

looks good

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

Successfully merging this pull request may close these issues.

2 participants