Changeset 176

Timestamp:

03/22/12 15:17:20 (12 years ago)

Author:

pinsard

Message:

parametrization of era-i new utility

File:

• trunk/src/get_erai.py (modified) (7 diffs)

trunk/src/get_erai.py

@@ -2 +2 @@
 # -*- coding: iso-8859-1 -*-
 
+__docformat__ = "restructuredtext en"
+__version__ = '$Id$'
+
 """
 
@@ -7 +10 @@
 ===========
 
-get ERA-Interim reference files
+get ERA-Interim reference files between two dates
 
 SYNOPSIS
@@ -14 +17 @@
 ::
 
- $ python get_erai.py
-
+ $ get_erai.py --begin yyyymmdd --end yyyymmdd
+
+EXAMPLES
+========
+
+::
+
+ $ get_erai.py --verbose --begin 20010101 --end 20010330
 
 SEE ALSO
@@ -22 +31 @@
 :ref:`ecmwf`
 
+for interactif retrieval
+http://data-portal.ecmwf.int/data/d/interim_daily/
+
+http://data-portal.ecmwf.int/data/d/interim_full_daily
+
+For python script syntax
+http://data-portal.ecmwf.int/data/d/token/interim_daily/
+
+"Please note that these batch scripts are available for download in grib format only. " ...
+
+MARS language
+http://www.ecmwf.int/publications/manuals/mars/guide/index.html
+
 TIPS
 ====
 
 To convet a GRIB file to a NetCDF file ::
- 
+
   $ cdo -f nc copy data.grib data.nc
 
@@ -32 +54 @@
 ====
 
-make it work for all needed variables 
+make it work for all needed variables
 
 produce netcdf files : now grib dispite name of the file target
 
 comment
+
+nb : le champ str (wind stress pourra/devra être récupéré pour faire des comparaisons (cf. 20c3m_erai_str_TROP_1989_2009.nc)
+
+nb :Mean sea level pressure en a-t-on vraiment besoin ? jusqu'à présent on ne sait pas si pkb l'a utilisé.
 
 EVOLUTIONS
@@ -45 +71 @@
 $URL$
 
+- fplod 20120314
+
+  * add parameters
+
+- fplod 20120306
+
+  * add all param for analyse (it might be better to do a loop to produce
+    one file/field
+  * add retrieve for forecast (only one parameter here)
+
 - fplod 20111130T153500Z cratos (Linux)
 
@@ -52 +88 @@
 """
 
+import os
+from optparse import OptionParser
+import time
+import datetime
+from dateutil.parser import *
+
 from ecmwf import ECMWFDataServer
 
-server = ECMWFDataServer(
-       'http://data-portal.ecmwf.int/data/d/dataserver/',
-              'd558d095ce964389eeda5e908e9462f3',
-             'francoise.pinsard@locean-ipsl.upmc.fr'
-         )
-
-server.retrieve({
-    'dataset' : "interim_full_daily",
-    'date'    : "19790101/to/19790110",
-    'time'    : "00:00:00/06:00:00/12:00:00/18:00:00",
-    'grid'    : "1/1",
-    'step'    : "0",
-    'levtype' : "sfc",
-    'type'    : "an",
-    'param'   : "167.128",
-    'area'    : "60/-120/30/-60",
-    'target'  : "t2m.nc"
-        })
+def get_option_parser ():
+    """parse CLI arguments
+
+    :returns: parser
+    :rtype: :class:`optparse.OptionParser`
+    """
+
+    parser = OptionParser('%prog [--verbose] --begin yyyymmdd --end yyymmdd ')
+
+    parser = OptionParser(conflict_handler="resolve")
+
+    parser.add_option ('-V', '--verbose', help='produce verbose output',
+          action='store_true',
+          default=False,
+          dest='is_verbose')
+
+    parser.add_option('--version',
+         action='store_true',
+         dest='version',
+         help='show version number and exit')
+
+    parser.add_option('--begin',
+           help='date (yyyymmdd) of the first date',
+           dest='yyyymmddb',
+           action='store')
+
+    parser.add_option('--end',
+           help='date (yyyymmdd) of the last date',
+           dest='yyyymmdde',
+           action='store')
+
+    return parser
+
+def get_erai():
+    """main
+    """
+    try:
+       parser = get_option_parser ()
+       (fromcli, args) = parser.parse_args()
+    except IOError, msg:
+       parser.error(str(msg))
+
+    yyyymmddb = parse(fromcli.yyyymmddb, ignoretz=True)
+    yyyymmdde = parse(fromcli.yyyymmdde, ignoretz=True)
+
+    # ++ check date validity and interval validity
+    if fromcli.is_verbose == True:
+       print ('yyyymmddb %s' % (yyyymmddb))
+       print ('yyyymmdde %s' % (yyyymmdde))
+
+    # build the value of date parameter yyyymmdd/to/yyyymmdd
+    date_param = fromcli.yyyymmddb + '/to/' + fromcli.yyyymmdde
+
+    if fromcli.is_verbose == True:
+       print ('date for erai %s' % (date_param))
+
+    server = ECMWFDataServer(
+           'http://data-portal.ecmwf.int/data/d/dataserver/',
+                  'd558d095ce964389eeda5e908e9462f3',
+                 'Francoise.Pinsard@locean-ipsl.upmc.fr'
+             )
+
+    # Pour les analyses de
+    #        10 metre U wind component
+    #        10 metre V wind component
+    #        2 metre temperature
+    #        2 metre dewpoint temperature
+    #        Sea surface temperature
+    #        Mean sea level pressure
+    param_param = []
+    param_param.append("34.128") #Sea surface temperature sstk
+    #param_param.append("151.128") Mean sea level pressure
+    #param_param.append("165.128") 10 metre U wind component u10
+    #param_param.append("166.128") 10 metre V wind component v10
+    #param_param.append("167.128") 2 metre temperature t2
+    #param_param.append("168.128") 2 metre dewpoint temperature d2
+
+    for one_param in param_param:
+
+        if fromcli.is_verbose == True:
+            print('recup de %s ' %  one_param)
+
+        server.retrieve({
+            'dataset' : "interim_full_daily",
+            'date'    : date_param,
+            'time'    : "00:00:00/06:00:00/12:00:00/18:00:00",
+            'grid'    : "1/1",
+            'step'    : "0",
+            'levtype' : "sfc",
+            'type'    : "an",
+            'param'   : one_param,
+            'area'    : "60/-120/30/-60",
+            'target'  : "toto.grib"
+            })
+
+    # Pour la "Surface thermal radiation downwards":
+    param_param = []
+    param_param.append("175.128")  # Surface thermal radiation downwards
+    for one_param in param_param:
+
+       if fromcli.is_verbose == True:
+           print('recup de %s'  % one_param)
+
+       server.retrieve({
+           'dataset' : "interim_full_daily",
+           'date'    : date_param,
+           'time'    : "00:00:00/12:00:00",
+           'grid'    : "1/1",
+           'step'    : "12",
+           'levtype' : "sfc",
+           'type'    : "fc",
+           'param'   : param_param,
+           'area'    : "60/-120/30/-60",
+           'target'  : "fc.grib"
+             })
+
+# Run main, if called from the command line
+if __name__ == '__main__':
+   get_erai()
+