Changeset 584 in lmdz_wrf

Timestamp:

Jul 8, 2015, 1:06:14 PM (10 years ago)

Author:

lfita

Message:

Fixing non-comprensible error!?

File:

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

trunk/tools/validation_sim.py

@@ -1401 +1401 @@
                 stdv = np.std(sovs)
 
-                fsovalues.append([meanv, np.array(sovalues[intv[0]:intv[1]])[0,1],    \
+                fsovalues.append([meanv, np.array(sovalues[intv[0]:intv[1]])[0,1],   \
                   minv, maxv, stdv])
             else:
-                fsovalues.append([fillValueF, np.ones((2,2),dtype=np.float)*fillValueF\
-                  , fillValueF, fillValueF, fillValueF]) 
-
+                fsovalues.append([fillValueF, fillValueF, fillValueF, fillValueF,    \
+                  fillValueF, fillValueF]) 
             if ds.has_key('Z'):
                 if okind == 'trajectory':
@@ -2053 +2052 @@
 
 # Re-arranging values...
-        arrayvals = np.array(simobsvalues)
+# For an incomprensible reason it is not working?
+#        arrayvals = np.array(simobsvalues)
+        arrayvals = np.zeros((len(simobsvalues),2), dtype=np.float)
+        for it in range(len(simobsvalues)):
+            arrayvals[it,:] = simobsvalues[it][0:1]
+        
         if len(valvars[ivar]) > 2:
             const=np.float(valvars[ivar][3])
@@ -2079 +2083 @@
 
 # statisics sim
-        print 'Lluis shapes: ',arrayvals.shape,':',len(simobsvalues)
         simstats[kst,0] = np.min(arrayvals[:,0])
         simstats[kst,1] = np.max(arrayvals[:,0])