add workaround for flopy ncelldim shape issue in exg dfns

MODFLOW-USGS · Dec 18, 2023 · 9b82be9 · 9b82be9
1 parent facc832
commit 9b82be9
Show file tree

Hide file tree

Showing 3 changed files with 4 additions and 4 deletions.
diff --git a/doc/mf6io/mf6ivar/dfn/exg-gwfgwf.dfn b/doc/mf6io/mf6ivar/dfn/exg-gwfgwf.dfn
@@ -237,7 +237,6 @@ description
 block exchangedata
 name cellidm1
 type integer
-shape (ncelldim)
 in_record true
 tagged false
 reader urword
@@ -249,7 +248,6 @@ numeric_index true
 block exchangedata
 name cellidm2
 type integer
-shape (ncelldim)
 in_record true
 tagged false
 reader urword

diff --git a/doc/mf6io/mf6ivar/dfn/exg-gwtgwt.dfn b/doc/mf6io/mf6ivar/dfn/exg-gwtgwt.dfn
@@ -198,7 +198,6 @@ description
 block exchangedata
 name cellidm1
 type integer
-shape (ncelldim)
 in_record true
 tagged false
 reader urword
@@ -210,7 +209,6 @@ numeric_index true
 block exchangedata
 name cellidm2
 type integer
-shape (ncelldim)
 in_record true
 tagged false
 reader urword

diff --git a/utils/idmloader/scripts/dfn2f90.py b/utils/idmloader/scripts/dfn2f90.py
@@ -332,6 +332,10 @@ def _set_blk_param_strs(self, blockname, component, subcomponent):
 
             shape = ""
             shapelist = []
+            # workaround for Flopy shape issue with exg dfns:
+            if c.upper() == "EXG":
+                if vn == "CELLIDM1" or vn == "CELLIDM2":
+                    v["shape"] = "(ncelldim)"
             if "shape" in v:
                 shape = v["shape"]
                 shape = shape.replace("(", "")