def climate_tasks(gdirs): """Prepare the climate data.""" # Only global tasks tasks.distribute_climate_data(gdirs) tasks.compute_ref_t_stars(gdirs) tasks.distribute_t_stars(gdirs)
# gdirs = [gd for gd in gdirs if gd.terminus_type == 'Land-terminating'] # Basic tasks task_list = [ itmix.glacier_masks_itmix, tasks.compute_centerlines, tasks.catchment_area, tasks.initialize_flowlines, tasks.catchment_width_geom, tasks.catchment_width_correction, ] for task in task_list: execute_entity_task(task, gdirs) # Climate related tasks tasks.distribute_climate_data(gdirs) tasks.compute_ref_t_stars(gdirs) tasks.distribute_t_stars(gdirs) # Inversion execute_entity_task(tasks.prepare_for_inversion, gdirs) itmix.optimize_thick(gdirs) execute_entity_task(tasks.volume_inversion, gdirs) # Write out glacier statistics df = utils.glacier_characteristics(gdirs) fpath = os.path.join(cfg.PATHS["working_dir"], "glacier_char.csv") df.to_csv(fpath) if do_itmix:
# For calibration if do_calib: # gdirs = [gd for gd in gdirs if gd.glacier_type != 'Ice cap'] # gdirs = [gd for gd in gdirs if gd.terminus_type == 'Land-terminating'] # Basic tasks task_list = [ itmix.glacier_masks_itmix, tasks.compute_centerlines, tasks.catchment_area, tasks.initialize_flowlines, tasks.catchment_width_geom, tasks.catchment_width_correction ] for task in task_list: execute_entity_task(task, gdirs) # Climate related tasks tasks.distribute_climate_data(gdirs) tasks.compute_ref_t_stars(gdirs) tasks.distribute_t_stars(gdirs) # Inversion execute_entity_task(tasks.prepare_for_inversion, gdirs) itmix.optimize_thick(gdirs) execute_entity_task(tasks.volume_inversion, gdirs) # Write out glacier statistics df = utils.glacier_characteristics(gdirs) fpath = os.path.join(cfg.PATHS['working_dir'], 'glacier_char.csv') df.to_csv(fpath) if do_itmix: