# Length change with varying ELA
for slope in slopes:
glacier.slope = slope
length_change = [
glacier.linear_equilibrium_length(elevation)
for elevation in ela_elevations
]
ax1.plot(ela_elevations, length_change, label=slope.round(2))
ax1.set_title('Length change with varying ELA')
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')
