Changeset 241 for trunk/MESOSCALE_DEV/PLOT/PYTHON/mylib

Timestamp:

Jul 23, 2011, 1:38:02 AM (14 years ago)

Author:

aslmd

Message:

MESOSCALE: python graphics: contour plots + various small but convenient improvements. fixed broken domain.py

File:

• trunk/MESOSCALE_DEV/PLOT/PYTHON/mylib/myplot.py (modified) (7 diffs)

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

@@ -71 +71 @@
     return output, error
 
-def latinterv (area):
-        if   area == "Europe": 
-                wlat = [20.,80.]
-                wlon = [-50.,50.]
-        elif area == "Central_America":
-                wlat = [-10.,40.]
-                wlon = [230.,300.]
-        elif area == "Africa":
-                wlat = [-20.,50.]
-                wlon = [-50.,50.]
-        elif area == "Whole":
-                wlat = [-90.,90.]
-                wlon = [-180.,180.]
-        elif area == "Southern_Hemisphere":
-                wlat = [-90.,60.]
-                wlon = [-180.,180.]
-        elif area == "Northern_Hemisphere":
-                wlat = [-60.,90.]
-                wlon = [-180.,180.]
-        elif area == "Tharsis":
-                wlat = [-30.,60.]
-                wlon = [-170.,-10.]
-        elif area == "Whole_No_High":
-                wlat = [-60.,60.]
-                wlon = [-180.,180.]
-        elif area == "Chryse":
-                wlat = [-60.,60.]
-                wlon = [-60.,60.]
-        elif area == "North_Pole":
-                wlat = [50.,90.]
-                wlon = [-180.,180.]
-        elif area == "Close_North_Pole":
-                wlat = [75.,90.]
-                wlon = [-180.,180.]
-        elif area == "South_Pole":
-                wlat = [-90.,-50.]
-                wlon = [-180.,180.]
-        elif area == "Close_South_Pole":
-                wlat = [-90.,-75.]
-                wlon = [-180.,180.]
-        return wlon,wlat
+#def latinterv (area):
+#       if   area == "Europe":                [wlat,wlon] = [[ 20., 80.],[- 50.,  50.]]
+#       elif area == "Central_America":       [wlat,wlon] = [[-10., 40.],[ 230., 300.]]
+#       elif area == "Africa":                [wlat,wlon] = [[-20., 50.],[- 50.,  50.]]
+#       elif area == "Whole":                 [wlat,wlon] = [[-90., 90.],[-180.,-180.]]
+#       elif area == "Southern_Hemisphere":   [wlat,wlon] = [[-90., 60.],[-180.,-180.]]
+#        elif area == "Northern_Hemisphere":   [wlat,wlon] = [[-60., 90.],[-180.,-180.]]
+#        elif area == "Tharsis":               [wlat,wlon] = [[-30., 60.],[-170.,- 10.]]
+#       elif area == "Whole_No_High":         [wlat,wlon] = [[-60., 60.],[-180., 180.]]
+#       elif area == "Chryse":                [wlat,wlon] = [[-60., 60.],[- 60.,  60.]]
+#        elif area == "North_Pole":            [wlat,wlon] = [[ 50., 90.],[-180., 180.]]
+#        elif area == "Close_North_Pole":      [wlat,wlon] = [[ 75., 90.],[-180., 180.]]
+#        elif area == "Far_South_Pole":        [wlat,wlon] = [[-90.,-40.],[-180., 180.]]
+#        elif area == "South_Pole":            [wlat,wlon] = [[-90.,-50.],[-180., 180.]]
+#        elif area == "Close_South_Pole":      [wlat,wlon] = [[-90.,-75.],[-180., 180.]]
+#       return wlon,wlat
 
 def definesubplot ( numplot, fig ):
@@ -142 +117 @@
     elif numplot == 1:
         sub = 99999
-    elif numplot < 0: 
+    elif numplot <= 0: 
         sub = 99999
     else:
@@ -182 +157 @@
         zeday = int(zetime[8]+zetime[9]) - int(zetimestart[8]+zetimestart[9])
         if zeday < 0:    lschar=""  ## might have crossed a month... fix soon
-        else:            lschar="_Ls"+str( int( sol2ls ( getattr( nc, 'JULDAY' ) + zeday ) ) )
+        else:            lschar="_Ls"+str( int( 10. * sol2ls ( getattr( nc, 'JULDAY' ) + zeday ) ) / 10. )
         ###
         zetime2 = nc.variables['Times'][1]
@@ -472 +447 @@
     return what_I_plot
 
+def nolow(what_I_plot):
+    from mymath import max,min
+    lim = 0.15*0.5*(abs(max(what_I_plot))+abs(min(what_I_plot)))
+    print "on vire en dessous de ", lim
+    what_I_plot [ abs(what_I_plot) < lim ] = 1.e40 
+    return what_I_plot
+
 def zoomset (wlon,wlat,zoom):
     dlon = abs(wlon[1]-wlon[0])/2.
@@ -492 +474 @@
              "TAU_ICE":      "%.2f",\
              "anomaly":      "%.1f",\
+             "W":            "%.1f",\
                     }
     if whichvar not in fmtvar:
@@ -509 +492 @@
              "ICETOT":       "YlGnBu",\
              "TAU_ICE":      "Blues",\
+             "W":            "jet",\
              "anomaly":      "RdBu_r",\
                      }
+#W --> spectral ou jet
 #spectral BrBG RdBu_r
+    print "predefined colorbars"
     if whichone not in whichcolorb:
         whichone = "def"
@@ -567 +553 @@
         return whichlink
 
+def latinterv (area="Whole"):
+    list =    { \
+        "Europe":                [[ 20., 80.],[- 50.,  50.]],\
+        "Central_America":       [[-10., 40.],[ 230., 300.]],\
+        "Africa":                [[-20., 50.],[- 50.,  50.]],\
+        "Whole":                 [[-90., 90.],[-180.,-180.]],\
+        "Southern_Hemisphere":   [[-90., 60.],[-180.,-180.]],\
+        "Northern_Hemisphere":   [[-60., 90.],[-180.,-180.]],\
+        "Tharsis":               [[-30., 60.],[-170.,- 10.]],\
+        "Whole_No_High":         [[-60., 60.],[-180., 180.]],\
+        "Chryse":                [[-60., 60.],[- 60.,  60.]],\
+        "North_Pole":            [[ 50., 90.],[-180., 180.]],\
+        "Close_North_Pole":      [[ 75., 90.],[-180., 180.]],\
+        "Far_South_Pole":        [[-90.,-40.],[-180., 180.]],\
+        "South_Pole":            [[-90.,-50.],[-180., 180.]],\
+        "Close_South_Pole":      [[-90.,-75.],[-180., 180.]],\
+              }
+    if area not in list:   area = "Whole"
+    [olat,olon] = list[area]
+    return olon,olat
+