From 0dc1d7409f272aae1fa0773d1b91089b50a088b9 Mon Sep 17 00:00:00 2001 From: Carolyn Begeman Date: Tue, 9 Dec 2025 16:09:05 -0600 Subject: [PATCH] Fix sea ice time series cache masking bug --- mpas_analysis/sea_ice/time_series.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/mpas_analysis/sea_ice/time_series.py b/mpas_analysis/sea_ice/time_series.py index 9e4552072..8c4c9e0b7 100644 --- a/mpas_analysis/sea_ice/time_series.py +++ b/mpas_analysis/sea_ice/time_series.py @@ -695,8 +695,8 @@ def _compute_area_vol(self): ds.iceThick <= maxAllowedSeaIceThickness) if os.path.exists(outFileNames[hemisphere]): dsCache = xr.open_dataset(outFileNames[hemisphere]) - timeMask = ds.startTime > dsCache.startTime.isel(Time=-1) - ds = ds.isel(Time=timeMask) + timeMask = np.argwhere(ds.startTime.values > dsCache.startTime.isel(Time=-1).values) + ds = ds.isel(Time=timeMask[:,0]) dsAreaSum = (ds.where(mask) * dsMesh.areaCell).sum('nCells') dsAreaSum = dsAreaSum.rename(