Changeset 2238 in lmdz_wrf for trunk/tools/read_grib.py

Timestamp:

Nov 17, 2018, 2:56:47 PM (6 years ago)

Author:

lfita

Message:

Working version with arguments of `read_grib.py'

Fixing error message in `datetimeStr_conversion'

File:

• trunk/tools/read_grib.py (modified) (3 diffs)

trunk/tools/read_grib.py

@@ -1 +1 @@
 # FROM: https://pypi.org/project/pupygrib/
+from optparse import OptionParser
 import numpy as np
 import pupygrib
@@ -8 +9 @@
 
 errormsg = 'ERROR -- error -- ERROR -- error'
-
-ifile = 'ERAI_pl200503_130.grib'
-
-ofile = 'ERAI_pl_130.nc'
 
 def grib128_ncname(varid):
@@ -53 +50 @@
     return [varn,Lvarn,units,cfvalues[0], cfvalues[1], cfvalues[4].replace('|', ' ')]
 
-timeu = 'minutes since 1949-12-01 00:00:00'
+parser = OptionParser()
+parser.add_option("-g", "--grib_file", dest="gribfile", help="GRIB file to use", 
+  metavar="FILE")
+parser.add_option("-o", "--output_file", dest="ncfile", help="name of NC file to generate", 
+  metavar="FILE")
+parser.add_option("-t", "--timeunits", dest="timeunits", 
+  help="CF reference of time values [Tunits] since [DateRef] ('!', for spaces minutes!since!1949-12-01!00:00:00 as default)", 
+  metavar="STRING")
+
+(opts, args) = parser.parse_args()
 
 #######    #######
 ## MAIN
     #######
+
+if opts.gribfile is None:
+    print errormsg
+    print "  no grib input file name provided !!"
+    print "     you should provide one file name as -g [FileName]"
+    quit(-1)
+
+if not os.path.isfile(opts.gribfile):
+    print errormsg
+    print "  no grib file '" + opts.gribfile + "' found !!"
+    quit(-1)
+else:
+    ifile = opts.gribfile
+
+if opts.ncfile is None:
+    print errormsg
+    print "  no netcdf output file name provided !!"
+    print "     you should provide one file name as -o [FileName]"
+    quit(-1)
+else:
+    ofile = opts.ncfile
+
+if opts.timeunits is None:
+    print errormsg
+    print "  no time units provided !!"
+    print "     using the default ones as 'minutes since 1949-12-01 00:00:00'"
+    quit(-1)
+else:
+  timeu = opts.timeunits.replace('!', ' ')
 
 with open(ifile, 'rb') as stream: