From 13d78e849d6cf70a4f4367ded1aa4981570361bb Mon Sep 17 00:00:00 2001
From: Wouter Peters <wouter.peters@wur.nl>
Date: Fri, 11 May 2012 14:24:26 +0000
Subject: [PATCH] use of obspack

---
 da/examples/das.py | 10 ++++++----
 1 file changed, 6 insertions(+), 4 deletions(-)

diff --git a/da/examples/das.py b/da/examples/das.py
index 6da70898..73ec734a 100755
--- a/da/examples/das.py
+++ b/da/examples/das.py
@@ -41,16 +41,16 @@ from da.tools.pipeline import EnsembleSmootherPipeline
 from da.platform.maunaloa import MaunaloaPlatForm 
 from da.ct.dasystem import CtDaSystem 
 from da.ct.statevector import CtStateVector 
-#from da.ct.obspack import ObsPackObservations 
-from da.ct.obs import CtObservations 
+from da.ct.obspack import ObsPackObservations 
+#from da.ct.obs import CtObservations 
 from da.tm5.observationoperator import TM5ObservationOperator 
 from da.ct.optimizer import CtOptimizer
 
 PlatForm    = MaunaloaPlatForm()
 DaSystem    = CtDaSystem(DaCycle['da.system.rc'])
 ObsOperator = TM5ObservationOperator(DaCycle['da.obsoperator.rc'])
-#Samples     = ObsPackObservations()
-Samples     = CtObservations()
+Samples     = ObsPackObservations()
+#Samples     = CtObservations()
 StateVector = CtStateVector()
 Optimizer   = CtOptimizer()
 
@@ -75,10 +75,12 @@ msg          = header+"Starting analysis"+footer      ; logging.info(msg)
 from da.analysis.expand_fluxes import SaveWeeklyAvg1x1Data
 from da.analysis.expand_fluxes import SaveWeeklyAvgStateData
 from da.analysis.expand_fluxes import SaveWeeklyAvgTCData
+from da.analysis.expand_mixingratios import WriteMixingRatios
 
 savedas      = SaveWeeklyAvg1x1Data(DaCycle, StateVector)
 savedas      = SaveWeeklyAvgStateData(DaCycle, StateVector)
 savedas      = SaveWeeklyAvgTCData(DaCycle, StateVector)
+savedas      = WriteMixingRatios(DaCycle)
 
 sys.exit(0)
 
-- 
GitLab