Changeset 3901 for TOOLS/MOSAIX/ComputeNemoCoast.py
- Timestamp:
- 06/07/18 15:59:11 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TOOLS/MOSAIX/ComputeNemoCoast.py
r3725 r3901 27 27 28 28 def usage () : 29 texte = """%(prog)s usage : 30 python %(prog)s [-d] [-i <orca grid file>] [-n <perio>] 29 texte = """Compute a mask for coastal point of NEMO domain 30 31 %(prog)s usage : 32 python %(prog)s [-h|--help] [-d|--debug] [-i <orca grid file>] [-n <perio>] 31 33 -d | --debug : debug 32 34 -i <file> | --input=<file> : input file (default: none) … … 49 51 ## Command line options 50 52 try: 51 myopts, myargs = getopt.getopt ( sys.argv[1:], 'i:n:h', [ 'input=', 'nperio=', 'debug=', ' --help' ] )53 myopts, myargs = getopt.getopt ( sys.argv[1:], 'i:n:h', [ 'input=', 'nperio=', 'debug=', 'help' ] ) 52 54 except getopt.GetoptError as cmdle : 53 55 print ( "Command line error : "+str(cmdle)+"\n" ) … … 66 68 sys.exit(-1) 67 69 68 print ("Input file :" + GridFile)70 print ("Input file : " + GridFile) 69 71 70 72 ## Open grid file … … 115 117 print ("Maximum number of neighbors : "+str(CoastCrit) ) 116 118 117 118 Temp = nemo.lbc ( Temp, nperio=4, cd_type='T' ) 119 Temp = nemo.lbc ( Temp, nperio=nperio, cd_type='T' ) 119 120 120 121 OceCoastal = np.where (OceMask == 1, True, False) * np.where (Temp < CoastCrit, True, False)
Note: See TracChangeset
for help on using the changeset viewer.