PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800

◆ get_PDDs()

void pism::surface::PDDMassBalance::get_PDDs ( double  dt_series,
const std::vector< double > &  S,
const std::vector< double > &  T,
std::vector< double > &  PDDs 
)
virtual

Compute the expected number of positive degree days from the input temperature time-series.

Use the rectangle method for simplicity.

Parameters
Sstandard deviation for air temperature excursions
dt_serieslength of the step for the time-series
Tair temperature (array of length N)
Nlength of the T array
[out]PDDspointer to a pre-allocated array with N-1 elements

Implements pism::surface::LocalMassBalance.

Reimplemented in pism::surface::PDDrandMassBalance.

Definition at line 115 of file localMassBalance.cc.

References pism::surface::CalovGreveIntegrand(), pism::k, pism::surface::LocalMassBalance::m_seconds_per_day, pdd_threshold_temp, and S().