Changeset 4186 for TOOLS/MOSAIX/CreateOasisGrids.bash

Timestamp:

12/12/18 11:26:20 (5 years ago)

Author:

omamce

Message:

O.M. :

Runoff weight for both ICO and LMDZ
Correct corc grid

File:

• TOOLS/MOSAIX/CreateOasisGrids.bash (modified) (9 diffs)

TOOLS/MOSAIX/CreateOasisGrids.bash

@@ -54 +54 @@
 PROGRAM=$(basename ${0})
 
-case ${arch} in
-    ( curie | irene )
-    set +vx
-    module unload cdo nco ferret
-    module unload netcdf hdf5
-    module load python # /2.7.12
-    module load netcdf/4.3.3.1_hdf5_parallel # Version for XIOS
-    module load nco
-    R_IN=$(ccc_home -u igcmg --cccwork)/IGCM
-    TMPDIR=${SCRATCHDIR}/TMP
-    SUBMIT_DIR=${BRIDGE_MSUB_PWD:-${SUBMIT_DIR}}
-    ;;
-    ( spip )
-    R_IN=${HOME}/Scratch/IGCM
-    TMPDIR=${HOME}/Scratch/TMP
-    SUBMIT_DIR=$(pwd)
-    MPIRUN=/opt/local/bin/mpirun-openmpi-gcc49 -n 2
-    ;; 
-    ( * ) exit -1 ;;
-esac
+# case ${arch} in
+#     ( curie | irene )
+#     set +vx
+#     module unload cdo nco ferret
+#     module unload netcdf hdf5
+#     module load python # /2.7.12
+#     module load netcdf/4.3.3.1_hdf5_parallel # Version for XIOS
+#     module load nco
+#     R_IN=$(ccc_home -u igcmg --cccwork)/IGCM
+#     TMPDIR=${SCRATCHDIR}/TMP
+#     SUBMIT_DIR=${BRIDGE_MSUB_PWD:-${SUBMIT_DIR}}
+#     ;;
+#     ( spip )
+#     R_IN=${HOME}/Scratch/IGCM
+#     TMPDIR=${HOME}/Scratch/TMP
+#     SUBMIT_DIR=$(pwd)
+#     MPIRUN=/opt/local/bin/mpirun-openmpi-gcc49 -n 2
+#     ;; 
+#     ( * ) exit -1 ;;
+# esac
 
 ##
@@ -185 +185 @@
 
 ncks --fl_fmt=${FL_FMT3} --history -C --variable OceMask         ${ATM}_grid_maskFrom_${OCE}_${FL_FMT3}.nc OceMask_tmp.nc
+ncwa --history --overwrite --average time_counter OceMask_tmp.nc OceMask_tmp2.nc # Suppress time dimension
 ncatted --history \
+        --attribute cell_methods,OceMask,d,c,       \
         --attribute coordinates,OceMask,d,,         \
         --attribute online_operation,OceMask,d,,    \
@@ -191 +193 @@
         --attribute missing_value,OceMask,d,,       \
         OceMask_tmp.nc
-ncap2 --fl_fmt=${FL_FMT3} --history --append --script "OceMask=int(1-OceMask)"            OceMask_tmp.nc masks_${CplModel}.nc # For OASIS, ocean=0, land=1
+ncap2 --fl_fmt=${FL_FMT3} --history --append --script "OceMask=int(1-OceMask)" OceMask_tmp2.nc masks_${CplModel}.nc # For OASIS, ocean=0, land=1
 rm OceMask_tmp.nc
 ncatted --history \
@@ -216 +218 @@
 
 #
-echo ${Titre}"DYNAMICO grids : 1) Set correct order of dimensions, 2) Add a dimension"${Norm}
-# ----------------------------------------------------------------------------
 if [[ ${atm} = dynamico* || ${atm} = ico* ]] ; then
+    echo ${Titre}"DYNAMICO grids : 1) Set correct order of dimensions, 2) Add a dimension"${Norm}
+    # -------------------------------------------------------------------------------------------
     mv grids_${CplModel}.nc grids_${CplModel}_tmp.nc
     ncpdq --permute nvertex,ycell,xcell  grids_${CplModel}_tmp.nc grids_${CplModel}.nc
-    
+
     mv grids_${CplModel}.nc grids_${CplModel}_tmp.nc
     cat <<EOF >add_dim.nco
@@ -232 +234 @@
 //
 EOF
-    
+   
     ncap2 --fl_fmt=${FL_FMT3} --history --append --script-file add_dim.nco  grids_${CplModel}_tmp.nc grids_${CplModel}.nc
     
     mv masks_${CplModel}.nc masks_${CplModel}_tmp.nc
-    ncap2 --fl_fmt=${FL_FMT3} --history --append --script 'defdim("xcell",1) ; OceMask[time_counter,ycell,xcell]=OceMask(:,:)'  masks_${CplModel}_tmp.nc masks_${CplModel}.nc
+    ncap2 --fl_fmt=${FL_FMT3} --history --append --script 'defdim("xcell",1) '   masks_${CplModel}_tmp.nc masks_${CplModel}.nc
+    mv masks_${CplModel}.nc masks_${CplModel}_tmp.nc
+    ncap2 --fl_fmt=${FL_FMT3} --history --append --script 'OceMask[ycell,xcell]=OceMask'  masks_${CplModel}_tmp.nc masks_${CplModel}.nc
 
     mv areas_${CplModel}.nc areas_${CplModel}_tmp.nc
-    ncap2 --fl_fmt=${FL_FMT3} --history --append --script 'defdim("xcell",1) ; aire[ycell,xcell]=aire(:) ; '                    areas_${CplModel}_tmp.nc areas_${CplModel}.nc
+    ncap2 --fl_fmt=${FL_FMT3} --history --append --script 'defdim("xcell",1) ; aire[ycell,xcell]=aire ; ' areas_${CplModel}_tmp.nc areas_${CplModel}.nc
 
     rm -f grids_${CplModel}_tmp.nc areas_${CplModel}_tmp.nc masks_${CplModel}_tmp.nc
 fi
-
+ 
 #
 echo ${Titre}"Some changes specific to LMDZ lon/lat grid"${Norm}
@@ -292 +296 @@
 
 #
+ 
 echo ${Titre}"Generates grid o${atm}, 'o' meaning 'one'"${Norm}
 # same as t${atm} grid, with surfaces set to 1
@@ -308 +313 @@
 rm grids_${CplModel}_tmp.nc areas_${CplModel}_tmp.nc masks_${CplModel}_tmp.nc
 
+
+echo ${Titre}"Generates grid o${oce}, 'o' meaning 'one'"${Norm}
+# same as t${oce} grid, with surfaces set to 1
+# and mask to 0 (ocean everywhere, to compute integral over the whole grid))
+# This grid is used when field are quantities instead of fluxes (i.e river flow)
+# --------------------------------------------------------------------------------------------------------
+mv grids_${CplModel}.nc grids_${CplModel}_tmp.nc
+ncap2 --fl_fmt=${FL_FMT3} --history --script "nav_lon_grid_O=nav_lon_grid_T; nav_lat_grid_O=nav_lat_grid_T; bounds_lon_grid_O=bounds_lon_grid_T; bounds_lat_grid_O=bounds_o${oce}_lat=bounds_lat_grid_T; " grids_${CplModel}_tmp.nc grids_${CplModel}.nc
+
+mv areas_${CplModel}.nc areas_${CplModel}_tmp.nc
+ncap2 --fl_fmt=${FL_FMT3} --history --script "area_grid_O=area_grid_T*0.0d+1.0d;"   areas_${CplModel}_tmp.nc areas_${CplModel}.nc
+
+ncks -C --history --overwrite -v maskutil_T ${OCE}_coordinates_mask_${FL_FMT3}.nc maskutil_T_${FL_FMT3}.nc
+ncap2 --fl_fmt=${FL_FMT3} --history --append --script "mask_O=maskutil_T; " maskutil_T_${FL_FMT3}.nc masks_${CplModel}.nc
+
+rm grids_${CplModel}_tmp.nc areas_${CplModel}_tmp.nc 
+
+
+echo ${Titre}"Creates OCEAN C grid : redundant points removed to compute proper integrals"${Norm}
+# --------------------------------------------------------------------------------------------------------
+
+mv grids_${CplModel}.nc grids_${CplModel}_tmp.nc
+ncap2 --fl_fmt=${FL_FMT3} --history --script "nav_lon_grid_C=nav_lon_grid_T; nav_lat_grid_C=nav_lat_grid_T; bounds_lon_grid_C=bounds_lon_grid_T; bounds_lat_grid_C=bounds_o${oce}_lat=bounds_lat_grid_T; " grids_${CplModel}_tmp.nc grids_${CplModel}.nc
+
+mv areas_${CplModel}.nc areas_${CplModel}_tmp.nc
+ncap2 --fl_fmt=${FL_FMT3} --history --script "area_grid_C=area_grid_T ; "   areas_${CplModel}_tmp.nc areas_${CplModel}.nc
+
+#mv masks_${CplModel}.nc masks_${CplModel}_tmp.nc
+ncap2 --fl_fmt=${FL_FMT3} --history --append --script "mask_C=maskutil_T; " maskutil_T_${FL_FMT3}.nc masks_${CplModel}.nc
+
+rm grids_${CplModel}_tmp.nc areas_${CplModel}_tmp.nc maskutil_T_${FL_FMT3}.nc
+
 #
 echo ${Titre}"Final renaming"${Norm}
 # ----------------------------------------------------------------------------
-for OCEGRID in T U V
+for OCEGRID in T U V O C
 do
     ocegrid=${OCEGRID~} # To lowercase
@@ -342 +379 @@
 ncrename --history --variable bounds_o${atm}_lat,o${atm}.cla  grids_${CplModel}.nc
 
-for atmgrid in t o ; do
+for ATMGRID in T O ; do
+    atmgrid=${ATMGRID~} # To lowercase
     ncatted --history \
             --attribute bounds,${atmgrid}${atm}.lon,c,c,"${atmgrid}${atm}.clo"   \
@@ -414 +452 @@
 
 ## ===========================================================================
-echo ${Titre}"exit from CreateOasisGrids "${Norm}
+echo ${Titre}"Exit from CreateOasisGrids "${Norm}
 
 ## ===========================================================================