Changeset 197

Timestamp:

Jul 5, 2011, 8:56:29 AM (14 years ago)

Author:

aslmd

Message:

MESOSCALE: python graphics. fix the program so that it guesses the correct local time.

Location:

trunk/MESOSCALE/PLOT/PYTHON

Files:

• mylib/myplot.py (modified) (1 diff)
• scripts/winds.py (modified) (7 diffs)

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

@@ -50 +50 @@
         if zeday < 0:    lschar=""  ## might have crossed a month... fix soon
         else:            lschar="_Ls"+str( int( sol2ls ( getattr( nc, 'JULDAY' ) + zeday ) ) )
+        ###
+        zetime2 = nc.variables['Times'][1]
+        one  = int(zetime[11]+zetime[12]) + int(zetime[14]+zetime[15])/37.
+        next = int(zetime2[11]+zetime2[12]) + int(zetime2[14]+zetime2[15])/37. 
+        zehour    = one
+        zehourin  = abs ( next - one )
     else:
         lschar=""
-    return lschar
+        zehour = 0
+        zehourin = 1  
+    return lschar, zehour, zehourin
 
 def api_onelevel (  path_to_input   = None, \

trunk/MESOSCALE/PLOT/PYTHON/scripts/winds.py

@@ -18 +18 @@
            colorb=None,\
            winds=True,\
-           addchar=None):
+           addchar=None,
+           interv=[0,1]):
 
     ####################################################################################################################
@@ -75 +76 @@
     else:                    [wlon,wlat] = simplinterv(lon2d,lat2d)
 
-    ##################
-    ### Get local time
-    if typefile in ['mesoapi','meso']:  
-        ltst = int(getattr(nc, 'GMT') + 0.5*(wlon[0]+wlon[1])/15.)
-        #zetime = nc.variables['Times'][0]
-        #zetimestart = getattr(nc, 'START_DATE')
-        #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 ) ) )
-    elif typefile in ['gcm']:           
-        ltst = 0
+    ###################
+    #### Get local time
+    #if typefile in ['mesoapi','meso']:  
+    #    ltst = interv[0]
+    #    #ltst = int(getattr(nc, 'GMT') + 0.5*(wlon[0]+wlon[1])/15.)
+    #elif typefile in ['gcm']:           
+    #    ltst = 0
 
     ##############################################################################
@@ -185 +182 @@
     for i in tabrange:
 
+       ### Which local time ?
+       ltst = ( interv[0] + 0.5*(wlon[0]+wlon[1])/15.) + i*interv[1]
+       ltst = int (ltst * 10) / 10.
+       ltst = ltst % 24
+
        ### General plot settings
        if numplot > 1: 
@@ -193 +195 @@
             ### If only one local time is requested (numplot < 0)
        elif numplot <= 0: 
-           if (ltst+i)%24 + numplot != 0:   continue
+           if ltst + numplot != 0:   continue
            else:                            found_lct = True
 
@@ -226 +228 @@
        
        ### Next subplot
-       plottitle = basename + "_LT"+str((ltst+i)%24)
+       plottitle = basename + "_LT"+str(ltst)
        if addchar:  plottitle = plottitle + addchar
        ptitle( plottitle )
@@ -287 +289 @@
     zelevel = opt.nvert   
     stralt = None
-    lschar = getlschar ( zefile )  ## getlschar from wrfout (or simply return "" if another file)
+    [lschar,zehour,zehourin] = getlschar ( zefile )  ## getlschar from wrfout (or simply return "" if another file)
 
     #####################################################
@@ -313 +315 @@
     name = winds (zefile,int(zelevel),\
            proj=opt.proj,back=opt.back,target=opt.target,stride=opt.stride,var=opt.var,numplot=opt.numplot,colorb=opt.colorb,winds=opt.winds,\
-           addchar=lschar)
+           addchar=lschar,interv=[zehour,zehourin])
 
     #########################################################