From 4d69d4b138b29d6f6fb955c8112b3f691455f560 Mon Sep 17 00:00:00 2001 From: Xylar Asay-Davis Date: Sun, 25 Sep 2022 14:30:45 -0500 Subject: [PATCH] Don't decode time as part of transect remapping xarray seems to have trouble with masking a dataset with a `cftime._cftime.DatetimeNoLeap` time coordinated. To avoid this, we'll just not decode time since we don't need to. --- mpas_analysis/ocean/compute_transects_subtask.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/mpas_analysis/ocean/compute_transects_subtask.py b/mpas_analysis/ocean/compute_transects_subtask.py index 16276f7e4..58f4e7552 100644 --- a/mpas_analysis/ocean/compute_transects_subtask.py +++ b/mpas_analysis/ocean/compute_transects_subtask.py @@ -262,7 +262,7 @@ def run_task(self): remappedFileName = self.get_remapped_file_name( season, comparisonGridName=self.transectCollectionName) - ds = xr.open_dataset(remappedFileName) + ds = xr.open_dataset(remappedFileName, decode_times=False) transectNames = list(obsDatasets.keys()) for transectIndex, transectName in enumerate(transectNames): self.logger.info(' {}'.format(transectName))