Changeset 527 in lmdz_wrf for trunk/tools

Timestamp:

Jun 22, 2015, 4:07:17 PM (10 years ago)

Author:

lfita

Message:

Initializing statistics matrices at the right place
On `tbackwardSmean' giving back the right matrices

File:

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

trunk/tools/validation_sim.py

@@ -1194 +1194 @@
         dimtf = len(uniqt)
         print '     initially we got',dimt,'values which will become',dimtf
+        print '     ==== NOT READY === '
+        quit(-1)
     else:
         print errormsg
@@ -1398 +1400 @@
             fsoTtvalues[it,1] = soTtvalues[intv[1],0]
 
-        return sovalues, soSvalues, soTtvalues, trjs
+        return fsovalues, fsoSvalues, fsoTtvalues, ftrjs
 
     elif kvals == 'tbackwardOmean':
@@ -1405 +1407 @@
         dimtf = len(uniqt)
         print '     initially we got',dimt,'values which will become',dimtf
-
 
     return 
@@ -1938 +1939 @@
         oFillValue = None
 
+    simstats = np.zeros((Nstsim,5), dtype=np.float)
+    obsstats = np.zeros((Nstsim,5), dtype=np.float)
+    simobsstats = np.zeros((Nstsim,13), dtype=np.float)
+    aroundstats = np.zeros((5,dimt), dtype=np.float)
+
     for kst in range(Nstsim):
         timedn = prestdescsim[kst] + 'time'
@@ -2024 +2030 @@
 
 # statisics sim
-        simstats = np.zeros((Nstsim,5), dtype=np.float)
         simstats[kst,0] = np.min(arrayvals[:,0])
         simstats[kst,1] = np.max(arrayvals[:,0])
@@ -2039 +2044 @@
         obsmask2 = obsmask*obsmask
 
-        obsstats = np.zeros((Nstsim,5), dtype=np.float)
         obsstats[kst,0] = obsmask.min()
         obsstats[kst,1] = obsmask.max()
@@ -2047 +2051 @@
 
 # Statistics sim-obs
-        simobsstats = np.zeros((Nstsim,13), dtype=np.float)
         diffvals = np.zeros((dimt), dtype=np.float)
 
@@ -2102 +2105 @@
 
 # Statistics around sim values
-        aroundstats = np.zeros((5,dimt), dtype=np.float)
         for it in range(dimt):
             aroundstats[0,it] = np.min(simobsSvalues[it,])