From 15e2b976d10f0d86d04e74e3a7c623243c94d1e3 Mon Sep 17 00:00:00 2001
From: Wouter Peters <wouter.peters@wur.nl>
Date: Tue, 4 Sep 2012 07:17:55 +0000
Subject: [PATCH] fix for huygens old numpy version

---
 da/analysis/tools_transcom.py | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/da/analysis/tools_transcom.py b/da/analysis/tools_transcom.py
index 590aaeb3..b7ef22f0 100755
--- a/da/analysis/tools_transcom.py
+++ b/da/analysis/tools_transcom.py
@@ -7,7 +7,7 @@ rootdir     = os.getcwd().split('da/')[0]
 analysisdir = os.path.join(rootdir,'da/analysis')
 
 from string import join, split
-from numpy import array, identity, zeros, arange
+from numpy import array, identity, zeros, arange, dot
 import da.tools.io4 as io
 
 # Get masks of different region definitions
@@ -248,7 +248,10 @@ def ExtendedTCRegions(data,cov=False):
     if not cov:
         return dot(array(data).squeeze(),M)
     else:
-        return M.transpose().dot(data).dot(M)
+	try:
+        	return M.transpose().dot(data).dot(M)
+	except:
+        	return dot(dot(M.transpose(),data),M) #Huygens fix
 
 def cov2corr(A):
     b=1./sqrt(A.diagonal())
-- 
GitLab