Changeset 558 for trunk

Timestamp:

Mar 4, 2012, 9:24:59 PM (13 years ago)

Author:

aslmd

Message:

UTIL PYTHON: moved corrected ECMWF stuff for simple use in the main PYTHON directory. corrected planetoplot so that output files from ECMWF can be used. added a more general treatment of lat/lon prescription.

Location:

trunk/UTIL/PYTHON

Files:

• create_ecmwf.py (added)
• ecmwf.py (copied) (copied from trunk/UTIL/PYTHON/examples/mylib/ecmwf.py)
• examples/map.py (copied) (copied from trunk/UTIL/PYTHON/examples/mylib/map.py)
• examples/map_ecmwf.py (copied) (copied from trunk/UTIL/PYTHON/examples/mylib/map_ecmwf.py)
• examples/mylib (deleted)
• myecmwf.py (copied) (copied from trunk/UTIL/PYTHON/examples/mylib/myecmwf.py) (2 diffs)
• myplot.py (modified) (3 diffs)
• planetoplot.py (modified) (1 diff)

trunk/UTIL/PYTHON/myecmwf.py

@@ -7 +7 @@
         elif var == "130": name="Atmospheric Temperature (K)"           # ["130", "an", "0", "pl" ]
         elif var == "167": name="2m Atmospheric Temperature (K)"        # ["167", "fc", "3", "sfc"]
-        elif var == "78": name="Total column liquid water (kg/kg)"
+        elif var == "78":  name="Total column liquid water (kg/kg)"
         elif var == "137": name="Total column water vapour (kg/m2)"
+        else:              name=""
         return name
 
@@ -22 +23 @@
         return varchar
 
-def get_ecmwf (var, dataset, wlat, wlon, lev, date, tim):
+def get_ecmwf (var, fieldtype, wlat, wlon, lev, date, tim):
         from ecmwf      import ECMWFDataServer
         from os         import system
         from netCDF4    import Dataset
         gbfile = 'output.grib'
-        ncfile = 'output.nc'
+        #ncfile = str(date[2]) + str(date[1]) + str(date[0]) + str(date[5]) + str(date[4]) + str(date[3])+'.nc'
+        ncfile = 'output.nc'
+        ######
+        if fieldtype == "3d":    
+            dataset = ["an","pl","0"]
+        elif fieldtype == "2d":  
+            dataset = ["fc","sfc","3"]
+            lev = [9999.]
         ######
         timchar = split_char (tim)

trunk/UTIL/PYTHON/myplot.py

@@ -43 +43 @@
     elif hasattr(nc,'START_DATE'):             typefile = 'meso' 
     elif 'HGT_M' in nc.variables:              typefile = 'geo'
+    elif hasattr(nc,'institution'):
+      if "European Centre" in getattr(nc,'institution'):  typefile = 'ecmwf'
     else:                                      typefile = 'gcm' # for lslin-ed files from gcm
     return typefile
@@ -449 +451 @@
     elif typefile in ['gcm']: 
         [lon2d,lat2d] = getcoord2d(nc,nlat="latitude",nlon="longitude",is1d=True)
-    elif typefile in ['earthgcm']: 
+    elif typefile in ['earthgcm','ecmwf']: 
         [lon2d,lat2d] = getcoord2d(nc,nlat="lat",nlon="lon",is1d=True)
     elif typefile in ['geo']:
@@ -886 +888 @@
                 "jupiter":     "http://www.mmedia.is/~bjj/data/jupiter_css/jupiter_css.jpg",\
                 "jupiter_voy": "http://www.mmedia.is/~bjj/data/jupiter/jupiter_vgr2.jpg",\
+                #"bw":          domain+"EarthElevation_2500x1250.jpg",\
                 "bw":          "http://users.info.unicaen.fr/~karczma/TEACH/InfoGeo/Images/Planets/EarthElevation_2500x1250.jpg",\
                 "contrast":    "http://users.info.unicaen.fr/~karczma/TEACH/InfoGeo/Images/Planets/EarthMapAtmos_2500x1250.jpg",\

trunk/UTIL/PYTHON/planetoplot.py

@@ -160 +160 @@
          else:                             [wlon,wlat] = simplinterv(lon2d,lat2d)
          if zoom:                          [wlon,wlat] = zoomset(wlon,wlat,zoom)
-         elif zarea is not None: [wlon,wlat] = latinterv(area=zarea)
+         elif zarea is not None:           [wlon,wlat] = latinterv(area=zarea)
 
 ##########################################################
 ############ LOAD 4D DIMENSIONS : x, y, z, t #############
 ##########################################################
-      if typefile in ["gcm","earthgcm"]:
-          ### SPACE
-          if typefile == "gcm":         lat = nc.variables["latitude"][:] ; lon = nc.variables["longitude"][:] ; vert = nc.variables["altitude"][:]
-          elif typefile == "earthgcm":  lat = nc.variables["lat"][:] ; lon = nc.variables["lon"][:] ; vert = nc.variables["Alt"][:]
-          if "aire" in nc.variables:      area = nc.variables["aire"][:,:]  #JL to weight means with the area
+      if typefile in ["gcm","earthgcm","ecmwf"]:
+          ### LAT
+          if "latitude" in nc.variables:  lat = nc.variables["latitude"][:]
+          elif "lat" in nc.variables:     lat = nc.variables["lat"][:]
+          ### LON
+          if "longitude" in nc.variables: lon = nc.variables["longitude"][:]
+          elif "lon" in nc.variables:     lon = nc.variables["lon"][:]
+          ### ALT
+          if "altitude" in nc.variables:  vert = nc.variables["altitude"][:]
+          elif "Alt" in nc.variables:     vert = nc.variables["Alt"][:]
+          elif "lev" in nc.variables:     vert = nc.variables["lev"][:]
+          else:                           vert = [0.]
+          ### AIRE (to weight means with the area)
+          if "aire" in nc.variables:      area = nc.variables["aire"][:,:]  
           else:                           area = None
           ### TIME