os.path.join(cost_dir, 'smoothing_cost_classes.map'))
dem = pcr.readmap(os.path.join(ref_map_dir, 'dem.map'))
# Read input distributions
road_distr = ec.read_distribution(cost_dir, 'roads')
smoothing_distr = ec.read_distribution(cost_dir, 'smoothing')
building_distr = ec.read_distribution(cost_dir, 'building_tot')
dike_raise50_distr = ec.read_distribution(cost_dir, 'dike_raise50')
dike_raise100_distr = ec.read_distribution(cost_dir, 'dike_raise100')
dike_reloc_distr = ec.read_distribution(cost_dir, 'dike_reloc')
land_acq_distr = ec.read_distribution(cost_dir, 'land_acq')
minemb_ext_use_distr = ec.read_distribution(cost_dir, 'minemb_ext_use')
minemb_int_use_distr = ec.read_distribution(cost_dir, 'minemb_int_use')
minemb_polluted_distr = ec.read_distribution(cost_dir, 'minemb_polluted')
groyne = measures.read_dike_maps(ref_map_dir, 'groyne')
minemb = measures.read_dike_maps(ref_map_dir, 'minemb')
main_dike = measures.read_dike_maps(ref_map_dir, 'main_dike')
# Set earthwork values
cost_input_ew = pd.read_csv(os.path.join(cost_dir, 'cost_ew.csv'), index_col=0)
cost_correction = pd.read_csv(os.path.join(cost_dir, 'cost_correction.csv'),
index_col=0,
comment='#')
#%% Instantiate cost classes
### Smoothing
c_sm = ec.CostSmoothing('dummy', smoothing_distr)
### Earth work
c_ew = ec.CostEarthwork('dummy', minemb_ext_use_distr, minemb_int_use_distr,
os.path.join(cost_dir, 'smoothing_cost_classes.map'))
dem = pcr.readmap(os.path.join(ref_map_dir, 'dem.map'))
# Read input distributions
road_distr = read_distribution(cost_dir, 'roads')
smoothing_distr = read_distribution(cost_dir, 'smoothing')
building_distr = read_distribution(cost_dir, 'building_tot')
dike_raise50_distr = read_distribution(cost_dir, 'dike_raise50')
dike_raise100_distr = read_distribution(cost_dir, 'dike_raise100')
dike_reloc_distr = read_distribution(cost_dir, 'dike_reloc')
land_acq_distr = read_distribution(cost_dir, 'land_acq')
minemb_ext_use_distr = read_distribution(cost_dir, 'minemb_ext_use')
minemb_int_use_distr = read_distribution(cost_dir, 'minemb_int_use')
minemb_polluted_distr = read_distribution(cost_dir, 'minemb_polluted')
groyne = read_dike_maps(ref_map_dir, 'groyne')
minemb = read_dike_maps(ref_map_dir, 'minemb')
main_dike = read_dike_maps(ref_map_dir, 'main_dike')
# Set earthwork values
cost_input_ew = pd.read_csv(os.path.join(cost_dir, 'cost_ew.csv'),
index_col=0)
cost_correction = pd.read_csv(os.path.join(cost_dir,
'cost_correction.csv'),
index_col=0,
comment='#')
#%% Test the cost calculation
msr_name = 'sidechannel_everywhere'
label = 'everywhere'
msr_name = 'lowering_everywhere'