set OC to save only specific layers #2104
-
Im currently trying to run a transient MF6 simulation where I want to save only head data from specific layers, for example only layer 1 and 2 (the model has 34 layers) to the .HDS file. I'm using flopy.mf6.ModflowGwfoc to set the OC-package and have tried to set the output to "Head 1 2" to save only layer 1 and 2. However, this results in a unsuccesful run, generating the follow error report:
Does anyone have an idea on how to go about this? Thanks in advance! The oc-package: _package_name = oc Block optionsbudget_filerecord head_filerecord Block periodsaverecord printrecord |
Beta Was this translation helpful? Give feedback.
Replies: 1 comment 1 reply
-
Hi @RdWitte, This is unfortunately possible in the OC package. From the MF6 input/output guide:
So the text must be "HEAD" and cannot contain any layer information. A workaround could be to create an Observations package that saves the heads for the cells of the layers you are interested in. For more information see the section on "Observation (OBS) Utility". This observation package can be added to your model with I hope this helps! Cheers, Davíd |
Beta Was this translation helpful? Give feedback.
Hi @RdWitte,
This is unfortunately possible in the OC package. From the MF6 input/output guide:
So the text must be "HEAD" and cannot contain any layer information.
A workaround could be to create an Observations package that saves the heads for the cells of the layers you are interested in. For more information see the section on "Observation (OBS) Utility". This observation package can be added to your model with
flopy.mf6.ModflowUtlobs
. This example in the flopy docs also shows how to create one of these …