Changeset 194 for trunk/MESOSCALE/PLOT/PYTHON

Timestamp:

Jul 4, 2011, 6:13:23 PM (14 years ago)

Author:

aslmd

Message:

MESOSCALE: corrected graphics script to display Ls. and added a figure in the first page of the user manual.

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) (3 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

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

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

@@ -180 +180 @@
     else:                  kcolor=color
     if key: p = plt.quiverkey(q,posx,posy,scale,\
-                   str(int(scale)),coordinates='data',color=kcolor,labelpos='S')
+                   str(int(scale)),coordinates='data',color=kcolor,labelpos='S',labelsep = 0.03)
     return 
 

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

@@ -26 +26 @@
     from netCDF4 import Dataset
     from myplot import getcoord2d,define_proj,makeplotpng,simplinterv,vectorfield,ptitle,latinterv,getproj,wrfinterv,dumpbdy
+    from mymath import deg
     from matplotlib.pyplot import contourf, subplot, figure, rcParams, savefig, colorbar
     from matplotlib.cm import get_cmap
+    from timestuff import sol2ls
     import numpy as np
 
@@ -73 +75 @@
     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.)
-    elif typefile in ['gcm']:           ltst = 0
+    #########################
+    ### Get local time and Ls
+    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
+        lschar = ""
 
     ##############################################################################
@@ -217 +227 @@
        
        ### Next subplot
-       ptitle( basename + "_LT"+str((ltst+i)%24) )
+       ptitle( basename + "_LT"+str((ltst+i)%24) + lschar )
        sub += 1