oggm.core.massbalance.MultipleFlowlineMassBalance

class oggm.core.massbalance.MultipleFlowlineMassBalance(gdir, fls=None, mu_star=None, mb_model_class=<class 'oggm.core.massbalance.PastMassBalance'>, use_inversion_flowlines=False, input_filesuffix='', bias=None, **kwargs)[source]

Handle mass-balance at the glacier level instead of flowline level.

Convenience class doing not much more than wrapping a list of mass-balance models, one for each flowline.

This is useful for real-case studies, where each flowline might have a different mu*.

Attributes:
fls : list

list of flowline objects

mb_models : list

list of mass-balance objects

__init__(self, gdir, fls=None, mu_star=None, mb_model_class=<class 'oggm.core.massbalance.PastMassBalance'>, use_inversion_flowlines=False, input_filesuffix='', bias=None, **kwargs)[source]

Initialize.

Parameters:
gdir : GlacierDirectory

the glacier directory

mu_star : float or list of floats, optional

set to the alternative value of mu* you want to use (the default is to use the calibrated value). Give a list of values for flowline-specific mu*

fls : list, optional

list of flowline objects to use (defaults to ‘model_flowlines’, and if not available, to ‘inversion_flowlines’)

mb_model_class : class, optional

the mass-balance model to use (e.g. PastMassBalance, ConstantMassBalance…)

use_inversion_flowlines: bool, optional

if True ‘inversion_flowlines’ instead of ‘model_flowlines’ will be used.

input_filesuffix : str

the file suffix of the input climate file

bias : float, optional

set to the alternative value of the calibration bias [mm we yr-1] you want to use (the default is to use the calibrated value) Note that this bias is substracted from the computed MB. Indeed: BIAS = MODEL_MB - REFERENCE_MB.

kwargs : kwargs to pass to mb_model_class

Methods

__init__(self, gdir[, fls, mu_star, …]) Initialize.
get_annual_mb(self, heights[, year, fl_id]) Like self.get_monthly_mb(), but for annual MB.
get_annual_mb_on_flowlines(self[, fls, year]) Get the MB on all points of the glacier at once.
get_ela(self[, year]) Compute the equilibrium line altitude for this year
get_monthly_mb(self, heights[, year, fl_id]) Monthly mass-balance at given altitude(s) for a moment in time.
get_specific_mb(self[, heights, widths, …]) Specific mb for this year and a specific glacier geometry.

Attributes

bias Residual bias to apply to the original series.
prcp_bias Precipitation factor to apply to the original series.
temp_bias Temperature bias to add to the original series.