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

Add support for discharge scaling in drainage package #1012

Draft
wants to merge 13 commits into
base: master
Choose a base branch
from

Conversation

Huite
Copy link
Contributor

@Huite Huite commented May 6, 2024

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 specific scaling_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 in get_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.

@Huite Huite marked this pull request as draft May 6, 2024 14:11
@Huite Huite requested a review from JoerivanEngelen May 6, 2024 14:11
Huite added 9 commits May 7, 2024 14:50
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.
@Huite
Copy link
Contributor Author

Huite commented May 8, 2024

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.

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.

Add support for discharge scaling depth in drainage package
1 participant