Index: /trunk/UTIL/PYTHON/myplot.py
===================================================================
--- /trunk/UTIL/PYTHON/myplot.py	(revision 456)
+++ /trunk/UTIL/PYTHON/myplot.py	(revision 457)
@@ -864,4 +864,7 @@
         "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"
Index: /trunk/UTIL/PYTHON/myscript.py
===================================================================
--- /trunk/UTIL/PYTHON/myscript.py	(revision 456)
+++ /trunk/UTIL/PYTHON/myscript.py	(revision 457)
@@ -31,4 +31,5 @@
     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]
Index: /trunk/UTIL/PYTHON/planetoplot.py
===================================================================
--- /trunk/UTIL/PYTHON/planetoplot.py	(revision 456)
+++ /trunk/UTIL/PYTHON/planetoplot.py	(revision 457)
@@ -56,5 +56,6 @@
            mrate=None,\
            mquality=False,\
-           trans=1):
+           trans=1,\
+           zarea=None):
 
 
@@ -189,5 +190,6 @@
          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,9 +218,10 @@
         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,4 +247,5 @@
     fig = figure()
     subv,subh = definesubplot( numplot, fig ) 
+    if ope in ['-','-%']: subv,subh = 2,2
  
     #################################
Index: /trunk/UTIL/PYTHON/pp.py
===================================================================
--- /trunk/UTIL/PYTHON/pp.py	(revision 456)
+++ /trunk/UTIL/PYTHON/pp.py	(revision 457)
@@ -141,5 +141,5 @@
                 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")