convert single to double precision hdf5 data

webgeodyn/inout/pygeodyn_hdf5.py

@@ -33,9 +33,10 @@ def load(dataDirectory, dataModel, keepRealisations, state_type='analysed'):
     with h5py.File(hdf_filename) as hdf_file:
         computed_data = hdf_file[state_type]
 
-        times = np.array(computed_data['times'])[firstpoint:]
+        times = np.array(computed_data['times'], dtype="float64")[firstpoint:]
 
         for measureName, measureData in computed_data.items():
+            measureData = np.float64(measureData) # Convert to float64
             if measureName not in measures_to_load:
                 continue
             else: