Cast back to float64 when reading from the HDF5 file

Author: cphyc

yt/frontends/dyablo/data_structures.py

@@ -302,7 +302,7 @@ def fill(self, fields, selector):
                 field_path = f"/{fname}"
                 if field_path in f:
                     # Apply cell-level mask
-                    data[ftype, fname] = f[field_path][:][indices]
+                    data[ftype, fname] = f[field_path][:][indices].astype("float64")
                 else:
                     raise KeyError(
                         f"Field {field_path} not found in HDF5 file {self.ds._hydro_filename}"