Changeset 407

Timestamp:

Nov 22, 2011, 11:17:24 AM (13 years ago)

Author:

acolaitis

Message:

Corrected small bug related to all_time

File:

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

trunk/UTIL/PYTHON/planetoplot.py

@@ -225 +225 @@
        print "********** nplot", nplot, "itime",itime,"error",error
        
-       time = all_time[index_f]
+#       time = all_time[index_f]
 
        ### Which local time ?
@@ -250 +250 @@
            m = define_proj(proj,wlon,wlat,back=back,blat=blat)
            x, y = m(lon2d, lat2d)
-
 ####################################################################
        ## get all indexes to be taken into account for this subplot and then reduce field
@@ -257 +256 @@
        indexlat  = getsindex(slat,((nplot-1)//nlon)%nlat,lat)
        indexvert = getsindex(svert,((nplot-1)//(nlon*nlat))%nvert,vert) 
+
+       if fileref is not None:      index_f = ((nplot-1)//(nlon*nlat*nvert*ntime))%(len(namefiles)+2)  ## OK only 1 var, see test in the beginning
+       else:                        yeah = len(namefiles)*len(var) ; index_f = ((nplot-1)//(nlon*nlat*nvert*ntime))%yeah
+       time = all_time[index_f]
+
        if mapmode == 1 and stime is None: 
            indextime = itime
@@ -264 +268 @@
            if typefile in ['mesoapi','meso'] and indextime is not None:  ltst = localtime ( interv[0]+indextime*interv[1], 0.5*(wlon[0]+wlon[1]) )  
 
-       if fileref is not None:      index_f = ((nplot-1)//(nlon*nlat*nvert*ntime))%(len(namefiles)+2)  ## OK only 1 var, see test in the beginning
-       else:                        yeah = len(namefiles)*len(var) ; index_f = ((nplot-1)//(nlon*nlat*nvert*ntime))%yeah
 ####################################################################