Changeset 385 for trunk/UTIL

Timestamp:

Nov 14, 2011, 9:56:08 PM (14 years ago)

Author:

acolaitis

Message:

PYTHON. Minor commit to avoid loosing stuff when gnome is down tomorrow.

Location:

trunk/UTIL/PYTHON

Files:

• gcm.py (modified) (1 diff)
• meso.py (modified) (1 diff)
• myplot.py (modified) (2 diffs)
• myscript.py (modified) (1 diff)
• planetoplot.py (modified) (2 diffs)
• zrecast_wrapper.py (modified) (1 diff)

trunk/UTIL/PYTHON/gcm.py

@@ -175 +175 @@
                 outputname=opt.out,resolution=opt.res,\
                 ope=opt.operat,fileref=reffile,minop=opt.vminope,maxop=opt.vmaxope,titleref=opt.titref,\
-                invert_y=opt.inverty,xaxis=zexaxis,yaxis=zeyaxis,ylog=opt.logy,yintegral=opt.column)
+                invert_y=opt.inverty,xaxis=zexaxis,yaxis=zeyaxis,ylog=opt.logy,yintegral=opt.column,\
+                blat=opt.blat)
         print 'Done: '+name
         system("rm -f to_be_erased")

trunk/UTIL/PYTHON/meso.py

@@ -172 +172 @@
                 outputname=opt.out,resolution=opt.res,\
                 ope=opt.operat,fileref=reffile,minop=opt.vminope,maxop=opt.vmaxope,titleref=opt.titref,\
-                invert_y=opt.inverty,xaxis=zexaxis,yaxis=zeyaxis,ylog=opt.logy,yintegral=opt.column)
+                invert_y=opt.inverty,xaxis=zexaxis,yaxis=zeyaxis,ylog=opt.logy,yintegral=opt.column,\
+                blat=opt.blat)
         print 'Done: '+name
         system("rm -f to_be_erased")

trunk/UTIL/PYTHON/myplot.py

@@ -444 +444 @@
 
 ## Author: AS
-def define_proj (char,wlon,wlat,back=None,blat=False):
+def define_proj (char,wlon,wlat,back=None,blat=None):
     from    mpl_toolkits.basemap            import Basemap
     import  numpy                           as np
@@ -451 +451 @@
     meanlon = 0.5*(wlon[0]+wlon[1])
     meanlat = 0.5*(wlat[0]+wlat[1])
-    if not blat:
+    if blat is None:
         if   wlat[0] >= 80.:   blat =  40. 
         elif wlat[1] <= -80.:  blat = -40.

trunk/UTIL/PYTHON/myscript.py

@@ -38 +38 @@
     parser.add_option('-e', '--itstep', action='store',dest='it',        type="int",     default=None,  help='stride time [4]')
     parser.add_option('-F', '--first',  action='store',dest='frt',       type="int",     default=1,     help='first subscript to plot [1]')
+    parser.add_option('--blat',         action='store',dest='blat',      type="int",     default=None,  help='bounding latitude for stereographic plots [computed]')
 
     ### SPECIFIC FOR SLICING [MAPMODE 0]

trunk/UTIL/PYTHON/planetoplot.py

@@ -49 +49 @@
            yaxis=[None,None],\
            ylog=False,\
-           yintegral=False):
+           yintegral=False,\
+           blat=None):
 
 
@@ -232 +233 @@
        ### Map projection                    
        if mapmode == 1:
-           m = define_proj(proj,wlon,wlat,back=back)
+           m = define_proj(proj,wlon,wlat,back=back,blat=blat)
            x, y = m(lon2d, lat2d)
 

trunk/UTIL/PYTHON/zrecast_wrapper.py

@@ -28 +28 @@
              system("echo 1 >> zrecast.auto.def")
              system("echo yes >> zrecast.auto.def")
-             system("echo 16.58 0.11 >> zrecast.auto.def")  #I put that randomly! (a.c.)
-             system("echo 32 >> zrecast.auto.def")
+             system("echo 370 0.1 >> zrecast.auto.def")  #I put that randomly! (a.c.)
+             system("echo 20 >> zrecast.auto.def")
         else:
              print "zrecast interp option unsupported for now. Exiting."