ax1.legend(title='Slope')
# Change of length over change of ELA
for slope in slopes:
lengths = [
glacier.length_change(ela_change, slope)
for ela_change in ela_changes
]
ax3.plot(ela_changes, lengths, label=slope.round(2))
ax3.set_title('Length change with varying ELA change')
ax3.legend(title='Slope')
# Fraction of change of length over change of ELA
for slope in slopes:
lengths = [
glacier.fractional_length_change(4500, ela_change, slope)
for ela_change in ela_changes
]
ax5.plot(ela_changes, lengths, label=slope.round(2))
ax5.set_title(
'Fractional length change with varying ELA change, ELA = 4500 m')
ax5.legend(title='Slope')
# Change of length over varying slope
for ela_elevation in ela_elevations:
lengths = []
for slope in slopes:
glacier.slope = slope
lengths.append(glacier.linear_equilibrium_length(ela_elevation))
ax2.plot(slopes, lengths, label=ela_elevation)
ax2.set_title('Length change with varying slope')
