now routine ensemble instead of variance in StateToGrid

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):
         """