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())
static_thickness.append(glacier.static_mean_thickness())
fig, (ax1, ax2, ax3) = plt.subplots(nrows=1, ncols=3, figsize=(14, 6))
x_axis = [1, 2, 3, 4]
labels = dict(textcoords='offset points',
