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

Timestamp:

Dec 1, 2011, 11:10:25 AM (13 years ago)

Author:

acolaitis

Message:

PYTHON. Contour plots now working in movies. I had to modify the call_contour (get out the reducefield part, introduce a new variable what_I_plot_contour dedicated to var2) in order to conserve the speed of movie rendering. (otherwise there were calls to reducefield for each frame in call_contour)

File:

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

trunk/UTIL/PYTHON/planetoplot.py

@@ -279 +279 @@
 
        #### Contour plot -- to be activated with movies
-       if var2 and mrate is None:
-           call_contour(lon,lat,vert,time,vertmode,all_var2[index_f],indextime,indexlon,indexlat,indexvert,yintegral,mapmode,typefile, var2, ticks)
+       if var2:
+           what_I_plot_contour, error = reducefield(all_var2[index_f], d4=indextime, d1=indexlon, d2=indexlat , d3=indexvert, yint=yintegral, alt=vert)
+           call_contour(what_I_plot_contour,error,lon,lat,vert,time,vertmode,all_var2[index_f],indextime,indexlon,indexlat,indexvert,yintegral,mapmode,typefile, var2, ticks)
 
        #### Shaded plot
@@ -329 +330 @@
                     if flagnolow:    what_I_plot_frame = nolow(what_I_plot_frame)
 
-                    # Renew axis directives for movie frames which are not the first one.
                     if imov > 0:
+                       # Renew axis directives for movie frames which are not the first one.
                        zxmin, zxmax = xaxis ; zymin, zymax = yaxis
                        if zxmin is not None: mpl.pyplot.xlim(xmin=zxmin)
@@ -338 +339 @@
                        if invert_y:     lima,limb = mpl.pyplot.ylim() ; mpl.pyplot.ylim(limb,lima)
                        if ylog:         mpl.pyplot.semilogy()
-                   
+                       # Renew contours for movie frames which are not the first one
+                       if var2:
+                         call_contour(what_I_plot_contour[imov,:,:],error,lon,lat,vert,time,vertmode,all_var2[index_f],imov,indexlon,indexlat,indexvert,yintegral,mapmode,typefile, var2, ticks)
                     if not tile:
                         #zelevels = np.linspace(zevmin*(1. + 1.e-7),zevmax*(1. - 1.e-7)) #,num=20)