aquifer_summary valueerror when more than two layers and topboundary semi

[tomvansteijn]

Example code

ml = tmf.steady.Model3D(z=[10., 0., -10.], topboundary="semi")
ml.aquifer_summary()

Output

---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
Cell In[277], [line 2](vscode-notebook-cell:?execution_count=277&line=2)
      1 ml = tmf.steady.Model3D(z=[10., 0., -10.], topboundary="semi")
----> [2](vscode-notebook-cell:?execution_count=277&line=2) ml.aquifer_summary()

File c:\Temp\conda-envs\bw_basis\Lib\site-packages\timflow\steady\model.py:592, in Model.aquifer_summary(self)
    590 aqs = {}
    591 if not isinstance(self.aq, SimpleAquifer):
--> [592](file:///C:/Temp/conda-envs/bw_basis/Lib/site-packages/timflow/steady/model.py:592)     aqs["background"] = self.aq.summary()
    593 for i, iaq in enumerate(self.aq.inhomlist):
    594     aqs[f"inhom{i}"] = iaq.summary()

File c:\Temp\conda-envs\bw_basis\Lib\site-packages\timflow\steady\aquifer.py:183, in AquiferData.summary(self)
    181     else:
    182         summary.iloc[~maskaq, 2] = self.Hll
--> [183](file:///C:/Temp/conda-envs/bw_basis/Lib/site-packages/timflow/steady/aquifer.py:183)         summary.iloc[~maskaq, 4] = self.c
    184 summary.loc[:, "layer"] = self.layernumber
    185 return summary

File c:\Temp\conda-envs\bw_basis\Lib\site-packages\pandas\core\indexing.py:912, in _LocationIndexer.__setitem__(self, key, value)
    909 self._has_valid_setitem_indexer(key)
    911 iloc = self if self.name == "iloc" else self.obj.iloc
--> [912](file:///C:/Temp/conda-envs/bw_basis/Lib/site-packages/pandas/core/indexing.py:912) iloc._setitem_with_indexer(indexer, value, self.name)
...
   2004 elif lplane_indexer == 0 and len(value) == len(self.obj.index):
   2005     # We get here in one case via .loc with a all-False mask
   2006     pass

ValueError: Must have equal len keys and value when setting with an iterable