Changeset 4186 for TOOLS/MOSAIX/CalvingWeights.py

Timestamp:

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

Author:

omamce

Message:

O.M. :

Runoff weight for both ICO and LMDZ
Correct corc grid

File:

• TOOLS/MOSAIX/CalvingWeights.py (modified) (4 diffs)

TOOLS/MOSAIX/CalvingWeights.py

@@ -43 +43 @@
 parser.add_argument ('--repartition_var' , help='Variable name for iceshelf'      , type=str, default=None)
 parser.add_argument ('--output'  , help='output rmp file name', default='rmp_tlmd_to_torc_calving_64bit.nc' )
-parser.add_argument ('--fmt'     , help='NetCDF file format, using nco syntax', default='64bits', choices=['classic', 'netcdf3', '64bit', '64bit_data', '64bit_data', 'netcdf4', 'netcdf4_classsic'] )
+parser.add_argument ('--fmt'     , help='NetCDF file format, using nco syntax', default='64bit', choices=['classic', 'netcdf3', '64bit', '64bit_data', '64bit_data', 'netcdf4', 'netcdf4_classsic'] )
 
 # Parse command line
@@ -308 +308 @@
 f_calving.Program         = "Generated by " + sys.argv[0] + " with flags " + str(sys.argv[1:])
 f_calving.repartitionType = myargs.type
-if myargs.type in ['iceberg', 'iceshelf' ] :
+if myargs.type in [ 'iceberg', 'iceshelf' ] :
     f_calving.repartitionFile = myargs.repartition_file
     f_calving.repartitionVar  = repartitionVar
@@ -333 +333 @@
 f_calving.SVN_HeadURL     = "$HeadURL$"
 
-d_nb_zone          = f_calving.createDimension ('nb_zone'         ,   nb_zone )
-d_num_links        = f_calving.createDimension ('num_links'       , num_links )
-d_num_wgts         = f_calving.createDimension ('num_wgts'        ,         1 )
+d_nb_zone   = f_calving.createDimension ('nb_zone'         ,   nb_zone )
+d_num_links = f_calving.createDimension ('num_links'       , num_links )
+d_num_wgts  = f_calving.createDimension ('num_wgts'        ,         1 )
 
 d_src_grid_size    = f_calving.createDimension ('src_grid_size'   , src_grid_size )
 d_src_grid_corners = f_calving.createDimension ('src_grid_corners', src_clo.shape[0]  )
-D_src_grid_rank    = f_calving.createDimension ('src_grid_rank'   ,        2  )
+d_src_grid_rank    = f_calving.createDimension ('src_grid_rank'   ,        2  )
 
 d_dst_grid_size    = f_calving.createDimension ('dst_grid_size'   , dst_grid_size )
@@ -398 +398 @@
 v_dst_repartition = f_calving.createVariable ( 'dst_repartition' , 'f8', ('nb_zone', 'dst_grid_size',)  )
 
+v_dst_southLimit = f_calving.createVariable ('dst_southLimit', 'f4', ('nb_zone',) )
+v_dst_northLimit = f_calving.createVariable ('dst_northLimit', 'f4', ('nb_zone',) )
+v_dst_southLimit[:] = np.min(limit_lat, axis=(1,) ) 
+v_dst_northLimit[:] = np.max(limit_lat, axis=(1,) ) 
+
 v_dst_grid_dims      [:] = ( dst_jpi, dst_jpi ) 
 v_dst_grid_center_lon[:] = dst_lon[:].ravel()