diff --git a/da/examples/das.py b/da/examples/das.py
index 6da70898b6a9136419ecf65976bd9e4a2f24b91b..73ec734a6379e1548dbddcc8ccba9091b5f61cea 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)