Commit 4acb2f61 authored by Luijkx, Ingrid's avatar Luijkx, Ingrid
Browse files

changes new structure analysis cteco2

parent 9bcf0114
......@@ -332,15 +332,15 @@ class CycleControl(dict):
strippedname = os.path.split(self['jobrcfilename'])[-1]
self['jobrcfilename'] = os.path.join(self['dir.exec'], strippedname)
shutil.copy(os.path.join(self.dasystem['regionsfile']),os.path.join(self['dir.exec'],'da','analysis','copied_regions.nc'))
shutil.copy(os.path.join(self.dasystem['regionsfile']),os.path.join(self['dir.exec'],'da','analysis','cteco2','copied_regions.nc'))
logging.info('Copied regions file to the analysis directory: %s'%os.path.join(self.dasystem['regionsfile']))
if 'extendedregionsfile' in self.dasystem:
shutil.copy(os.path.join(self.dasystem['extendedregionsfile']),os.path.join(self['dir.exec'],'da','analysis','copied_regions_extended.nc'))
shutil.copy(os.path.join(self.dasystem['extendedregionsfile']),os.path.join(self['dir.exec'],'da','analysis','cteco2','copied_regions_extended.nc'))
logging.info('Copied extended regions file to the analysis directory: %s'%os.path.join(self.dasystem['extendedregionsfile']))
else:
shutil.copy(os.path.join(self['dir.exec'],'da','analysis','olson_extended.nc'),os.path.join(self['dir.exec'],'da','analysis','copied_regions_extended.nc'))
logging.info('Copied extended regions within the analysis directory: %s'%os.path.join(self['dir.exec'],'da','analysis','olson_extended.nc'))
for filename in glob.glob(os.path.join(self['dir.exec'],'da','analysis','*.pickle')):
shutil.copy(os.path.join(self['dir.exec'],'da','analysis','cteco2','olson_extended.nc'),os.path.join(self['dir.exec'],'da','analysis','cteco2','copied_regions_extended.nc'))
logging.info('Copied extended regions within the analysis directory: %s'%os.path.join(self['dir.exec'],'da','analysis','cteco2','olson_extended.nc'))
for filename in glob.glob(os.path.join(self['dir.exec'],'da','analysis','cteco2','*.pickle')):
logging.info('Deleting pickle file %s to make sure the correct regions are used'%os.path.split(filename)[1])
os.remove(filename)
for filename in glob.glob(os.path.join(self['dir.exec'],'*.pickle')):
......
......@@ -142,10 +142,10 @@ def forward_pipeline(dacycle, platform, dasystem, samples, statevector, obsopera
def analysis_pipeline(dacycle, platform, dasystem, samples, statevector):
""" Main entry point for analysis of ctdas results """
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.summarize_obs import summarize_obs
from da.analysis.time_avg_fluxes import time_avg
from da.analysis.cteco2.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.cteco2.expand_molefractions import write_mole_fractions
from da.analysis.cteco2.summarize_obs import summarize_obs
from da.analysis.cteco2.time_avg_fluxes import time_avg
logging.info(header + "Starting analysis" + footer)
......
!!! Info for the CarbonTracker data assimilation system
datadir : /projects/0/ctdas/input/ctdas_2012
datadir : /mnt/beegfs/user/gkoren/ctdas_2012
! XCO2 samples
obs.column.input.dir : ${datadir}/xco2/OCO2_L2_Lite_FP_v10r_10sec/
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment