Changeset 5915 for TOOLS

Timestamp:

09/13/21 10:24:46 (3 years ago)

Author:

snguyen

Message:

Ajout de l'option --ocePerio à RunoffWeights?.py et CalvingWeights?.py

Location:

TOOLS/MOSAIX

Files:

• CalvingWeights.py (modified) (2 diffs)
• RunoffWeights.py (modified) (3 diffs)

TOOLS/MOSAIX/CalvingWeights.py

@@ -44 +44 @@
 parser.add_argument ('--output'  , help='output rmp file name', default='rmp_tlmd_to_torc_calving.nc' )
 parser.add_argument ('--fmt'     , help='NetCDF file format, using nco syntax', default='netcdf4', choices=['classic', 'netcdf3', '64bit', '64bit_data', '64bit_data', 'netcdf4', 'netcdf4_classsic'] )
+parser.add_argument ('--ocePerio'   , help='periodicity of ocean grid', type=int, default=0 )
 
 # Parse command line
@@ -110 +111 @@
 print ('dst_mskutil : ' + str(np.sum(dst_mskutil)))
 
+# Ocean grid periodicity
+nperio_dst=myargs.ocePerio
+
 # Periodicity masking for NEMO
-if dst_Name == 'ORCA2.3'    : nperio_dst = 4
-if dst_Name == 'eORCA1.2'   : nperio_dst = 6
-if dst_Name == 'ORCA025'    : nperio_dst = 6
-if dst_Name == 'eORCA025'   : nperio_dst = 6
-if dst_Name == 'eORCA025.1' : nperio_dst = 6
+if nperio_dst == 0 :
+    if dst_Name == 'ORCA2.3'    : nperio_dst = 4
+    if dst_Name == 'eORCA1.2'   : nperio_dst = 6
+    if dst_Name == 'ORCA025'    : nperio_dst = 6
+    if dst_Name == 'eORCA025'   : nperio_dst = 6
+    if dst_Name == 'eORCA025.1' : nperio_dst = 6
 print ('nperio_dst: ' + str(nperio_dst) )
 dst_mskutil = nemo.lbc_mask (dst_mskutil, nperio=nperio_dst, cd_type='T' )

TOOLS/MOSAIX/RunoffWeights.py

@@ -68 +68 @@
 parser.add_argument ('--output', help='output rmp file name', default='rmp_tlmd_to_torc_runoff.nc' )
 parser.add_argument ('--fmt'   , help='NetCDF file format, using nco syntax', default='netcdf4', choices=['classic', 'netcdf3', '64bit', '64bit_data', '64bit_data', 'netcdf4', 'netcdf4_classsic'] )
+parser.add_argument ('--ocePerio'   , help='periodicity of ocean grid', type=int, default=0 )
 
 # Parse command line
@@ -101 +102 @@
 print ('atmQuantity : ' + str (myargs.atmQuantity) )
 print ('oceQuantity : ' + str (myargs.oceQuantity) )
-    
+
+# Ocean grid periodicity
+oce_perio=myargs.ocePerio
+
 ### Read coordinates of all models
 ###
@@ -141 +145 @@
 oce_grid_imask      = 1-maskFile['torc.msk'][:].ravel()
 oce_grid_dims       = gridFile['torc.lat'][:].shape
-if oce_jpi ==  182 : oce_perio = 4 # ORCA 2
-if oce_jpi ==  362 : oce_perio = 6 # ORCA 1
-if oce_jpi == 1442 : oce_perio = 6 # ORCA 025
+if oce_perio == 0 :
+    if oce_jpi ==  182 : oce_perio = 4 # ORCA 2
+    if oce_jpi ==  362 : oce_perio = 6 # ORCA 1
+    if oce_jpi == 1442 : oce_perio = 6 # ORCA 025
+print(f"oce_perio = {oce_perio}")
 oce_grid_pmask = nemo.lbc_mask (np.reshape(oce_grid_imask, (oce_jpj,oce_jpi)), 'T', oce_perio).ravel()
 oce_address = np.arange(oce_jpj*oce_jpi)