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::PDDrandMassBalance::get_PDDs ( double  dt_series,
const std::vector< double > &  S,
const std::vector< double > &  T,
std::vector< double > &  PDDs 
)
virtual

Computes

\[ \text{PDD} = \sum_{i=0}^{N-1} h_{\text{days}} \cdot \text{max}(T_i-T_{\text{threshold}}, 0). \]

Parameters
S\(\sigma\) (standard deviation for daily temperature excursions)
dt_seriestime-series step, in seconds
Tair temperature
Nnumber of points in the temperature time-series, each corresponds to a sub-interval
PDDspointer to a pre-allocated array of length N

Reimplemented from pism::surface::PDDMassBalance.

Definition at line 350 of file localMassBalance.cc.

References pism::k, m_impl, pism::surface::LocalMassBalance::m_seconds_per_day, pism::surface::PDDMassBalance::pdd_threshold_temp, pism::surface::PDDrandMassBalance::Impl::rng, and S().