Changeset 457 for trunk/UTIL

Timestamp:

Dec 7, 2011, 11:56:29 AM (13 years ago)

Author:

acolaitis

Message:

PYTHON.

• added possibility to specify area of plotting, using --area "name_of_area" example --area Sirenum_Crater_small list of areas are in myplot: latinterv (feel free to add more)

• added operation "-%": difference in percent between two files, with respect to the reference file. [ 100.*(file-ref)/ref ]

• changed default subplot tile for operations "-" and "-%" to 2x2 instead of 3x1

Location:

trunk/UTIL/PYTHON

Files:

• myplot.py (modified) (1 diff)
• myscript.py (modified) (1 diff)
• planetoplot.py (modified) (4 diffs)
• pp.py (modified) (1 diff)

trunk/UTIL/PYTHON/myplot.py

@@ -864 +864 @@
         "South_Pole":            [[-90.,-50.],[-180., 180.]],\
         "Close_South_Pole":      [[-90.,-75.],[-180., 180.]],\
+        "Sirenum_Crater_large":  [[-46.,-34.],[-166., -151.]],\
+        "Sirenum_Crater_small":  [[-36.,-26.],[-168., -156.]],\
+
               }
     if area not in list:   area = "Whole"

trunk/UTIL/PYTHON/myscript.py

@@ -31 +31 @@
     parser.add_option('--res',          action='store',dest='res',       type="float",   default=200.,  help='Resolution for png outputs. --save png needed. [200.]')
     parser.add_option('--trans',        action='store',dest='trans',     type="float",   default=1.,    help='shaded plots transparency, from 0 (transparent) to 1 (opaque) [1]')
+    parser.add_option('--area',         action='store',dest='area',       type="string",   default=None,  help='Area on the map to be plot [None]')
 
     ### SPECIFIC FOR MAPPING [MAPMODE 1]

trunk/UTIL/PYTHON/planetoplot.py

@@ -56 +56 @@
            mrate=None,\
            mquality=False,\
-           trans=1):
+           trans=1,\
+           zarea=None):
 
 
@@ -189 +190 @@
          elif proj in ["lcc","laea"]:      [wlon,wlat] = wrfinterv(lon2d,lat2d)
          else:                             [wlon,wlat] = simplinterv(lon2d,lat2d)
-         if zoom:                          [wlon,wlat] = zoomset(wlon,wlat,zoom) 
+         if zoom:                          [wlon,wlat] = zoomset(wlon,wlat,zoom)
+         elif zarea is not None: [wlon,wlat] = latinterv(area=zarea) 
 
       all_varname[k] = varname
@@ -216 +218 @@
         if "var" not in ope:
              if len(var) > 1: errormess("for this operation... please set only one var !")
-             if ope in ["-","+"]:
+             if ope in ["-","+","-%"]:
                 if fileref is not None:   all_var[k] = getfield(Dataset(fileref),all_varname[k-1]) ; all_varname[k] = all_varname[k-1] ; all_time[k] = all_time[k-1] ; all_namefile[k] = all_namefile[k-1]
                 else:                     errormess("fileref is missing!") 
                 if ope == "-":     all_var[k+1]= all_var[k-1] - all_var[k]
                 elif ope == "+":   all_var[k+1]= all_var[k-1] + all_var[k]
+                elif ope == "-%":  all_var[k+1]= 100.*(all_var[k-1] + all_var[k])/all_var[k]
                 all_varname[k+1] = all_varname[k] ; all_time[k+1] = all_time[k] ; all_namefile[k+1] = all_namefile[k] ; numplot = numplot+2
              elif ope in ["cat"]:
@@ -244 +247 @@
     fig = figure()
     subv,subh = definesubplot( numplot, fig ) 
+    if ope in ['-','-%']: subv,subh = 2,2
  
     #################################

trunk/UTIL/PYTHON/pp.py

@@ -141 +141 @@
                 invert_y=opt.inverty,xaxis=zexaxis,yaxis=zeyaxis,ylog=opt.logy,yintegral=opt.column,\
                 blat=opt.blat,blon=opt.blon,tsat=opt.tsat,flagnolow=opt.nolow,\
-                mrate=opt.rate,mquality=opt.quality,trans=opt.trans)
+                mrate=opt.rate,mquality=opt.quality,trans=opt.trans,zarea=opt.area)
         print 'DONE: '+name
         system("rm -f to_be_erased")