Changeset 180 for trunk/MESOSCALE/PLOT/PYTHON

Timestamp:

Jun 29, 2011, 5:44:33 PM (14 years ago)

Author:

aslmd

Message:

MESOSCALE: PYTHON: ajouts librairies maison: myplot mymath myecmwf. ajouts scripts maison: domain winds

Location:

trunk/MESOSCALE/PLOT/PYTHON

Files:

• map_ecmwf.py (modified) (1 diff)
• mylib (added)
• mylib/ecmwf.py (added)
• mylib/myecmwf.py (added)
• mylib/mymath.py (added)
• mylib/myplot.py (added)
• scripts (added)
• scripts/domain.py (moved) (moved from trunk/MESOSCALE/PLOT/PYTHON/domain.py) (2 diffs)
• scripts/winds.py (added)

trunk/MESOSCALE/PLOT/PYTHON/map_ecmwf.py

@@ -11 +11 @@
 ##########################################################################
 proj = "cyl"  #"moll" "ortho" "lcc"
-proj = "ortho"
+#proj = "ortho"
 area = "Europe"
-area = "Central_America"
-area = "Southern_Hemisphere"
-area = "Northern_Hemisphere"
+#area = "Central_America"
+#area = "Southern_Hemisphere"
+#area = "Northern_Hemisphere"
 ##########################################################################
 

trunk/MESOSCALE/PLOT/PYTHON/scripts/domain.py

@@ -8 +8 @@
     print 'Example:  domain.py /tmp7/aslmd/HELLAS/wrfout_d02_2024-10-03_06:00:00'
 
-def domain (namefile): 
+def domain (namefile,proj="ortho",back="molabw"): #vishires
     from netCDF4 import Dataset
     from myplot import getcoord2d,define_proj,makeplotpng,simplinterv
@@ -15 +15 @@
     [lon2d,lat2d] = getcoord2d(nc)
     [wlon,wlat] = simplinterv(lon2d,lat2d)
-    m = define_proj("ortho",wlon,wlat,back="vishires")
+    m = define_proj(proj,wlon,wlat,back=back)
     x, y = m(lon2d, lat2d)
-    contourf(x, y, nc.variables['HGT'][0,:,:] / 1000.,50)
-    makeplotpng(namefile+".domain",pad_inches_value=0.15)
+    contourf(x, y, nc.variables['HGT'][0,:,:] / 1000., 50)
+    makeplotpng(namefile+".domain",pad_inches_value=0.05)
 
 if __name__ == "__main__":