oggm_edu.GlacierBed

class oggm_edu.GlacierBed(top=None, bottom=None, width=None, altitudes=None, widths=None, slopes=None, slope_sections=None, nx=200, map_dx=100)

The glacier bed used to construct a oggm_edu.Glacier.

Attributes

bed_harray(float)

Array of the bed surface heights. [m]

bottomint or float

Bottom altitude of the glacier domain. [m]

distance_along_glacierarray(float)

Horisontal distance along the glacier bed. [km]

map_dxint

Grid resolution. [m]

nxint

Number of gridpoints.

topint or float

Top altitude of the glacier domain. [m]

widthint or list(int)

Width of the glacier (single scalar) or list of widths defining the width-profile of the glacier.

widthsarray(float)

Interpolated widths along the distance of the glacier.

__init__(top=None, bottom=None, width=None, altitudes=None, widths=None, slopes=None, slope_sections=None, nx=200, map_dx=100)

Initialise the bed. Pass single scalars for top, bottom and width to create a square glacier bed. For finer control of the width pass altitudes and widths as lists or tuples for custom geometry. This will linearly intepolate between the altitude/width pairs.

To contol the slope of the glacier provide a single value to slopes. For even finer contol pass a sequence of values to slopes along with a sequence of altitudes to slope_sections.

Parameters

topint or float

Elevation at the top of the bed domain in meters.

bottomint or float

Elevation at the bottom of the bed domain in meters.

widthint or float

Width of the bed in meters. Generates a square bed.

altitudesarray_like (ints or floats)

list of values corresponding to the altitude/width distribution. First and last value will be the top and bottom. Length should match widths.

widthsarray_like, (ints of floats)

List of values defining the widths along the glacier. Length should match altitudes.

slopesarray_like(int)

Define the slope of the bed, degrees. One or multiple values. For the latter, slope sections are required. Possible values range from -90 to 90 degrees. To create a hole, there has to be an increase in the altitude in slope_section.

slope_sectionsarray_like(int)

Defines the altitude spans for the slopes. Should start with the top altitude and end with the bottom altitude. The number of altitudes need to be one greater then the number of slopes.

nxint

Number of grid points. Will be overriden if slope is provided.

map_dxint

Grid point spacing in meters.

Methods

__init__([top, bottom, width, altitudes, ...])	Initialise the bed.
plot()	Plot the bed