-
Notifications
You must be signed in to change notification settings - Fork 6
Issue #598 layered hfb error #1089
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
Issue #598 layered hfb error #1089
Conversation
not more than 1 layer assigned to LayeredHFB. Also remove unused validation to_mf6_pkg.
luitjansl
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
looks good, some comments
imod/mf6/hfb.py
Outdated
| compound_dataframe = pd.concat(dataframe_ls, ignore_index=True) | ||
|
|
||
| return LayeredHorizontalFlowBarrierResistance(compound_dataframe) | ||
| return cls(compound_dataframe) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
there is also a "LayeredHorizontalFlowBarrierHydraulicCharacteristic".
Should the changes in this PR also be applied there?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I think not. The hydraulic characteristic is a MODFLOW 6 concept: it is the inverse of the resistance. In the end, both the Resistance as well as the HydraulicCharacteristic variant are converted to a hydraulic characteristic for the Modflow6 HFB package. iMOD5 only allows you to specify a hfb with a resistance, so I don't see a real use case for adding a from_imod5_data method to the HydraulicCharacteristic class.
Fixes #598
Description
There was some confusion (also by me) about the use of the
LayeredHorizontalFlowBarrierclasses, so @HendrikKok explained me they were intended for single layers. I therefore changed the following to clarify:SingleLayerHorizontalFlowBarrierclassesValidationErrorif multiple layers are added to datasetThis means multiple SingleLayerHorizontalFlowBarrier need to be added to GroundwaterFlowModel object, whereas MODFLOW 6 only accepts one single HFB. So these have to be merged together. I think the best opportunity for this is to merge the low level
Mf6HorizontalFlowBarrierobjects.Checklist
Issue #nr, e.g.Issue #737