Bedmachine Grids

Description: Gridded data in EPSG 3431 projection. Topographic layers, uncertainties, and masks.

Requires an EarthData login, see link for how to configure this.

Dataset:

National Snow and Ice Data Center. “IceBridge BedMachine Greenland, Version 5,” October 21, 2020. https://nsidc.org/data/idbmg4/versions/5.

Associated paper:

N.A.

%%capture
%load_ext autoreload
%autoreload 2

import pygmt

from polartoolkit import fetch, maps, regions, utils

layer_names = [
    # topographies
    "surface",
    "icebase",
    "bed",
    # uncertainties
    "errbed",
    # others
    "ice_thickness",
    "dataid",
    "geoid",
    "mask",
    "source",
]

grids = []
for name in layer_names:
    data = fetch.bedmachine(
        layer=name,
        hemisphere="north",
        # available options
        # reference ("eigen-6c4" or "ellipsoid")
        # region,
        # spacing,
        # registration ("g" for gridline or "p" for pixel),
    )
    grids.append(data)
    print(f"Info for {name}")
    _ = utils.get_grid_info(data, print_info=True)
    print("##########")

Info for surface
grid spacing: 150.0 m
grid region: (-653000.0, 879700.0, -3384350.0, -632750.0)
grid zmin: 0.0
grid zmax: 3673.38549805
grid registration: p
##########
Info for icebase
grid spacing: 150.0 m
grid region: (-653000.0, 879700.0, -3384350.0, -632750.0)
grid zmin: -1913.28369141
grid zmax: 3673.34838867
grid registration: p
##########
Info for bed
grid spacing: 150.0 m
grid region: (-653000.0, 879700.0, -3384350.0, -632750.0)
grid zmin: -5571.67285156
grid zmax: 3673.34838867
grid registration: p
##########
Info for errbed
grid spacing: 150.0 m
grid region: (-653000.0, 879700.0, -3384350.0, -632750.0)
grid zmin: 0.0
grid zmax: 1000.0
grid registration: p
##########
Info for ice_thickness
grid spacing: 150.0 m
grid region: (-653000.0, 879700.0, -3384350.0, -632750.0)
grid zmin: 0.0
grid zmax: 3409.73779297
grid registration: p
##########
Info for dataid
grid spacing: 150.0 m
grid region: (-653000.0, 879700.0, -3384350.0, -632750.0)
grid zmin: 0.0
grid zmax: 10.0
grid registration: p
##########
Info for geoid
grid spacing: 150.0 m
grid region: (-653000.0, 879700.0, -3384350.0, -632750.0)
grid zmin: 6.0
grid zmax: 64.0
grid registration: p
##########
Info for mask
grid spacing: 150.0 m
grid region: (-653000.0, 879700.0, -3384350.0, -632750.0)
grid zmin: 0.0
grid zmax: 4.0
grid registration: p
##########
Info for source
grid spacing: 150.0 m
grid region: (-653000.0, 879700.0, -3384350.0, -632750.0)
grid zmin: 1.0
grid zmax: 50.0
grid registration: p
##########

subplot_dimensions = utils.square_subplots(len(layer_names))

region = regions.greenland
_, _, fig_width, fig_height = utils.set_proj(
    region,
    hemisphere="north",
    fig_height=15,
)

fig = pygmt.Figure()

with fig.subplot(
    nrows=subplot_dimensions[0],
    ncols=subplot_dimensions[1],
    subsize=(fig_width, fig_height),
    frame="f",
    margins="1.5c",
):
    for i, (layer, grid) in enumerate(zip(layer_names, grids)):
        with fig.set_panel(panel=i):
            # plot the grids
            maps.plot_grd(
                grid,
                fig=fig,
                title=layer,
                origin_shift=None,
                region=region,
                coast=True,
                robust=True,
                hemisphere="north",
            )
fig.show(dpi=200)

The Kernel crashed while executing code in the current cell or a previous cell.

Please review the code in the cell(s) to identify a possible cause of the failure.

Click <a href='https://aka.ms/vscodeJupyterKernelCrash'>here</a> for more info.

View Jupyter <a href='command:jupyter.viewOutput'>log</a> for further details.