From fb87545b74b73bbeb4c1061bd2fb41cbfc279469 Mon Sep 17 00:00:00 2001
From: Wouter Peters <wouter.peters@wur.nl>
Date: Mon, 1 Oct 2012 08:02:22 +0000
Subject: [PATCH] now routine ensemble instead of variance in StateToGrid
---
da/baseclasses/statevector.py | 8 +++++---
1 file changed, 5 insertions(+), 3 deletions(-)
diff --git a/da/baseclasses/statevector.py b/da/baseclasses/statevector.py
index 8c7efafa..bfa14d78 100755
--- a/da/baseclasses/statevector.py
+++ b/da/baseclasses/statevector.py
@@ -622,10 +622,12 @@ class StateVector(object):
# And now the covariance, first create covariance matrix (!), and then multiply
deviations = np.array([mem.ParameterValues*fluxvector-ensemblemean for mem in ensemble])
- variance = deviations.std(axis=0)**2
- gridvar = self.VectorToGrid(variance)
+ ensemble = []
+ for mem in deviations:
+ ensemble.append ( self.VectorToGrid(mem) )
- return (gridmean,gridvar,)
+
+ return (gridmean,np.array(ensemble),)
def StateToTC(self, fluxvector = None, lag=1):
"""
--
GitLab