[meteo] account for descending lats (netcdf) + 0-360 lons (grib/netcdf)

tomsail · Jan 11, 2024 · bebe775 · bebe775
1 parent 9635d91
commit bebe775
Showing 1 changed file with 13 additions and 1 deletion.
diff --git a/pyposeidon/meteo.py b/pyposeidon/meteo.py
@@ -184,6 +184,10 @@ def cfgrib(
             v10 = data.v10[~mask_]
         data = xr.merge([msl, u10, v10])
 
+    convert360 = kwargs.get("convert360", False)
+    if convert360:
+        lon_min, lon_max = 0, 360
+
     if not lon_min:
         lon_min = data.longitude.data.min()
     if not lon_max:
@@ -441,7 +445,15 @@ def netcdf(
         }
     )
 
-    data = data.sel(longitude=slice(lon_min, lon_max)).sel(latitude=slice(lat_min, lat_max))
+    convert360 = kwargs.get("convert360", False)
+    if convert360:
+        lon_min, lon_max = 0, 360
+
+    descending_lats = data.latitude.values[0] > data.latitude.values[1]
+    if descending_lats:
+        data = data.sel(longitude=slice(lon_min, lon_max), latitude=slice(lat_max, lat_min))
+    else:
+        data = data.sel(longitude=slice(lon_min, lon_max), latitude=slice(lat_min, lat_max))
 
     try:
         data = data.sel(time=slice(start_date, end_date))