Skip to content
Snippets Groups Projects
Commit 4343ca4e authored by Peters, Wouter's avatar Peters, Wouter
Browse files

extra analysis added for time averaged file creation

parent f72f98d6
No related branches found
No related tags found
No related merge requests found
...@@ -28,6 +28,7 @@ from da.tm5.observationoperator import TM5ObservationOperator ...@@ -28,6 +28,7 @@ from da.tm5.observationoperator import TM5ObservationOperator
from da.analysis.expand_fluxes import save_weekly_avg_1x1_data, save_weekly_avg_state_data, save_weekly_avg_tc_data, save_weekly_avg_ext_tc_data, save_weekly_avg_agg_data from da.analysis.expand_fluxes import save_weekly_avg_1x1_data, save_weekly_avg_state_data, save_weekly_avg_tc_data, save_weekly_avg_ext_tc_data, save_weekly_avg_agg_data
from da.analysis.expand_molefractions import write_mole_fractions from da.analysis.expand_molefractions import write_mole_fractions
from da.analysis.summarize_obs import summarize_obs from da.analysis.summarize_obs import summarize_obs
from da.analysis.time_avg_fluxes import time_avg
################################################################################################# #################################################################################################
# Parse and validate the command line options, start logging # Parse and validate the command line options, start logging
...@@ -73,14 +74,31 @@ ensemble_smoother_pipeline(dacycle, platform, dasystem, samples, statevector, ob ...@@ -73,14 +74,31 @@ ensemble_smoother_pipeline(dacycle, platform, dasystem, samples, statevector, ob
logging.info(header + "Starting analysis" + footer) logging.info(header + "Starting analysis" + footer)
save_weekly_avg_1x1_data(dacycle, statevector) try:
save_weekly_avg_state_data(dacycle, statevector) save_weekly_avg_1x1_data(dacycle, statevector)
save_weekly_avg_tc_data(dacycle, statevector) save_weekly_avg_state_data(dacycle, statevector)
save_weekly_avg_ext_tc_data(dacycle) except:
save_weekly_avg_agg_data(dacycle,region_aggregate='olson') pass
save_weekly_avg_agg_data(dacycle,region_aggregate='transcom') try:
write_mole_fractions(dacycle) save_weekly_avg_tc_data(dacycle, statevector)
summarize_obs(dacycle) save_weekly_avg_ext_tc_data(dacycle)
except:
pass
try:
save_weekly_avg_agg_data(dacycle,region_aggregate='transcom')
save_weekly_avg_agg_data(dacycle,region_aggregate='olson')
except:
pass
try:
time_avg(dacycle,'transcom')
time_avg(dacycle,'olson')
except:
pass
try:
write_mole_fractions(dacycle)
summarize_obs(dacycle)
except:
pass
sys.exit(0) sys.exit(0)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment