Changeset 610 in lmdz_wrf

Timestamp:

Jul 20, 2015, 1:18:53 PM (10 years ago)

Author:

lfita

Message:

Fixing wrong calculation of statistics
Adding statistics of the values `between'

File:

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

trunk/tools/validation_sim.py

@@ -1937 +1937 @@
 iobst = CFtimes_datetime_NOfile([valdimobs['T'][0]], obstunits, 'standard')
 fobst = CFtimes_datetime_NOfile([valdimobs['T'][dimtobs-1]], obstunits, 'standard')
-print 'Lluis:', simobstimes[0],'tobs:',obstunits
 isimt = CFtimes_datetime_NOfile([simobstimes[0]], obstunits, 'standard')
 fsimt = CFtimes_datetime_NOfile([simobstimes[dimtsim-1]], obstunits, 'standard')
@@ -1963 +1962 @@
 
 exacttvalues = np.array(exacttvalues0, dtype=np.float)
-print 'Lluis: shapes exactvalues:',exacttvalues.shape
 
 if len(exacttvalues) == 0:
@@ -2218 +2216 @@
         arrayvals = np.zeros((len(simobsvalues),2), dtype=np.float)
         for it in range(len(simobsvalues)):
-            arrayvals[it,:] = simobsvalues[it][0:1]
+            arrayvals[it,:] = simobsvalues[it][0:2]
         
         if len(valvars[ivar]) > 2:
@@ -2238 +2236 @@
                 print '  ' + fname + ": operation '"+valvars[ivar][2]+"' not ready!!"
                 quit(-1)
+
+#        for it in range(len(arrayvals[:,0])):
+#            print it,arrayvals[it,:],':',simobsvalues[it]
+#        quit()
 
         if kst == 0:
@@ -2355 +2357 @@
             newvar[:] = arrayvals[:,1]
 
+# Including statistics of the between simulated values
+        if stdescsim[kst] == 'B': 
+            if not searchInlist(onewnc.dimensions, 'Bstats'):
+                newdim = onewnc.createDimension('Bstats',3)
+
+                newvar = onewnc.createVariable('Bstats','c',('Bstats','StrLength'))
+                descvar = prestdescsim[kst] + ' time simulated: ' + valvars[ivar][0]
+                basicvardef(newvar, 'Bstats', 'Between values statistics',           \
+                  ovobs.getncattr('units'))
+                Bstsvals=['minimum','maximum','stdandard deviation'] 
+                writing_str_nc(newvar,Bstsvals,StringLength)
+
+            newvar = onewnc.createVariable(valvars[ivar][0] + '_Bstats_sim', 'f',    \
+              (timedn,'Bstats'), fill_value=fillValueF)
+            descvar = prestdescsim[kst] + ' time simulated: ' + valvars[ivar][0] +   \
+              ' statistics'
+            basicvardef(newvar, valvars[ivar][0], descvar, ovobs.getncattr('units'))
+            newvar[0,:] = [fillValueF, fillValueF, fillValueF]
+            for it in range(dimt-1):
+                newvar[it+1,:] = simobsvalues[it][2:5]
+
 # Around values
         if not onewnc.variables.has_key(valvars[ivar][0] + 'around'):