Changeset 6064 for TOOLS/MOSAIX/CreateOasisGrids.bash

Timestamp:

02/24/22 09:58:49 (2 years ago)

Author:

omamce

Message:

O.M. : changes in MOSAIX

• Correct ocean mask to remove periodic duplicated points, to have full conservation of run-off
• Suppress creation of corc mask, which is not used
• Put SVN keywords in make_mosaic
• Update nemo.py with somùe additionnal utilities
• Adapt RunOffWeights?.py to new nemo module
• Adapt CalvingWeights?.py to new ORCA configurations

File:

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

TOOLS/MOSAIX/CreateOasisGrids.bash

@@ -51 +51 @@
 if [[ $(hostname) = irene* ]]    ; then arch=irene ; center=tgcc ; fi
 if [[ $(hostname) = lsce3005* ]] ; then arch=spip  ; center=spip ; fi
+if [[ $(hostname) = lsce5138* ]] ; then arch=spip  ; center=spip ; fi
 
 PROGRAM=$(basename ${0})
@@ -152 +153 @@
     ncks -C --history --append --variable area_grid_${OCEGRID} ${OCE}_coordinates_mask.nc areas_${CplModel}.nc
     # Inverts mask values and switch to integer
-    ncks --history -C --variable mask_${OCEGRID} ${OCE}_coordinates_mask.nc mask_${OCEGRID}_tmp.nc
-    ncatted --history \
-            --attribute coordinates,mask_${OCEGRID},d,,         \
-            --attribute online_operation,mask_${OCEGRID},d,,    \
-            --attribute cell_measures,mask_${OCEGRID},d,,       \
+    ncks --history -C --variable maskutil_${OCEGRID} ${OCE}_coordinates_mask.nc mask_${OCEGRID}_tmp.nc
+    ncatted --history \
+            --attribute coordinates,maskutil_${OCEGRID},d,,         \
+            --attribute online_operation,maskutil_${OCEGRID},d,,    \
+            --attribute cell_measures,maskutil_${OCEGRID},d,,       \
     mask_${OCEGRID}_tmp.nc
     
-    ncap2 --history --append --script "mask_${OCEGRID}=int(1-mask_${OCEGRID});"    mask_${OCEGRID}_tmp.nc masks_${CplModel}.nc
+    ncap2 --history --append --script "maskutil_${OCEGRID}=int(1-maskutil_${OCEGRID});"    mask_${OCEGRID}_tmp.nc masks_${CplModel}.nc
     rm mask_${OCEGRID}_tmp.nc
     ncatted --history \
-            --attribute long_name,mask_${OCEGRID},o,c,"Land-sea mask" \
-            --attribute units,mask_${OCEGRID},o,c,"Land:1, Ocean:0"    masks_${CplModel}.nc
+            --attribute long_name,maskutil_${OCEGRID},o,c,"Land-sea mask" \
+            --attribute units,maskutil_${OCEGRID},o,c,"Land:1, Ocean:0"    masks_${CplModel}.nc
     # Change order of dimensions
     mv grids_${CplModel}.nc grids_${CplModel}_tmp.nc
@@ -327 +328 @@
 
 
-echo ${Titre}"Generates grid c${atm}, 'o' meaning 'one'"${Norm}
-# same as t${atm} grid, with surfaces set to grid area
-# and mask to 0 (ocean everywhere, to compute integral over the whole grid))
-# ----------------------------------------------------------------------------
-mv grids_${CplModel}.nc grids_${CplModel}_tmp.nc
-ncap2 --history --script "c${atm}_lon=alon ; c${atm}_lat=alat ; bounds_c${atm}_lon=bounds_lon ; bounds_c${atm}_lat=bounds_lat ; " grids_${CplModel}_tmp.nc grids_${CplModel}.nc
-
-mv areas_${CplModel}.nc areas_${CplModel}_tmp.nc
-ncap2 --history --script "c${atm}_aire=aire"   areas_${CplModel}_tmp.nc areas_${CplModel}.nc
-
-mv masks_${CplModel}.nc masks_${CplModel}_tmp.nc
-ncap2 --history --script "c${atm}_mask=int(Oce2AtmMask)*0+0 ; " masks_${CplModel}_tmp.nc masks_${CplModel}.nc
-
-rm grids_${CplModel}_tmp.nc areas_${CplModel}_tmp.nc masks_${CplModel}_tmp.nc
+# echo ${Titre}"Generates grid c${atm}, 'o' meaning 'one'"${Norm}
+# # same as t${atm} grid, with surfaces set to grid area
+# # and mask to 0 (ocean everywhere, to compute integral over the whole grid))
+# # ----------------------------------------------------------------------------
+# mv grids_${CplModel}.nc grids_${CplModel}_tmp.nc
+# ncap2 --history --script "c${atm}_lon=alon ; c${atm}_lat=alat ; bounds_c${atm}_lon=bounds_lon ; bounds_c${atm}_lat=bounds_lat ; " grids_${CplModel}_tmp.nc grids_${CplModel}.nc
+
+# mv areas_${CplModel}.nc areas_${CplModel}_tmp.nc
+# ncap2 --history --script "c${atm}_aire=aire"   areas_${CplModel}_tmp.nc areas_${CplModel}.nc
+
+# mv masks_${CplModel}.nc masks_${CplModel}_tmp.nc
+# ncap2 --history --script "c${atm}_mask=int(Oce2AtmMask)*0+0 ; " masks_${CplModel}_tmp.nc masks_${CplModel}.nc
+
+# 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)
+# This grid is used when field are quantities instead of fluxes
 # --------------------------------------------------------------------------------------------------------
 mv grids_${CplModel}.nc grids_${CplModel}_tmp.nc
@@ -354 +355 @@
 
 ncks -C --history --overwrite -v maskutil_T ${OCE}_coordinates_mask.nc maskutil_T.nc
-ncap2 --history --append --script "mask_O=maskutil_T; " maskutil_T.nc masks_${CplModel}.nc
+ncap2 --history --append --script "maskutil_O=maskutil_T; " maskutil_T.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 --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 --history --script "area_grid_C=area_grid_T ; "   areas_${CplModel}_tmp.nc areas_${CplModel}.nc
-
-ncap2 --history --append --script "mask_C=maskutil_T; " maskutil_T.nc masks_${CplModel}.nc
-
-rm grids_${CplModel}_tmp.nc areas_${CplModel}_tmp.nc maskutil_T.nc
+# echo ${Titre}"Creates OCEAN C grid : redundant points removed to compute proper integrals"${Norm}
+# # --------------------------------------------------------------------------------------------------------
+
+# mv grids_${CplModel}.nc grids_${CplModel}_tmp.nc
+# ncap2 --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 --history --script "area_grid_C=area_grid_T ; "   areas_${CplModel}_tmp.nc areas_${CplModel}.nc
+
+# ncap2 --history --append --script "mask_C=maskutil_T; " maskutil_T.nc masks_${CplModel}.nc
+
+# rm grids_${CplModel}_tmp.nc areas_${CplModel}_tmp.nc maskutil_T.nc
 
 #
 echo ${Titre}"Final renaming"${Norm}
 # ----------------------------------------------------------------------------
-for OCEGRID in T U V O C
+for OCEGRID in T U V O
 do
     ocegrid=${OCEGRID~} # To lowercase
@@ -382 +383 @@
     ncrename --history --variable bounds_lon_grid_${OCEGRID},${ocegrid}${oce}.clo  grids_${CplModel}.nc
     ncrename --history --variable bounds_lat_grid_${OCEGRID},${ocegrid}${oce}.cla  grids_${CplModel}.nc
-    ncrename --history --variable mask_${OCEGRID},${ocegrid}${oce}.msk             masks_${CplModel}.nc
+    ncrename --history --variable maskutil_${OCEGRID},${ocegrid}${oce}.msk         masks_${CplModel}.nc
     ncrename --history --variable area_grid_${OCEGRID},${ocegrid}${oce}.srf        areas_${CplModel}.nc
     
@@ -399 +400 @@
 ncrename --history --variable bounds_lat,t${atm}.cla          grids_${CplModel}.nc
 
-for ATMGRID in O C ; do
+for ATMGRID in O ; do
     atmgrid=${ATMGRID~} # To lowercase
     ncrename --history --variable ${atmgrid}${atm}_lon,${atmgrid}${atm}.lon         grids_${CplModel}.nc
@@ -409 +410 @@
 done
     
-for ATMGRID in T O C ; do
+for ATMGRID in T O ; do
     atmgrid=${ATMGRID~} # To lowercase
     ncatted --history \