Changeset 252 for trunk/MESOSCALE_DEV/PLOT

Timestamp:

Aug 2, 2011, 1:48:14 AM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: oops. forgot a file in the previous commit. all apologies.

File:

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

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

@@ -1 @@
-def errormess(text):
+def errormess(text,printvar=None):
     print text
+    if printvar: print printvar 
     exit()
     return
+
+def getname(var=False,winds=False,anomaly=False):
+    if var and winds:     basename = var + '_UV'
+    elif var:             basename = var
+    elif winds:           basename = 'UV'
+    else:                 errormess("please set at least winds or var",printvar=nc.variables)
+    if anomaly:           basename = 'd' + basename
+    return basename
+
+def localtime(utc,lon):
+    ltst = utc + lon / 15.
+    ltst = int (ltst * 10) / 10.
+    ltst = ltst % 24
+    return ltst
 
 def whatkindfile (nc):
@@ -21 +36 @@
 
 def reducefield (input,d4=None,d3=None,d2=None,d1=None):
-    ### do it the reverse way to be compliant with netcdf "t z y x" or "t y x" or "y x"
+    ### we do it the reverse way to be compliant with netcdf "t z y x" or "t y x" or "y x"
     ### it would be actually better to name d4 d3 d2 d1 as t z y x
     import numpy as np
@@ -71 +86 @@
     return output, error
 
-#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 ):
     from matplotlib.pyplot import rcParams
@@ -197 +195 @@
         else:
             proj="merc"
-    elif typefile in ['gcm']:  proj="cyl"  ## pb avec les autres (de trace derriere la sphere ?)
+    elif typefile in ['gcm']:  proj="cyl"    ## pb avec les autres (de trace derriere la sphere ?)
     else:                      proj="ortho"
     return proj    
@@ -205 +203 @@
     title(name)
     print name
+
+def polarinterv (lon2d,lat2d):
+    import numpy as np
+    wlon = [np.min(lon2d),np.max(lon2d)]
+    ind = np.array(lat2d).shape[0] / 2  ## to get a good boundlat and to get the pole
+    wlat = [np.min(lat2d[ind,:]),np.max(lat2d[ind,:])]
+    return [wlon,wlat]
 
 def simplinterv (lon2d,lat2d):
@@ -365 +370 @@
     elif wlat[1] >= 0.:    blat = wlat[0] 
     elif wlat[0] <= 0.:    blat = wlat[1]
-    print blat
+    print "blat ", blat
     h = 50.  ## en km
     radius = 3397200.
@@ -477 +482 @@
              "ICETOT":       "%.1e",\
              "TAU_ICE":      "%.2f",\
+             "VMR_ICE":      "%.1e",\
+             "MTOT":         "%.0f",\
              "anomaly":      "%.1f",\
              "W":            "%.1f",\
@@ -495 +502 @@
              "HFX":          "RdYlBu",\
              "ICETOT":       "YlGnBu",\
+             "MTOT":         "PuBu",\
              "TAU_ICE":      "Blues",\
+             "VMR_ICE":      "Blues",\
              "W":            "jet",\
              "anomaly":      "RdBu_r",\