Changeset 528 in lmdz_wrf for trunk/tools

Timestamp:

Jun 22, 2015, 4:26:14 PM (9 years ago)

Author:

lfita

Message:

Using the right values

File:

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

trunk/tools/validation_sim.py

@@ -14 +14 @@
 import numpy.ma as ma
 
-main = 'validarion_sim.py'
+main = 'validation_sim.py'
 errormsg = 'ERROR -- errror -- ERROR -- error'
 warnmsg = 'WARNING -- warning -- WARNING -- warning'
@@ -1773 +1773 @@
 Nexactt = len(exacttvalues[:,0])
 
-print main + ': found',Nexactt,'Temporal same values in simulation and observations'
+print main + ': found',Nexactt,'Temporal exact values in simulation and observations'
 
 # Sim Closest times
@@ -1939 +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'
@@ -1991 +1986 @@
               getting_ValidationValues(obskind, Ncoindt, dims, trajpos, ovsim, ovobs,\
               coindtvalues, oFillValue, Ngrid, 'tbackwardSmean')
-            dimt = Ncoindt
+            dimt = simobsSvalues.shape[0]
 
             if ivar == 0:
@@ -2006 +2001 @@
                 basicvardef(newvar, vname, 'observed between time', obstunits )
                 set_attribute(newvar, 'calendar', 'standard')
-                newvar[:] = np.unique(coindtvalues[:,2])
+                newvar[:] = np.unique(coindtvalues[:,3])
 
 # Re-arranging values...
@@ -2028 +2023 @@
                 print '  ' + fname + ": operation '"+valvars[ivar][2]+"' not ready!!"
                 quit(-1)
+
+        if kst == 0:
+            simstats = np.zeros((Nstsim,5), dtype=np.float)
+            obsstats = np.zeros((Nstsim,5), dtype=np.float)
+            simobsstats = np.zeros((Nstsim,13), dtype=np.float)
 
 # statisics sim
@@ -2105 +2105 @@
 
 # Statistics around sim values
+        aroundstats = np.zeros((5,dimt), dtype=np.float)
+        print 'Lluis dimt:',dimt,simobsSvalues.shape
         for it in range(dimt):
             aroundstats[0,it] = np.min(simobsSvalues[it,])