From cbe6e00a05be651009642832dfb81b920e3542b1 Mon Sep 17 00:00:00 2001 From: timothyas Date: Tue, 23 Jul 2024 18:11:30 +0000 Subject: [PATCH] interp at same precision as dataset --- ufs2arco/fv3dataset.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/ufs2arco/fv3dataset.py b/ufs2arco/fv3dataset.py index 316fdad..33087c0 100644 --- 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