Changeset 346 in lmdz_wrf

Timestamp:

Mar 10, 2015, 11:21:17 AM (10 years ago)

Author:

lfita

Message:

Splitting simulated and observed variables
Adding time_bnds

File:

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

trunk/tools/validation_sim.py

@@ -1020 +1020 @@
 gstatsn = ['bias', 'simobs_mean', 'sim_obsmin', 'sim_obsmax', 'sim_obsmean', 'mae',  \
   'rmse', 'r_pearsoncorr', 'p_pearsoncorr']
-ostatsn = ['number of points', 'start', 'end', 'minimum', 'maximum', 'mean',         \
-  'mean2', 'standard deviation']
+ostatsn = ['number of points', 'minimum', 'maximum', 'mean', 'mean2',                \
+  'standard deviation']
 
 parser = OptionParser()
@@ -1333 +1333 @@
 # Dimensions
 newdim = onewnc.createDimension('time',None)
+newdim = onewnc.createDimension('bnds',2)
 newdim = onewnc.createDimension('obstime',None)
 newdim = onewnc.createDimension('couple',2)
@@ -1340 +1341 @@
 newdim = onewnc.createDimension('gstats',9)
 newdim = onewnc.createDimension('stats',5)
-newdim = onewnc.createDimension('tstats',8)
+newdim = onewnc.createDimension('tstats',6)
 
 # Variable dimensions
@@ -1347 +1348 @@
 basicvardef(newvar, 'obstime', 'time observations', obstunits )
 set_attribute(newvar, 'calendar', 'standard')
+set_attribute(newvar, 'bounds', 'time_bnds')
 newvar[:] = coindtvalues[:,3]
 
@@ -1608 +1610 @@
 
 # Statistics around obs values
-    aroundostats = np.zeros((8,Ncoindt), dtype=np.float)
+    aroundostats = np.zeros((6,Ncoindt), dtype=np.float)
 
     for it in range(Ncoindt):
@@ -1615 +1617 @@
 
         aroundostats[0,it] = len(obsmask.flatten())
-        aroundostats[1,it] = simobsTtvalues[it,0]
-        aroundostats[2,it] = simobsTtvalues[it,1]
-        aroundostats[3,it] = obsmask.min()
-        aroundostats[4,it] = obsmask.max()
-        aroundostats[5,it] = obsmask.mean()
-        aroundostats[6,it] = obsmask2.mean()
-        aroundostats[7,it] = np.sqrt(aroundostats[6,it] - aroundostats[5,it]*        \
-          aroundostats[5,it])
-
-# Values to netCDF
-    newvar = onewnc.createVariable(varsimobs, 'f', ('time','couple'),                \
+        aroundostats[1,it] = obsmask.min()
+        aroundostats[2,it] = obsmask.max()
+        aroundostats[3,it] = obsmask.mean()
+        aroundostats[4,it] = obsmask2.mean()
+        aroundostats[5,it] = np.sqrt(aroundostats[4,it] - aroundostats[3,it]*        \
+          aroundostats[3,it])
+
+# sim Values to netCDF
+    newvar = onewnc.createVariable(valvars[ivar][0], 'f', ('time'),                  \
       fill_value=fillValueF)
-    descvar = 'couples of simulated: ' + valvars[ivar][0] + ' and observed ' +       \
-      valvars[ivar][1]
-    basicvardef(newvar, varsimobs, descvar, ovobs.getncattr('units'))
-    newvar[:] = arrayvals
+    descvar = 'simulated: ' + valvars[ivar][0]
+    basicvardef(newvar, valvars[ivar][0], descvar, ovobs.getncattr('units'))
+    newvar[:] = arrayvals[:,0]
+
+# obs Values to netCDF
+    newvar = onewnc.createVariable(valvars[ivar][1], 'f', ('time'),                  \
+      fill_value=fillValueF)
+    descvar = 'observed: ' + valvars[ivar][1]
+    basicvardef(newvar, valvars[ivar][1], descvar, ovobs.getncattr('units'))
+    newvar[:] = arrayvals[:,1]
 
 # Around values
@@ -1682 +1688 @@
         newvar[:] = aroundstats.transpose()
 
+    if not searchInlist(onewnc.variables, 'time_bnds'):
+        newvar = onewnc.createVariable('time_bnds','f8',('time','bnds'))
+        basicvardef(newvar, 'time_bnds', 'time', obstunits )
+        set_attribute(newvar, 'calendar', 'standard')
+        newvar[:] = simobsTtvalues
+
 # around obs Statistics
     if not searchInlist(onewnc.variables,valvars[ivar][1] + 'staround'):
@@ -1689 +1701 @@
         basicvardef(newvar, valvars[ivar][1] + 'staround', descvar,                  \
           ovobs.getncattr('units'))
+        set_attribute(newvar, 'cell_methods', 'statistics')
+
         newvar[:] = aroundostats.transpose()