fatiando/rockhound

"""
Bedmap2
=======

Bedmap2 is a suite of gridded products describing surface elevation,
ice-thickness, the sea floor and subglacial bed elevation of the Antarctic
south of 60°S [BEDMAP2]_. Each dataset is projected in Antarctic Polar
Stereographic projection, latitude of true scale -71 degrees south, datum
WGS84. All heights are in metres relative to sea level as defined by the g104c
geoid. The datasets are downloaded as ``tiff`` files and loaded into
a :class:`xarray.Dataset` object.
"""
import rockhound as rh
import matplotlib.pyplot as plt
import cmocean

# Load the ice thickness grid
bedmap = rh.fetch_bedmap2(datasets=["thickness"])
print(bedmap)

plt.figure(figsize=(8, 7))
ax = plt.subplot(111)
pc = bedmap.thickness.plot.pcolormesh(
    ax=ax, cmap=cmocean.cm.ice, cbar_kwargs=dict(pad=0.01, aspect=30)
)
ax.set_title("Bedmap2 Antarctica Ice Thickness")
plt.tight_layout()
plt.show()