Update ex-gwf-csub-p03.py

scripts/ex-gwf-csub-p03.py

@@ -507,10 +507,19 @@ def build_models(
     opth = f"{name}.csub.obs"
     csub_csv = opth + ".csv"
     obs = [
-        ("cunit", "interbed-compaction", "cunit"),
-        ("aquitard", "interbed-compaction", "aquitard"),
-        ("nodelay", "interbed-compaction", "nodelay"),
-        ("delay", "interbed-compaction", "delay"),
+        ("cunit1", "interbed-compaction", (0,)),
+        ("cunit2", "interbed-compaction", (1,)),
+        ("cunit3", "interbed-compaction", (2,)),
+        ("aquitard6", "interbed-compaction", (5,)),
+        ("aquitard7", "interbed-compaction", (6,)),
+        ("aquitard8", "interbed-compaction", (7,)),
+        ("nodelay4", "interbed-compaction", (3,)),
+        ("nodelay5", "interbed-compaction", (4,)),
+        ("nodelay9", "interbed-compaction", (8,)),
+        ("nodelay10", "interbed-compaction", (9,)),
+        ("delay11", "interbed-compaction", (10,)),
+        ("delay12", "interbed-compaction", (11,)),
+        ("delay13", "interbed-compaction", (12,)),
         ("es14", "estress-cell", (nlay - 1, 0, 0)),
     ]
     for k in (1, 2, 3, 4, 6, 7, 8, 9, 11, 13):
@@ -793,7 +802,9 @@ def process_csub_obs(fpth):
     # transfer data from temporary storage
     for key in pcomp:
         if key != "TOTAL" and key != "SKELETAL":
-            v[key] = tdata[key].copy()
+            for obs_key in tdata.dtype.names:
+                if key in obs_key:
+                    v[key] += tdata[key]
 
     # calculate skeletal
     for key in tdata.dtype.names[1:]: