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', ha='right',