fix for huygens old numpy version

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()