Changeset 608 for trunk/UTIL

Timestamp:

Mar 31, 2012, 2:12:55 PM (13 years ago)

Author:

aslmd

Message:

LMDZ.MARS: minor correction in solarlong. MESOSCALE: output mtot and icetot now same units, notes updated, improved save_this_simu so that all .def are saved. UTIL PYTHON: now in mesoscale mode, sections show lat/lon instead of nx or ny (to be improved), also corrected 1D labels that was not working for several time asked.

File:

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

trunk/UTIL/PYTHON/planetoplot.py

@@ -430 +430 @@
                    if len(what_I_plot.shape) == 3: itime=[0]
                    m = None ; x = None ; y = None
-                   what_I_plot, x, y = define_axis(lon,lat,vert,time,indexlon,indexlat,indexvert,\
+                   latyeah = lat ; lonyeah = lon
+                   if typefile in ['meso']:
+                       # now that the section is determined we can set the real lat
+                       # ... or for now, a temptative one.
+                       if slon is not None: latyeah = lat2d[:,0]
+                       if slat is not None: lonyeah = lon2d[0,:]
+                   what_I_plot, x, y = define_axis(lonyeah,latyeah,vert,time,indexlon,indexlat,indexvert,\
                          itime,what_I_plot, len(all_var[index_f].shape),vertmode)
                ###
@@ -491 +497 @@
                         if indextime is not None: lbl = lbl + " it" + str(indextime[0])
                         if lbl == "": lbl = namefiles[index_f]
-
-                        if lbls is not None: lbl=lbls[index_f]
+                        if lbls is not None: lbl=lbls[nplot-1]
 
                         if mrate is not None: x = y  ## because swapaxes...
@@ -548 +553 @@
                         zelevels = np.linspace(zevminc,zevmaxc,ticks/2) #20)
                         if var2 == 'HGT': zelevels = np.arange(-10000.,30000.,2000.)
-                        if mapmode == 0:    
+                        if mapmode == 0:   
                             what_I_plot_frame, x, y = define_axis( lon,lat,vert,time,indexlon,indexlat,indexvert,\
                                                               itime,what_I_plot_frame, len(all_var2[index_f].shape),vertmode )