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

Timestamp:

Jun 15, 2015, 5:38:25 PM (9 years ago)

Author:

lfita

Message:

Using mask pearson correlation

File:

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

trunk/tools/validation_sim.py

@@ -1672 +1672 @@
 
 # Observed values temporally around coincident times
+#    simobsvalues, simobsSvalues, simobsTvalues, simobsTtvalues, trjsim =             \
+#        getting_ValidationValues(obskind, Nexactt, dims, trajpos, ovsim, ovobs,      \
+#        exacttvalues, oFillValue, Ngrid)
+
+# Observed values temporally around coincident times
     simobsvalues, simobsSvalues, simobsTvalues, simobsTtvalues, trjsim =             \
         getting_ValidationValues(obskind, Ncoindt, dims, trajpos, ovsim, ovobs,      \
@@ -1706 +1711 @@
 
 # statisics obs
-    obsmask = ma.masked_equal(arrayvals[:,1], fillValueF)
+# Masking 'nan'
+    obsmask0 = np.where(arrayvals[:,1] != arrayvals[:,1], fillValueF, arrayvals[:,1]) 
+    obsmask = ma.masked_equal(obsmask0, fillValueF)
     obsmask2 = obsmask*obsmask
 
@@ -1715 +1722 @@
     obsstats[3] = obsmask2.mean()
     obsstats[4] = np.sqrt(obsstats[3] - obsstats[2]*obsstats[2])
-
+ 
 # Statistics sim-obs
     simobsstats = np.zeros((9), dtype=np.float)
@@ -1729 +1736 @@
     simobsstats[5] = np.mean(np.abs(diffvals))
     simobsstats[6] = np.sqrt(np.mean(diffvals*diffvals))
-    simobsstats[7], simobsstats[8] = sts.pearsonr(arrayvals[:,0], arrayvals[:,1])
+    simobsstats[7], simobsstats[8] = sts.mstats.pearsonr(arrayvals[:,0], arrayvals[:,1])
 
 # Statistics around sim values