From 73288741621ff2831a95373357407784211db601 Mon Sep 17 00:00:00 2001 From: Jiwoo Lee Date: Tue, 3 Oct 2023 11:17:55 -0700 Subject: [PATCH] resolve bug (error occuring in case coords does not have axis attributes) --- pcmdi_utils/land_sea_mask.py | 10 ++++++++-- 1 file changed, 8 insertions(+), 2 deletions(-) diff --git a/pcmdi_utils/land_sea_mask.py b/pcmdi_utils/land_sea_mask.py index e7a6d2c..d708694 100644 --- a/pcmdi_utils/land_sea_mask.py +++ b/pcmdi_utils/land_sea_mask.py @@ -472,10 +472,16 @@ def generate_land_sea_mask__global_land_mask(ds, maskname="lsmask"): """ for key in list(ds.coords.keys()): - if key in ["lat", "latitude"] or ds[key].attrs["axis"] == "Y": + if key in ["lat", "latitude"]: lat_key = key - elif key in ["lon", "longitude"] or ds[key].attrs["axis"] == "X": + elif key in ["lon", "longitude"]: lon_key = key + else: + if "axis" in list(ds["latitude"].attrs.keys()): + if ds["latitude"].attrs["axis"] == "Y": + lat_key = key + elif ds["latitude"].attrs["axis"] == "X": + lon_key = key lat = ds[lat_key].to_numpy() lon = ds[lon_key].to_numpy()