Index: trunk/UTIL/PYTHON/planetoplot.py
===================================================================
--- trunk/UTIL/PYTHON/planetoplot.py	(revision 437)
+++ trunk/UTIL/PYTHON/planetoplot.py	(revision 438)
@@ -288,5 +288,5 @@
               if var2 == 'HGT':  zelevels = np.arange(-10000.,30000.,2000.)
               if mapmode == 0:
-                  if typefile in ['mesoideal']:    what_I_plot = dumpbdy(what_I_plot,0,stag='W')
+                  #if typefile in ['mesoideal']:    what_I_plot = dumpbdy(what_I_plot,0,stag='W')
                   itime=indextime 
                   if len(what_I_plot.shape) is 3: itime=[0]
@@ -295,5 +295,5 @@
               ### If we plot a 2-D field
               if len(what_I_plot.shape) is 2:
-                  #zelevels=[-10.,-8.,-6.,-4.,-2.,0.,2.,4.,6.,8.,10.]
+                  #zelevels=[0.2,0.4,0.6,0.8,1.]
                   cs = contour(x,y,what_I_plot, zelevels, colors='k', linewidths = 1 ) #0.33 colors='w' )# , alpha=0.5)
                   #clabel(cs,zelevels,inline=3,fmt='%1.1f',fontsize=7)
@@ -310,4 +310,5 @@
                                              yint=yintegral, alt=vert, anomaly=anomaly )
            what_I_plot = what_I_plot*mult
+           #if typefile in ['mesoideal']:    what_I_plot = dumpbdy(what_I_plot,0,stag='W')
            if not error: 
                fvar = varname 
@@ -321,8 +322,9 @@
                    if zxmin is not None: mpl.pyplot.xlim(xmin=zxmin)
                    if zxmax is not None: mpl.pyplot.xlim(xmax=zxmax)
-                   if zymin is not None: mpl.pyplot.ylim(ymin=zymin) 
+                   if zymin is not None: mpl.pyplot.ylim(ymin=zymin)
                    if zymax is not None: mpl.pyplot.ylim(ymax=zymax)
                    if invert_y:     lima,limb = mpl.pyplot.ylim() ; mpl.pyplot.ylim(limb,lima)
                    if ylog:         mpl.pyplot.semilogy()
+
                if (fileref is not None) and (index_f is numplot-1):    zevmin, zevmax = calculate_bounds(what_I_plot,vmin=minop,vmax=maxop)
                else:                                                   zevmin, zevmax = calculate_bounds(what_I_plot,vmin=vmin,vmax=vmax)
@@ -348,4 +350,16 @@
                     else:            what_I_plot_frame = bounds(what_I_plot_frame,zevmin,zevmax)
                     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:
+                       zxmin, zxmax = xaxis ; zymin, zymax = yaxis
+                       if zxmin is not None: mpl.pyplot.xlim(xmin=zxmin)
+                       if zxmax is not None: mpl.pyplot.xlim(xmax=zxmax)
+                       if zymin is not None: mpl.pyplot.ylim(ymin=zymin)
+                       if zymax is not None: mpl.pyplot.ylim(ymax=zymax)
+                       if invert_y:     lima,limb = mpl.pyplot.ylim() ; mpl.pyplot.ylim(limb,lima)
+                       if ylog:         mpl.pyplot.semilogy()
+
+
                     if not tile:
                         #zelevels = np.linspace(zevmin*(1. + 1.e-7),zevmax*(1. - 1.e-7)) #,num=20)
@@ -362,7 +376,4 @@
                         colorbar( fraction=0.05,pad=0.03,format=daformat,\
                                   ticks=np.linspace(zevmin,zevmax,num=min([ticks/2+1,20])),extend='neither',spacing='proportional' ) 
-                        #mframe=mpl.pyplot.gci().to_rgba(mpl.pyplot.gci().get_array())
-                        #mframe=mpl.pyplot.imshow()
-                        #mframe=fig2data(mpl.pyplot.gcf())
                         mframe=fig2img(mpl.pyplot.gcf())
                     if ((mrate is not None) and (imov is 0)):