Changeset 511 in lmdz_wrf for trunk/tools/validation_sim.py

Timestamp:

Jun 17, 2015, 5:43:57 PM (9 years ago)

Author:

lfita

Message:

Fixing minor typos/issues

File:

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

trunk/tools/validation_sim.py

@@ -1656 +1656 @@
 
 print main + ': found',Nexactt,'Temporal same values in simulation and observations'
+
 # Sim Coincident times
 ##
+sumsimval = 0.
+sum2simval = 0.
+Nsimt = 0
 coindtvalues0 = []
+coindtvalues0st = []
+tsiminit = 0
+tsimend = 0
+
 for it in range(dimtsim):    
     ot = 0
@@ -1666 +1674 @@
             ot = ito
             tdist = simobstimes[it] - valdimobs['T'][ito]
-            coindtvalues0.append([it, ito, simobstimes[it], valdimobs['T'][ito],      \
+            coindtvalues0.append([it, ito, simobstimes[it], valdimobs['T'][ito],     \
               tdist])
+            Nsimt = Nsimt + 1
+            if tsiminit == 0: tsiminit=simobstimes[it]
+            tsimend = simobstimes[it]
+        elif simobstimes[it] > valdimobs['T'][ito+1]:
+            coindtvalues0st.append([Nsimt, ito, valdimobs['T'][ito],tsimend-tsiminit])
 
 coindtvalues = np.array(coindtvalues0, dtype=np.float)
+coindtvaluesst = np.array(coindtvalues0st, dtype=np.float)
 
 Ncoindt = len(coindtvalues[:,0])
@@ -1688 +1702 @@
 # Dimensions
 newdim = onewnc.createDimension('time',None)
+newdim = onewnc.createDimension('betweentime',None)
 newdim = onewnc.createDimension('bnds',2)
 newdim = onewnc.createDimension('obstime',None)
@@ -1701 +1716 @@
 # Variable dimensions
 ##
-newvar = onewnc.createVariable('obstime','f8',('time'))
-basicvardef(newvar, 'obstime', 'time observations', obstunits )
+newvar = onewnc.createVariable('exacttime','f8',('time'))
+basicvardef(newvar, 'obstime', 'exact coincident time observations and simulation', obstunits)
+set_attribute(newvar, 'calendar', 'standard')
+newvar[:] = exacttvalues[:,3]
+
+newvar = onewnc.createVariable('obstime','f8',('obstime'))
+basicvardef(newvar, 'obstime', 'time observations for between values', obstunits)
 set_attribute(newvar, 'calendar', 'standard')
 set_attribute(newvar, 'bounds', 'time_bnds')
 newvar[:] = coindtvalues[:,3]
+
+newvar = onewnc.createVariable('betweentime','f8',('betweentime'))
+basicvardef(newvar, 'obstime', 'time simulations for between values', simtunits )
+set_attribute(newvar, 'calendar', 'standard')
+set_attribute(newvar, 'bounds', 'time_bnds')
+newvar[:] = coindtvalues[:,2]
 
 newvar = onewnc.createVariable('couple', 'c', ('couple','StrLength'))
@@ -1966 +1992 @@
           aroundostats[3,it])
 
-# sim Values to netCDF
-    newvar = onewnc.createVariable(valvars[ivar][0] + '_sim', 'f', ('time'),         \
+# exact sim Values to netCDF
+    newvar = onewnc.createVariable(valvars[ivar][0] + '_Esim', 'f', ('time'),        \
       fill_value=fillValueF)
-    descvar = 'simulated: ' + valvars[ivar][0]
+    descvar = 'exact time simulated: ' + valvars[ivar][0]
+    basicvardef(newvar, valvars[ivar][0], descvar, ovobs.getncattr('units'))
+    newvar[:] = Earrayvals[:,0]
+
+# exact obs Values to netCDF
+    newvar = onewnc.createVariable(valvars[ivar][1] + '_Eobs', 'f', ('time'),        \
+      fill_value=fillValueF)
+    descvar = 'exact time observed: ' + valvars[ivar][1]
+    basicvardef(newvar, valvars[ivar][1], descvar, ovobs.getncattr('units'))
+    newvar[:] = Earrayvals[:,1]
+
+# between sim Values to netCDF
+    newvar = onewnc.createVariable(valvars[ivar][0] + '_sim', 'f', ('betweentime'),  \
+      fill_value=fillValueF)
+    descvar = 'between observed time 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] + '_obs', 'f', ('time'),         \
+# between obs Values to netCDF
+    newvar = onewnc.createVariable(valvars[ivar][1] + '_obs', 'f', ('obstime'),      \
       fill_value=fillValueF)
     descvar = 'observed: ' + valvars[ivar][1]