-
Notifications
You must be signed in to change notification settings - Fork 3
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 support for discharge scaling in drainage package #1012
base: master
Are you sure you want to change the base?
Conversation
MODFLOW 6 has changed a number of options in the IMS package (since quite some time). * This adds support for outer_csvfile and inner_csvfile instead of csv_output. * This simplifies the no_ptc option: instead of a bool and a possible additional entry, just the entry is required. (In the IMS file, NO_PTC is equal to NO_PCT ALL, i.e. the ALL is an assumed default value.) * This also adds support for the adaptive time stepping entry ats_outer_maximum_fraction. All the deprecated solution keyword arguments have been removed from the examples.
Something to consider when dealing with the auxiliary arrays: MODFLOW 6 has no info in the .dat/.bin files to see which column is which. Instead, (I'm pretty sure) the meaning of the columns depends on the order of the auxiliary entry in the options block. This means it's crucial to coordinate the order of the variables when converting to a dictionary of arrays with the auxiliary entry written in the options block. Otherwise you end up with MF6 interpreting the concentration as the factor or vice versa. |
Quality Gate passedIssues Measures |
05ac3f7
to
5694d30
Compare
Fixes #1011
Still a draft, this adds the drainage discharge scaling depth.
This is a rather pragmatic approach. We're currently only using the auxiliary variables for concentration, which are treated separately. The
_auxiliary_data = {"concentration": "species"}
cannot be meaningfully used, since this specificscaling_depth
does not contain any additional dimensions.Adding it to
period_data
ensures that it is written to the data files, but since it's optional, the presences must then be checked inget_period_names
.I think the name "scaling_depth" is relatively clear. "discharge_scaling_depth" is maybe more complete, but is a bit of a mouthful.
We don't support auxiliary variables, I'm not sure if it's worth adding general support. Most packages support a multiplication factor auxiliary variable, which we really don't care for.