new test sequence at end of script

da/ct/obspack.py

@@ -10,9 +10,10 @@ File created on 28 Jul 2010.
 """
 import os
 import sys
+sys.path.append(os.getcwd())
+sys.path.append('../../')
 import logging
 import datetime
-sys.path.append(os.getcwd())
 
 identifier = 'CarbonTracker CO2 mixing ratios'
 version = '0.0'
@@ -423,24 +424,41 @@ class MixingRatioList(list):
 
 if __name__ == "__main__":
 
-    from da.tools.initexit import StartLogger
-    from da.tools.pipeline import JobStart
-    from datetime import datetime
     import logging
-    import sys,os
+    from da.tools.initexit import CycleControl
+    from da.ct.dasystem import CtDaSystem 
+    from da.ctgridded.statevector import CtGriddedStateVector 
+
+    sys.path.append('../../')
+
+    logging.root.setLevel(logging.DEBUG)
+
+    DaCycle = CycleControl(args={'rc':'../../dagriddedjet.rc'})
+    DaCycle.Initialize()
+    DaCycle.ParseTimes()
+
+    DaSystem    = CtDaSystem('../rc/carbontracker.rc')
+    DaSystem.Initialize()
+
+    DaCycle.DaSystem    = DaSystem
+
+    obs = ObsPackObservations()
+    obs.DaCycle         = DaCycle
+
+    while DaCycle['time.start'] < DaCycle['time.finish']:
+
+        DaCycle['time.sample.start'] = DaCycle['time.start']
+        DaCycle['time.sample.end'] =  DaCycle['time.end']
 
-    sys.path.append(os.getcwd())
+        obs.Initialize()
+        obs.Validate()
+        obs.AddObs()
+        obs.AddModelDataMismatch()
 
-    StartLogger()
+        print(obs.Data.getvalues('obs'))
+        print(obs.Data.getvalues('mdm'))
 
-    obs = CtObservations()
+        DaCycle.AdvanceCycleTimes()
 
-    DaCycle = JobStart(['-v'],{'rc':'da.rc'})
-    DaCycle['time.sample.start'] = datetime(2000,1,1)
-    DaCycle['time.sample.end'] = datetime(2000,1,2)
 
-    obs.Initialize()
-    obs.Validate()
-    obs.AddObs()
-    print(obs.Data.getvalues('obs'))