Changeset 273 for trunk/MESOSCALE_DEV/PLOT

Timestamp:

Aug 30, 2011, 12:17:50 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: python improvements for graphics. nothing essential.

Location:

trunk/MESOSCALE_DEV/PLOT/PYTHON/mylib

Files:

• map_ecmwf.py (moved) (moved from trunk/MESOSCALE_DEV/PLOT/PYTHON/map_ecmwf.py) (2 diffs)
• myecmwf.py (modified) (1 diff)

trunk/MESOSCALE_DEV/PLOT/PYTHON/mylib/map_ecmwf.py

@@ -3 +3 @@
 ##########################################################################
 var = ["151","146","167"]
-#var = ["167"]
+var = ["167"]
+var = ["78"]
+var = ["137"]
 lev = ["700.","850."]
 tim = ["00"]
 date = ['10','08','2010','10','08','2010']
+date = ['01','09','2009','01','09','2009']
 dataset = ["an","pl","0"]
 dataset = ["fc","sfc","3"]
 ##########################################################################
 proj = "cyl"  #"moll" "ortho" "lcc"
-#proj = "ortho"
+proj = "ortho"
+#proj = "moll"
 area = "Europe"
+area = "Africa"
 #area = "Central_America"
 #area = "Southern_Hemisphere"
 #area = "Northern_Hemisphere"
+area = "Whole_No_High"
+area = "Whole"
+back="blue"
+#back="bw"
 ##########################################################################
 
@@ -29 +38 @@
 nc = mye.get_ecmwf (var, dataset, wlat, wlon, lev, date, tim) 
 ##########################################################################
-lat    = nc.variables['lat'             ][:]
-lon    = nc.variables['lon'             ][:]
-[lon2d,lat2d] = np.meshgrid(lon,lat)
+[lon2d,lat2d] = myp.getcoord2d (nc,nlat='lat',nlon='lon',is1d=True)
 step=10.
 ##########################################################################
 ntime = 0
 for i in range( np.array(var).size ):
-        for j in range( np.array(lev).size ):
-                for k in range( np.array(tim).size ):
-                        if dataset[1] == "pl":  field  = nc.variables['var'+var[i]      ][k,j,:,:]
-                        if dataset[1] == "sfc": field  = nc.variables['var'+var[i]      ][k,  :,:]
-                        m = myp.define_proj(proj,wlon,wlat)
-                        x, y = m(lon2d, lat2d)
-                        #m.drawmeridians(np.r_[wlon[0]:wlon[1]:step*2], labels=[0,0,0,1], color='grey')
-                        #m.drawparallels(np.r_[wlat[0]:wlat[1]:step  ], labels=[1,0,0,0], color='grey')
-                        m.bluemarble()
-                        #m.warpimage(myp.earthmap("nice"))
-                        zeplot = m.contour(x, y, field, 20)
-                        plt.clabel(zeplot, inline=1, inline_spacing=1, fontsize=7, fmt='%0i')
-                        #plt.colorbar(fraction=0.05,pad=0.1)
-                        plt.title(mye.ecmwf_title_field(var[i]))
-                        name = str(var[i])+str(lev[j])+str(tim[k])+".png"
-                        plt.savefig(name)
-                        plt.clf()
-                        yeah = myp.display(name)
+        for z in range( np.array(lev).size ):
+                for t in range( np.array(tim).size ):
+
+                    field, error = myp.reducefield( myp.getfield(nc,'var'+var[i]), d4=t, d3=z )
+                    if not error:
+                        ### Map projection
+                        m = myp.define_proj(proj,wlon,wlat,back=back)
+                        x, y = m(lon2d, lat2d)
+                        zeplot = m.contour(x, y, field, 15)
+                        plt.clabel(zeplot, inline=1, inline_spacing=1, fontsize=7, fmt='%0i')
+                        #plt.colorbar(fraction=0.05,pad=0.1)
+                        plt.title(mye.ecmwf_title_field(var[i]))
+                        myp.makeplotres(str(var[i])+str(lev[z])+str(tim[t]),res=100.,pad_inches_value=0.35)

trunk/MESOSCALE_DEV/PLOT/PYTHON/mylib/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 == "137": name="Total column water vapour (kg/m2)"
         return name