Index: trunk/UTIL/PYTHON/myplot.py
===================================================================
--- trunk/UTIL/PYTHON/myplot.py	(revision 431)
+++ trunk/UTIL/PYTHON/myplot.py	(revision 432)
@@ -39,5 +39,6 @@
     if 'controle' in nc.variables:             typefile = 'gcm'
     elif 'phisinit' in nc.variables:           typefile = 'gcm'
-    elif '9999' in getattr(nc,'START_DATE') :  typefile = 'mesoideal'
+    elif hasattr(nc,'START_DATE'):
+      if '9999' in getattr(nc,'START_DATE') :  typefile = 'mesoideal'
     elif 'vert' in nc.variables:               typefile = 'mesoapi'
     elif 'U' in nc.variables:                  typefile = 'meso'
Index: trunk/UTIL/PYTHON/myscript.py
===================================================================
--- trunk/UTIL/PYTHON/myscript.py	(revision 431)
+++ trunk/UTIL/PYTHON/myscript.py	(revision 432)
@@ -61,4 +61,5 @@
     ### SPECIAL
     parser.add_option('--tsat',         action='store_true',dest='tsat',               default=False,help='convert temperature field T in Tsat-T using pressure')
+    parser.add_option('--rate',         action='store'     ,dest='rate', type="int",   default=None, help='Output is a movie along Time dimension. --time must not be specified. [None]')
 
     return parser
Index: trunk/UTIL/PYTHON/planetoplot.py
===================================================================
--- trunk/UTIL/PYTHON/planetoplot.py	(revision 431)
+++ trunk/UTIL/PYTHON/planetoplot.py	(revision 432)
@@ -51,5 +51,5 @@
            tsat=False,\
            flagnolow=False,\
-           movie_axis=0):
+           mrate=None):
 
 
@@ -306,5 +306,7 @@
        varname = all_varname[index_f]
        if varname:
-           what_I_plot, error = reducefield(all_var[index_f], d4=indextime, d1=indexlon, d2=indexlat , d3=indexvert , yint=yintegral, alt=vert, anomaly=anomaly)
+           zindtime=indextime
+           if mrate is not None:zindtime=None
+           what_I_plot, error = reducefield(all_var[index_f], d4=zindtime, d1=indexlon, d2=indexlat , d3=indexvert , yint=yintegral, alt=vert, anomaly=anomaly)
            what_I_plot = what_I_plot*mult
            if not error: 
@@ -334,5 +336,5 @@
                if len(what_I_plot.shape) >= 2:
                  istart=0
-                 if indextime is None:iend=len(time)-1
+                 if mrate is not None:iend=len(time)-1
                  else:iend=istart
                  imov=istart
@@ -340,5 +342,5 @@
                  while imov <= iend:
                     what_I_plot_frame=what_I_plot
-                    if len(what_I_plot.shape) is 3:
+                    if mrate is not None:
                        what_I_plot_frame=what_I_plot[imov,:,:]
                        print "-> frame ",imov+1
@@ -369,14 +371,13 @@
                         mframe=fig2img(mpl.pyplot.gcf())
 
-                        if ((len(what_I_plot.shape) is 3) and (imov is 0)):
+                        if ((mrate is not None) and (imov is 0)):
                             W,H = mpl.pyplot.gcf().canvas.get_width_height()
-                            rate = 10
-                            video = VideoSink((H,W), "test", rate=rate, byteorder="rgba")
-                        if len(what_I_plot.shape) is 3: 
+                            video = VideoSink((H,W), "test", rate=mrate, byteorder="rgba")
+                        if mrate is not None: 
                             video.run(mframe)
                             mpl.pyplot.close()
 
                         imov=imov+1
-                 if len(what_I_plot.shape) is 3: video.close
+                 if mrate is not None: video.close
                ##### 1D field
                elif len(what_I_plot.shape) is 1:
Index: trunk/UTIL/PYTHON/pp.py
===================================================================
--- trunk/UTIL/PYTHON/pp.py	(revision 431)
+++ trunk/UTIL/PYTHON/pp.py	(revision 432)
@@ -137,5 +137,5 @@
                 ope=opt.operat,fileref=reffile,minop=opt.vminope,maxop=opt.vmaxope,titleref=opt.titref,\
                 invert_y=opt.inverty,xaxis=zexaxis,yaxis=zeyaxis,ylog=opt.logy,yintegral=opt.column,\
-                blat=opt.blat,tsat=opt.tsat,flagnolow=opt.nolow)
+                blat=opt.blat,tsat=opt.tsat,flagnolow=opt.nolow,mrate=opt.rate)
         print 'DONE: '+name
         system("rm -f to_be_erased")