This function calculates a simple static pressure head in units of meters water equivalent at the base of ice. The function used Bedmap2 data by default, or you can enter your own surface elevations, bed elevations, or ice thickness.
Inspired by: Antarctic Mapping Tools, Antarctic basal geothermal heat flux
Inspired: How to estimate subglacial water routes
Create scripts with code, output, and formatted text in a single executable document.