Changeset 468 in lmdz_wrf for trunk/tools

Timestamp:

Jun 9, 2015, 6:15:34 PM (10 years ago)

Author:

lfita

Message:

'Final'? correction for the `None' vardim for the simulation
Solved issue with the masked values from the observations

File:

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

trunk/tools/validation_sim.py

@@ -1422 +1422 @@
 
     ovobs = oobs.variables[valvars[ivar][1]]
+    if searchInlist(ovobs.ncattrs(),'_FillValue'): 
+        oFillValue = ovobs.getncattr('_FillValue')
 
 # Simulated values spatially around coincident times
@@ -1451 +1453 @@
         for it in range(Ncoindt):
             ito = int(coindtvalues[it,1])
-            slicev = dims['X'][0] + ':' + str(stationpos[1]) + '|' +                 \
-              dims['Y'][0] + ':' + str(stationpos[0]) + '|' +                        \
-              dims['T'][0] + ':' + str(int(coindtvalues[it][0]))
+            if valdimsim.has_key('X') and valdimsim.has_key('Y'):
+                slicev = dims['X'][0] + ':' + str(stationpos[1]) + '|' +             \
+                  dims['Y'][0] + ':' + str(stationpos[0]) + '|' +                    \
+                  dims['T'][0] + ':' + str(int(coindtvalues[it][0]))
+            else:
+                slicev = dims['T'][0] + ':' + str(int(coindtvalues[it][0]))
             slicevar, dimslice = slice_variable(ovsim, slicev)
-            if ovobs[int(ito)] != ovobs[int(ito)]: 
+            if ovobs[int(ito)] == oFillValue or ovobs[int(ito)] == '--':
                 simobsvalues.append([ slicevar, fillValueF])
+#            elif ovobs[int(ito)] != ovobs[int(ito)]: 
+#                simobsvalues.append([ slicevar, fillValueF])
             else:
                 simobsvalues.append([ slicevar, ovobs[int(ito)]])
-            slicev = dims['X'][0] + ':' + str(stationpos[1]-Ngrid) + '@' +           \
-              str(stationpos[1]+Ngrid+1) + '|' + dims['Y'][0] + ':' +                \
-              str(stationpos[0]-Ngrid) + '@' + str(stationpos[0]+Ngrid+1) + '|' +    \
-              dims['T'][0] + ':' + str(int(coindtvalues[it,0]))
+            if valdimsim.has_key('X') and valdimsim.has_key('Y'):
+                slicev = dims['X'][0] + ':' + str(stationpos[1]-Ngrid) + '@' +       \
+                  str(stationpos[1]+Ngrid+1) + '|' + dims['Y'][0] + ':' +            \
+                  str(stationpos[0]-Ngrid) + '@' + str(stationpos[0]+Ngrid+1) + '|' +\
+                  dims['T'][0] + ':' + str(int(coindtvalues[it,0]))
+            else:
+                slicev = dims['T'][0] + ':' + str(int(coindtvalues[it][0]))
             slicevar, dimslice = slice_variable(ovsim, slicev)
             simobsSvalues[it,:,:] = slicevar
@@ -1629 +1639 @@
 
 # sim Values to netCDF
-    newvar = onewnc.createVariable(valvars[ivar][0], 'f', ('time'),                  \
+    newvar = onewnc.createVariable(valvars[ivar][0] + '_sim', 'f', ('time'),         \
       fill_value=fillValueF)
     descvar = 'simulated: ' + valvars[ivar][0]
@@ -1636 +1646 @@
 
 # obs Values to netCDF
-    newvar = onewnc.createVariable(valvars[ivar][1], 'f', ('time'),                  \
+    newvar = onewnc.createVariable(valvars[ivar][1] + '_obs', 'f', ('time'),         \
       fill_value=fillValueF)
     descvar = 'observed: ' + valvars[ivar][1]