Alpine paleoglaciers

Plot Ehlers et al. (2011) global paleoglacier data on top of gridded Alpine topographic data, whose coordinate system is read from an netCDF attribute.

downloading http://static.us.elsevierhealth.com/ehlers_digital_maps/digital_maps_02_all_other_files.zip...

import matplotlib.pyplot as plt
import hyoga

# open demo data
with hyoga.open.example('pism.alps.in.boot.nc') as ds:

    # plot model output
    ds.hyoga.plot.bedrock_altitude(vmin=0)
    ax = ds.hyoga.plot.paleoglaciers(alpha=0.75)

# set title
ax.set_title('Last Glacial Maximum (Ehlers et al., 2011)')

# show
plt.show()