Changeset 519 in lmdz_wrf for trunk/tools

Timestamp:

Jun 19, 2015, 4:11:33 PM (10 years ago)

Author:

lfita

Message:

Right time-bounds

File:

• trunk/tools/validation_sim.py (modified) (2 diffs)

trunk/tools/validation_sim.py

@@ -1393 +1393 @@
                         ftrjs[ip,it] = np.mean(trjs[ip,intv[0]:intv[1]])
                 for iz in range(2*Ng+1):
-                    for iy in rangedd(2*Ng+1):
-                        for ix in rangedd(2*Ng+1):
-                            fsoSvalues[it,iz,iy,ix] = np.meanv(soSvalues[intv[0]:    \
+                    for iy in range(2*Ng+1):
+                        for ix in range(2*Ng+1):
+                            fsoSvalues[it,iz,iy,ix] = np.mean(soSvalues[intv[0]:     \
                               intv[1],iz,iy,ix])
             else:
@@ -1401 +1401 @@
                     for ip in range(3):
                         ftrjs[ip,it] = np.mean(trjs[ip,intv[0]:intv[1]])
-                for iy in rangedd(2*Ng+1):
-                    for ix in rangedd(2*Ng+1):
-                        fsoSvalues[it,iy,ix] = np.meanv(soSvalues[intv[0]:intv[1],   \
+                for iy in range(2*Ng+1):
+                    for ix in range(2*Ng+1):
+                        print iy,ix,':',soSvalues[intv[0]:intv[1],iy,ix]
+                        fsoSvalues[it,iy,ix] = np.mean(soSvalues[intv[0]:intv[1],    \
                           iy,ix])
             print fsoSvalues[it,]
-            print 'soTvalues________________'
-            print soTvalues[intv[0]:intv[1]]
-            print 'soTtvalues________________'
-            print soTtvalues[intv[0]:intv[1],:]
+            print 'soTvalues________________',type(soTvalues)
+            for itt in range(intv[0],intv[1]):
+                print soTvalues[str(itt)]
+            fsoTtvalues[it,0] = soTtvalues[intv[0],0]
+            fsoTtvalues[it,1] = soTtvalues[intv[1],0]
+            print fsoTtvalues[it,:]
             quit()
+
 
     elif kvals == 'tbackwardOmean':