interp at same precision as dataset

NOAA-PSL · Jul 23, 2024 · cbe6e00 · cbe6e00
1 parent 91a3cd4
commit cbe6e00
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/ufs2arco/fv3dataset.py b/ufs2arco/fv3dataset.py
@@ -96,7 +96,7 @@ def calc_pressure_interfaces(self, pressfc) -> xr.DataArray:
             prsi (xr.DataArray): 3D pressure field (or 4D if input has time) at vertical grid interfaces (Pa)
         """
 
-        return self.ak + pressfc*self.bk
+        return self.ak.astype(pressfc.dtype) + pressfc*self.bk.astype(pressfc.dtype)
 
     def calc_pressure_thickness(self, prsi) -> xr.DataArray:
         """Compute pressure thickness at each vertical level