if __name__ == '__main__':
glacier = Glacier()
measured_length = []
linear_equilibrium_length = []
max_surfaces = []
ela = []
critical_ela = []
measured_thickness = []
mean_thickness = []
static_thickness = []
for basin, data in BASINS.items():
glacier.max_bed_height = data['bed_height']
glacier.length = data['length_of_segment']
glacier.slope = calc_slope(data['min_surface'], data['max_surface'],
data['length_of_segment'])
max_surfaces.append(data['max_surface'])
ela.append(glacier.ela_from_length())
critical_ela.append(glacier.critical_ela())
measured_length.append(data['length_of_segment'] / 1000)
linear_equilibrium_length.append(
glacier.linear_equilibrium_length(
ela[-1], thickness=data['mean_ice_thickness']) / 1000)
measured_thickness.append(data['mean_ice_thickness'])
mean_thickness.append(glacier.mean_thickness())
