Changeset 195 for trunk/MESOSCALE/PLOT/PYTHON

Timestamp:

Jul 5, 2011, 12:08:27 AM (14 years ago)

Author:

aslmd

Message:

MESOSCALE: graphics interface. corrected a bug in displaying automatic Ls

Location:

trunk/MESOSCALE/PLOT/PYTHON

Files:

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

trunk/MESOSCALE/PLOT/PYTHON/mylib/api_g95.sh

@@ -9 +9 @@
 \rm loge
 f2py -c -m api ../../../LMD_MM_MARS/SRC/POSTPROC/api.F90 --fcompiler=g95 -L$NETCDF/lib -lnetcdf -lm -I$NETCDF/include --f90flags="-Wall -Wno=112,141,137,155 -fno-second-underscore -ffree-form" > logc 2> loge
-f2py -c -m timestuff time.F --fcompiler=g95
+f2py -c -m timestuff time.F --fcompiler=g95 >> logc 2>> loge

trunk/MESOSCALE/PLOT/PYTHON/mylib/api_pgf90.sh

@@ -13 +13 @@
 \rm loge
 f2py -c -m api ../../../LMD_MM_MARS/SRC/POSTPROC/api.F90 --fcompiler=pg -L$NETCDF/lib -lnetcdf -lm -I$NETCDF/include --f90flags="-Mfree" > logc 2> loge
-f2py -c -m timestuff time.F --fcompiler=pg
+f2py -c -m timestuff time.F --fcompiler=pg >> logc 2>> loge

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

@@ -39 +39 @@
 #    return
 
+def getlschar ( namefile ):
+    #### strangely enough this does not work for api or ncrcat results!
+    from netCDF4 import Dataset
+    from timestuff import sol2ls
+    nc  = Dataset(namefile)
+    if 'Times' in nc.variables:
+        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 ) ) )
+    else:
+        lschar=""
+    return lschar
+
 def api_onelevel (  path_to_input   = None, \
                     input_name      = 'wrfout_d0?_????-??-??_??:00:00', \

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

@@ -17 +17 @@
            var=None,\
            colorb=None,\
-           winds=True):
+           winds=True,\
+           addchar=None):
 
     ####################################################################################################################
@@ -29 +30 @@
     from matplotlib.pyplot import contourf, subplot, figure, rcParams, savefig, colorbar
     from matplotlib.cm import get_cmap
-    from timestuff import sol2ls
     import numpy as np
 
@@ -75 +75 @@
     else:                    [wlon,wlat] = simplinterv(lon2d,lat2d)
 
-    #########################
-    ### Get local time and Ls
+    ##################
+    ### 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 ) ) )
+        #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
-        lschar = ""
 
     ##############################################################################
@@ -227 +226 @@
        
        ### Next subplot
-       ptitle( basename + "_LT"+str((ltst+i)%24) + lschar )
+       plottitle = basename + "_LT"+str((ltst+i)%24)
+       if addchar:  plottitle = plottitle + addchar
+       ptitle( plottitle )
        sub += 1
 
@@ -261 +262 @@
     from optparse import OptionParser    ### to be replaced by argparse
     from api_wrapper import api_onelevel
+    from netCDF4 import Dataset
+    from myplot import getlschar
 
     #############################
@@ -284 +287 @@
     zelevel = opt.nvert   
     stralt = None
+    lschar = getlschar ( zefile )  ## getlschar from wrfout (or simply return "" if another file)
 
     #####################################################
@@ -308 +312 @@
     ### Main call
     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)
+           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)
 
     #########################################################