Changeset 428 for trunk/UTIL/PYTHON/planetoplot.py

Timestamp:

Nov 30, 2011, 11:22:27 AM (13 years ago)

Author:

acolaitis

Message:

PYTHON. Very simple modifications to allow for slices with LES files. Performance is not optimized and turbulence computations are not implemented. More to come in time.

File:

• trunk/UTIL/PYTHON/planetoplot.py (modified) (5 diffs)

trunk/UTIL/PYTHON/planetoplot.py

@@ -92 +92 @@
     numplot = zelen*nslices
     print "********** FILES, SLICES, VARS, TOTAL PLOTS: ", len(namefiles), nslices, len(var), numplot
-    print "********** MAPMODE: ", mapmode
     if ope is not None:
         if fileref is not None:       zelen = zelen + 2
@@ -116 +115 @@
       ### ... TYPEFILE
       typefile = whatkindfile(nc)                                  
+      if typefile in ['les']:
+          mapmode=0
+          winds=False
+      if firstfile: print "********** MAPMODE: ", mapmode
       if firstfile:                 typefile0 = typefile
       elif typefile != typefile0:   errormess("Not the same kind of files !", [typefile0, typefile])
@@ -130 +133 @@
       [lon2d,lat2d] = getcoorddef(nc)
       ### ... PROJECTION
-      if proj == None:   proj = getproj(nc)                 
+      if ((proj == None) and (typefile not in ['les'])):   proj = getproj(nc)                 
 
 ##########################################################
@@ -139 +142 @@
           else:                           errormess("no time axis found.")
           vert = nc.variables["altitude"][:]
-      elif typefile in ['meso','mesoapi','geo']:
+      elif typefile in ['meso','mesoapi','geo','les']:
           if vlon is not None or vlat is not None:   indices = bidimfind(lon2d,lat2d,vlon,vlat) ; print '********** INDICES: ', indices
           if slon is not None: slon[0][0] = indices[0] ; slon[0][1] = indices[0]
@@ -329 +332 @@
                      #zelevels = np.linspace(zevmin*(1. + 1.e-7),zevmax*(1. - 1.e-7)) #,num=20)
                      zelevels = np.linspace(zevmin,zevmax,num=ticks)
+                     print np.array(x).shape
+                     print np.array(y).shape
+                     print np.array(what_I_plot).shape
+
                      if mapmode == 1:       m.contourf( x, y, what_I_plot, zelevels, cmap = palette)
                      elif mapmode == 0:     contourf( x, y, what_I_plot, zelevels, cmap = palette)