debugging expand_fluxes for jet results

da/analysis/expand_fluxes.py

@@ -499,6 +499,7 @@ def SaveWeeklyAvgTCData(DaCycle, StateVector):
             data    = ncf_in.GetVariable(vname)[index]
 
             savedict            = ncf.StandardVar(varname=vname)
+            savedict['name']    = vname
 
             if vname=='latitude': continue
             elif vname=='longitude': continue
@@ -506,7 +507,8 @@ def SaveWeeklyAvgTCData(DaCycle, StateVector):
             elif vname=='idate': continue
             elif 'cov' in vname:
                 
-                cov                 = data.transpose().dot(data) 
+                #cov                 = data.transpose().dot(data) 
+                cov                 = np.dot(data.transpose(),data) 
                 tcdata              = StateVector.VectorToTC(vectordata=cov,cov=True) # vector to TC
 
                 savedict['units']   = '[mol/region/s]**2'
@@ -807,11 +809,11 @@ if __name__ == "__main__":
 
     logging.root.setLevel(logging.DEBUG)
 
-    DaCycle = CycleControl(args={'rc':'../../dagridded.rc'})
+    DaCycle = CycleControl(args={'rc':'../../dajet.rc'})
     DaCycle.Initialize()
     DaCycle.ParseTimes()
 
-    DaSystem    = CtDaSystem('../rc/carbontrackergridded.rc')
+    DaSystem    = CtDaSystem('../rc/carbontrackerjet.rc')
     DaSystem.Initialize()
 
     DaCycle.DaSystem    = DaSystem