def do_L6_batch(cf_level):
logger = pfp_log.change_logger_filename("pfp_log", "L6")
for i in cf_level.keys():
if not os.path.isfile(cf_level[i]):
msg = " Control file " + cf_level[i] + " not found"
logger.error(msg)
continue
cf_file_name = os.path.split(cf_level[i])
msg = "Starting L6 processing with " + cf_file_name[1]
logger.info(msg)
try:
cf = pfp_io.get_controlfilecontents(cf_level[i])
if "Options" not in cf:
cf["Options"] = {}
cf["Options"]["call_mode"] = "batch"
cf["Options"]["show_plots"] = "No"
infilename = pfp_io.get_infilenamefromcf(cf)
ds5 = pfp_io.nc_read_series(infilename)
ds6 = pfp_levels.l6qc(None, cf, ds5)
outfilename = pfp_io.get_outfilenamefromcf(cf)
ncFile = pfp_io.nc_open_write(outfilename)
pfp_io.nc_write_series(ncFile, ds6)
msg = "Finished L6 processing with " + cf_file_name[1]
logger.info(msg)
logger.info("")
except Exception:
msg = "Error occurred during L6 with " + cf_file_name[1]
logger.error(msg)
error_message = traceback.format_exc()
logger.error(error_message)
continue
return
def do_run_l6(main_gui, cfg=None):
"""
Purpose:
Top level routine for running the L6 gap filling.
Usage:
pfp_top_level.do_run_l6()
Side effects:
Creates an L6 netCDF file with NEE partitioned into GPP and ER.
Author: PRI
Date: Back in the day
Mods:
December 2017: rewrite for use with new GUI
"""
try:
logger.info("Starting L6 processing")
if not cfg:
cfg = pfp_io.load_controlfile(path='controlfiles')
if len(cfg) == 0:
logger.info("Quiting L6 processing (no control file)")
return
in_filepath = pfp_io.get_infilenamefromcf(cfg)
if not pfp_utils.file_exists(in_filepath):
in_filename = os.path.split(in_filepath)
logger.error("File " + in_filename[1] + " not found")
return
ds5 = pfp_io.nc_read_series(in_filepath)
#ds5.globalattributes['controlfile_name'] = cfg['controlfile_name']
sitename = ds5.globalattributes['site_name']
if "Options" not in cfg:
cfg["Options"] = {}
cfg["Options"]["call_mode"] = "interactive"
cfg["Options"]["show_plots"] = "Yes"
ds6 = pfp_levels.l6qc(main_gui, cfg, ds5)
if ds6.returncodes["value"] != 0:
logger.info("Quitting L6: " + sitename)
else:
logger.info("Finished L6: " + sitename)
out_filepath = pfp_io.get_outfilenamefromcf(cfg)
nc_file = pfp_io.nc_open_write(out_filepath)
pfp_io.nc_write_series(nc_file, ds6)
logger.info("Finished saving L6 gap filled data")
logger.info("")
except Exception:
msg = " Error running L6, see below for details ..."
logger.error(msg)
error_message = traceback.format_exc()
logger.error(error_message)
return