Add valid range before doing remapycon when processing lpjg tp avoid …

…problems with NaNs, Infs, etc (see subject of #823).
EC-Earth · Apr 15, 2024 · fdc8d79 · fdc8d79
1 parent 9acc3c6
commit fdc8d79
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/ece2cmor3/lpjg2cmor.py b/ece2cmor3/lpjg2cmor.py
@@ -518,10 +518,10 @@ def create_lpjg_netcdf(freq, inputfile, outname, outdims):
     cdo = Cdo()
 
     if target_grid_ == "T159":
-        cdo.remapycon('n80', input="-setgrid," + gridfile_ + " " + temp_ncfile,
+        cdo.remapycon('n80', input="-setvrange,-1e16,1e16 -setgrid," + gridfile_ + " " + temp_ncfile,
                       output=ncfile)
     else:
-        cdo.remapycon('n128',input="-setgrid," + gridfile_ + " " + temp_ncfile,
+        cdo.remapycon('n128',input="-setvrange,-1e16,1e16 -setgrid," + gridfile_ + " " + temp_ncfile,
                       output=ncfile)  # TODO: add remapping for possible other grids
 
     os.remove(temp_ncfile)