From 2ee20c3ca1dd72b97fd7f1651c07cc5758748c5e Mon Sep 17 00:00:00 2001 From: Wouter Peters <wouter.peters@wur.nl> Date: Tue, 4 Sep 2012 07:17:35 +0000 Subject: [PATCH] fix for huygens old numpy version --- da/analysis/expand_fluxes.py | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/da/analysis/expand_fluxes.py b/da/analysis/expand_fluxes.py index 132dbbf..0abd2ea 100755 --- a/da/analysis/expand_fluxes.py +++ b/da/analysis/expand_fluxes.py @@ -541,7 +541,10 @@ def SaveWeeklyAvgTCData(DaCycle, StateVector): tcdata.append( StateVector.VectorToTC(vectordata=member) ) tcdata = np.array(tcdata) - cov = tcdata.transpose().dot(tcdata)/(StateVector.nmembers-1) + try: + cov = tcdata.transpose().dot(tcdata)/(StateVector.nmembers-1) + except: + cov = np.dot(tcdata.transpose(),tcdata)/(StateVector.nmembers-1) # Huygens fix #print vname,cov.sum() -- GitLab