bfModel hdf5 - Githubissues

Describe the bug non-existent flowAcc data

To Reproduce run bfModel notebook udpated

Expected behavior vtk output

Desktop (please complete the following information):

OS: WSL

Additional context

EDrate
elev
erodep
fillFA
hdisp
rain
sedLoad
waterFill
---------------------------------------------------------------------------
KeyError                                  Traceback (most recent call last)
Cell In[9], line 2
      1 # Reading the final output generated by gospl
----> 2 output = rout.readOutput(filename='forward.yml',step=15)
      4 # Exporting the final output as a VTK mesh
      5 output.exportVTK('step15.vtk')

File /mnt/c/users/rscott/OneDrive - BHP/exploration2024/gospl/notebooks/bfModel/scripts/readOutput.py:33, in readOutput.__init__(self, filename, step)
     29 self._inputParser()
     31 self.nbCPUs = len(glob.glob1(self.outputDir + "/h5/", "topology.p*"))
---> 33 self._readElevationData(step)
     35 return

File /mnt/c/users/rscott/OneDrive - BHP/exploration2024/gospl/notebooks/bfModel/scripts/readOutput.py:298, in readOutput._readElevationData(self, step)
    296 def _readElevationData(self, step):
--> 298     self._getCoordinates(step)
    300     return

File /mnt/c/users/rscott/OneDrive - BHP/exploration2024/gospl/notebooks/bfModel/scripts/readOutput.py:105, in readOutput._getCoordinates(self, step)
    103 erodep = np.array((df2["/erodep"]))
    104 sedLoad = np.array((df2["/sedLoad"]))
--> 105 flowAcc = np.array((df2["/flowAcc"]))
    106 uplift = np.array((df2["/uplift"]))
    107 hdisp = np.array((df2["/hdisp"]))

File h5py/_objects.pyx:54, in h5py._objects.with_phil.wrapper()

File h5py/_objects.pyx:55, in h5py._objects.with_phil.wrapper()

File ~/miniconda3/envs/gospl-package/lib/python3.9/site-packages/h5py/_hl/group.py:357, in Group.__getitem__(self, name)
    355         raise ValueError("Invalid HDF5 object reference")
    356 elif isinstance(name, (bytes, str)):
--> 357     oid = h5o.open(self.id, self._e(name), lapl=self._lapl)
    358 else:
    359     raise TypeError("Accessing a group is done with bytes or str, "
    360                     "not {}".format(type(name)))

File h5py/_objects.pyx:54, in h5py._objects.with_phil.wrapper()

File h5py/_objects.pyx:55, in h5py._objects.with_phil.wrapper()

File h5py/h5o.pyx:241, in h5py.h5o.open()

KeyError: "Unable to synchronously open object (object 'flowAcc' doesn't exist)"

Geodels / gospl

bfModel hdf5 #274