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',
va='bottom',
bbox=dict(boxstyle='round,pad=0.4', fc='wheat', alpha=0.8))
ax1.scatter(x_axis, critical_ela, label='critical ELA', s=100)
ax1.scatter(x_axis, ela, label='ELA from Length', marker='*', s=100)
ax1.scatter(x_axis,
max_surfaces,
label='Max Surface Elevation',
