Changeset 225 for trunk/MESOSCALE_DEV/PLOT/PYTHON/mylib

Timestamp:

Jul 15, 2011, 2:47:21 PM (14 years ago)

Author:

aslmd

Message:

MESOSCALE: python graphics. now possible for multi-files and multi-var. very easy to make an atlas.

Location:

trunk/MESOSCALE_DEV/PLOT/PYTHON/mylib

Files:

• api_g95.sh (modified) (1 diff)
• api_pgf90.sh (modified) (2 diffs)
• api_wrapper.py (modified) (1 diff)
• myplot.py (modified) (5 diffs)

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

@@ -8 +8 @@
 \rm logc
 \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 api 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 >> logc 2>> loge

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

@@ -5 +5 @@
 #NETCDF=/distrib/local/netcdf-4.0.1/pgi_7.1-6_64/
 #NETCDF=/distrib/local/netcdf-3.6.0-p1/pgi_64bits/
-NETCDF=/donnees/aslmd/MODELES/MESOSCALE/NETCDF/pgf90_64_fpic/netcdf-3.6.1/
+NETCDF=/donnees/aslmd/MODELES/MESOSCALE_DEV/NETCDF/pgf90_64_netcdf_fpic-3.6.1/
+
 echo $NETCDF
 
@@ -12 +13 @@
 \rm logc
 \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 api 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 >> logc 2>> loge

trunk/MESOSCALE_DEV/PLOT/PYTHON/mylib/api_wrapper.py

@@ -26 +26 @@
         if interp_method == 4:    output_name = input_name+'_zabg'
 
-    print input_name, output_name
+    #print input_name, output_name
 
     if nocall:     pass

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

@@ -44 +44 @@
     from timestuff import sol2ls
     nc  = Dataset(namefile)
-    if 'Times' in nc.variables and 'vert' not in nc.variables:
+    if namefile[0]+namefile[1]+namefile[2] != "geo" \
+       and 'Times' in nc.variables \
+       and 'vert' not in nc.variables:
         zetime = nc.variables['Times'][0]
         zetimestart = getattr(nc, 'START_DATE')
@@ -100 +102 @@
     nx = len(lon2d[0,:])-1
     ny = len(lon2d[:,0])-1
-    return [[lon2d[0,0],lon2d[nx,ny]],[lat2d[0,0],lat2d[nx,ny]]]
+    lon1 = lon2d[0,0] 
+    lon2 = lon2d[nx,ny] 
+    lat1 = lat2d[0,0] 
+    lat2 = lat2d[nx,ny] 
+    wider = 0.5 * (abs(lon1)+abs(lon2)) * 0.1
+    if lon1 < lon2:  wlon = [lon1, lon2 + wider]  ## a tester en normal
+    else:            wlon = [lon2, lon1 + wider]
+    if lat1 < lat2:  wlat = [lat1, lat2]
+    else:            wlat = [lat2, lat1]
+    return [wlon,wlat]
 
 def makeplotpngres (filename,res,pad_inches_value=0.25,folder='',disp=True):
@@ -237 +248 @@
     if char in ["cyl","lcc","merc","nsper","laea"]:   step = findstep(wlon)
     else:                                             step = 10.
+    print step
     m.drawmeridians(np.r_[-180.:180.:step*2.], labels=[0,0,0,1], color='grey', fontsize=fontsizemer)
     m.drawparallels(np.r_[-90.:90.:step], labels=[1,0,0,0], color='grey', fontsize=fontsizemer)
@@ -252 +264 @@
              "HGT":          "%.1e",\
              "USTM":         "%.2f",\
+             "HFX":          "%.0f",\
                     }
     if whichvar not in fmtvar:
@@ -264 +277 @@
              "QH2O":         "PuBu",\
              "USTM":         "YlOrRd",\
-#"RdPu",\
+             "HFX":          "RdYlBu",\
                      }
     if whichone not in whichcolorb: