Changeset 3901 for TOOLS/MOSAIX/ComputeNemoCoast.py

Timestamp:

06/07/18 15:59:11 (6 years ago)

Author:

omamce

Message:

O.M. : add weight for coastal points only, to prepare computation

of weights for river run-off

File:

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

TOOLS/MOSAIX/ComputeNemoCoast.py

@@ -27 +27 @@
 
 def usage () :
-    texte = """%(prog)s usage :
-python %(prog)s [-d] [-i <orca grid file>] [-n <perio>]
+    texte = """Compute a mask for coastal point of NEMO domain
+
+    %(prog)s usage :
+python %(prog)s [-h|--help] [-d|--debug] [-i <orca grid file>] [-n <perio>]
  -d         | --debug         : debug
  -i <file>  | --input=<file>  : input file  (default: none)
@@ -49 +51 @@
 ## Command line options
 try:
-    myopts, myargs = getopt.getopt ( sys.argv[1:], 'i:n:h', [ 'input=', 'nperio=', 'debug=', '--help' ] )
+    myopts, myargs = getopt.getopt ( sys.argv[1:], 'i:n:h', [ 'input=', 'nperio=', 'debug=', 'help' ] )
 except getopt.GetoptError as cmdle :
     print ( "Command line error : "+str(cmdle)+"\n" )
@@ -66 +68 @@
     sys.exit(-1)
 
-print ("Input file  :" + GridFile)   
+print ("Input file : " + GridFile)   
 
 ## Open grid file
@@ -115 +117 @@
 print ("Maximum number of neighbors : "+str(CoastCrit) )
 
-
-Temp = nemo.lbc ( Temp,  nperio=4, cd_type='T' )
+Temp = nemo.lbc ( Temp,  nperio=nperio, cd_type='T' )
 
 OceCoastal = np.where (OceMask == 1, True, False) * np.where (Temp < CoastCrit, True, False)