Changeset 4249

Timestamp:

01/28/19 15:29:51 (5 years ago)

Author:

cetlod

Message:

CPLRestart: update cplrestart tools for New NEMO files

Location:

TOOLS/CPLRESTART

Files:

• CreateRestartOce4Oasis.bash (modified) (5 diffs)
• FillOceRestart.py (modified) (1 diff)

TOOLS/CPLRESTART/CreateRestartOce4Oasis.bash

@@ -39 +39 @@
 # Default values
 Comment="Preliminary attempt - Do not trust !"
-IceVar=ice_pres # siconc #> Variable containing sea ice fraction [0-1]
+##IceVar=ice_pres # siconc #> Variable containing sea ice fraction [0-1]
+IceVar=siconc #> Variable containing sea ice fraction [0-1]
 
 #
@@ -99 +100 @@
 ## Defines associated sea ice file
 ## ===========================================================================
-IceFile=$(echo ${OceFile} | sed 's=/OCE/=/ICE/=' | sed 's=grid_T=icemod=')
+###IceFile=$(echo ${OceFile} | sed 's=/OCE/=/ICE/=' | sed 's=grid_T=icemod=')
+IceFile=${OceFile}
 
 ##
@@ -114 +116 @@
 ## ===========================================================================
 dim_y=$(ncdump -h sstoce_fields_notime.nc | grep "y *=" | awk '{print $3}' )
-dim_x=$(ncdump -h sstoce_fields_notime.nc | grep "x *=" | awk '{print $3}' )
+dim_x=$(ncdump -h sstoce_fields_notime.nc | grep "x *=" | grep -v "nvertex" | awk '{print $3}' )
 echo ${dim_x} ${dim_y}
 
 [[ ${dim_x} = 182 && ${dim_y} = 149 ]] && OCE=ORCA2.3
 [[ ${dim_x} = 362 && ${dim_y} = 332 ]] && OCE=eORCA2.1
+[[ ${dim_x} = 1442 && ${dim_y} = 1207 ]] && OCE=eORCA025.1
 
 
@@ -126 +129 @@
 
 cat <<EOF > create_sstoce.nco
-OIceFrc        = double ( ${IceVar}(:,:) ) ; 
+OIceFrc[y,x]   = double ( ${IceVar}(:,:) ) ; 
 //tos.simple_fill_miss(tos)               ;
-O_SSTSST[y,x]  = double ( (tos   (:,:) + 273.15d) * (1.0d-OIceFrc) ) ;
+O_SSTSST[y,x]  = double ( (tos   (:,:) + 273.15d) * (1.0d-OIceFrc(:,:)) ) ;
 //O_AlbIce[y,x]  = double ( ialb   (:,:)            * OIceFrc(:,:)   ) ;
 //O_TepIce[y,x]  = double ( (tsice (:,:) + 273.15d) * OIceFrc(:,:)   ) ;
@@ -139 +142 @@
 
 ncap2 --overwrite --fl_fmt=${FL_FMT} --history --script-file create_sstoce.nco sstoce_fields_notime.nc tmp_sstoc.nc
-
 ncks --fl_fmt=${FL_FMT} --overwrite --history --variable OIceFrc,O_SSTSST,O_AlbIce,O_TepIce,O_OCurx1,O_OCury1,O_OCurz1 tmp_sstoc.nc sstoc.nc
-
 ncatted --history --attribute comment,O_SSTSST,o,c,"SST weighted by fraction of open ocean"          sstoc.nc
-ncatted --history --attribute comment,O_SSTSST,o,c,"Albedo weighted by fraction of sea ice"          sstoc.nc
-ncatted --history --attribute comment,O_SSTSST,o,c,"Ice temperature weighted by fraction of sea ice" sstoc.nc
+ncatted --history --attribute comment,O_AlbIce,o,c,"Albedo weighted by fraction of sea ice"          sstoc.nc
+ncatted --history --attribute comment,O_TepIce,o,c,"Ice temperature weighted by fraction of sea ice" sstoc.nc
 
 ##

TOOLS/CPLRESTART/FillOceRestart.py

@@ -148 +148 @@
         nperio = 6
         print ("eORCA1 grid found from dimensions, nperio=6" )
+    elif (jpoj, jpoi) == (1021, 1442) :
+        nperio = 4
+        print ("ORCA025 grid found from dimensions, nperio=6" )
+    elif (jpoj, jpoi) == (1207, 1442) :
+        nperio = 4
+        print ("eORCA025 grid found from dimensions, nperio=6" )
         
 if nperio == None :