Changeset 5915 for TOOLS/MOSAIX/RunoffWeights.py

Timestamp:

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

Author:

snguyen

Message:

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

File:

• TOOLS/MOSAIX/RunoffWeights.py (modified) (3 diffs)

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)