Changeset 581 in lmdz_wrf for trunk/tools

Timestamp:

Jul 7, 2015, 12:36:59 PM (10 years ago)

Author:

lfita

Message:

Adding what happens when 'Nexactt' == 0

File:

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

trunk/tools/validation_sim.py

@@ -1980 +1980 @@
         print '    ' + prestdescsim[kst] + ' ...'
 
-        if stdescsim[kst] == 'E' and Nexactt != 0:
+        if stdescsim[kst] == 'E':
 # Observed and simualted exact times
             simobsvalues, simobsSvalues, simobsTtvalues, trjsim =                    \
@@ -1992 +1992 @@
                   'simulation', obstunits)
                 set_attribute(newvar, 'calendar', 'standard')
-                newvar[:] = exacttvalues[:,3]
+                if Nexactt == 0:
+                    newvar[:] =  np.float64(0.)
+                    simobsSvalues = np.array(np.float(0.))
+                    simobsvalues = [np.float(0.)]
+                else:
+                    newvar[:] = exacttvalues[:,3]
 
             dimt = Nexactt
@@ -2002 +2007 @@
               ovobs, closesttvalues, oFillValue, Ngrid, 'instantaneous')
             dimt = Nclosest
-
 
             if ivar == 0:
@@ -2157 +2161 @@
         descvar = prestdescsim[kst] + ' time simulated: ' + valvars[ivar][0]
         basicvardef(newvar, valvars[ivar][0], descvar, ovobs.getncattr('units'))
-        newvar[:] = arrayvals[:,0]
+        if stdescsim[kst] == 'E' and Nexactt == 0:
+            newvar[:] = fillValueF
+        else:
+            newvar[:] = arrayvals[:,0]
 
 # obs Values to netCDF
@@ -2169 +2176 @@
         descvar = prestdescsim[kst] + ' time observed: ' + valvars[ivar][1]
         basicvardef(newvar, valvars[ivar][1], descvar, ovobs.getncattr('units'))
-        newvar[:] = arrayvals[:,1]
+
+        if stdescsim[kst] == 'E' and Nexactt == 0:
+            newvar[:] = fillValueF
+        else:
+            newvar[:] = arrayvals[:,1]
 
 # Around values
@@ -2189 +2200 @@
           valvars[ivar][0]
         basicvardef(newvar, vname, descvar, ovobs.getncattr('units'))
-        newvar[:] = simobsSvalues
+
+        if stdescsim[kst] == 'E' and Nexactt == 0:
+            newvar[:] = np.ones((0,Ngrid*2+1,Ngrid*2+1))*fillValueF
+        else:
+            newvar[:] = simobsSvalues
 
 
@@ -2205 +2220 @@
         basicvardef(newvar, vname, descvar, ovobs.getncattr('units'))
         set_attribute(newvar, 'cell_methods', 'time_bnds')
-        newvar[:] = aroundstats.transpose()
+        if stdescsim[kst] == 'E' and Nexactt == 0:
+            newvar[:] = np.ones((0,5))*fillValueF
+        else:
+            newvar[:] = aroundstats.transpose()
 
         if stdescsim[kst] == 'B':