<jupyter_start><jupyter_text>ETOPO1 Earth Relief===================ETOPO1 is a 1 arc-minute global relief model of Earth's surface that integratesland topography and ocean bathymetry [AmanteEakins2009]_. It's available in twoversions: "Ice Surface" (top of Antarctic and Greenland ice sheets) and"Bedrock" (base of the ice sheets). The grids are loaded into:class:`xarray.Dataset` which can be used to plot and make computations.<jupyter_code>import rockhound as rh
import matplotlib.pyplot as plt
import cmocean

# Load a version of the topography grid
grid = rh.fetch_etopo1(version="bedrock")
print(grid)

# Select a subset that corresponds to Africa to make plotting faster given the
# size of the grid.
africa = grid.sel(latitude=slice(-40, 45), longitude=slice(-20, 60))

# Plot the age grid.
# We're not using a map projection to speed up the plotting but this NOT
# recommended.
plt.figure(figsize=(9, 8))
ax = plt.subplot(111)
africa.bedrock.plot.pcolormesh(
    cmap=cmocean.cm.topo, cbar_kwargs=dict(pad=0.01, aspect=30), ax=ax
)
ax.set_title("ETOPO1")
plt.tight_layout()
plt.show()<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.11219946571682991, 'fraction_numerical': 0.02404274265360641, 'mean_word_length': 5.771084337349397, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 5, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '1', 'dir': 'jupyter-structured-clean-dedup', 'id': '492808', 'license': 'BSD-3-Clause', 'n_tokens_mistral': 397, 'n_tokens_neox': 360, 'n_words': 125, 'path': 'doc/gallery/etopo1.ipynb', 'repo_name': 'fatiando/rockery', 'size': '1123'}