PISM, A Parallel Ice Sheet Model 2.2.2-d6b3a29ca committed by Constantine Khrulev on 2025-03-28
|
◆ CalovGreveIntegrand()Compute the integrand in integral (6) in [CalovGreve05]. The integral is \[\mathrm{PDD} = \int_{t_0}^{t_0+\mathtt{dt}} dt\, \bigg[\frac{\sigma}{\sqrt{2\pi}}\,\exp\left(-\frac{T_{ac}(t)^2}{2\sigma^2}\right) + \frac{T_{ac}(t)}{2}\,\mathrm{erfc} \left(-\frac{T_{ac}(t)}{\sqrt{2}\,\sigma}\right)\bigg] \] This procedure computes the quantity in square brackets. The value \(T_{ac}(t)\) in the above integral is in degrees C. Here we think of the argument This integral is used for the expected number of positive degree days. The user can choose \(\sigma\) by option Definition at line 94 of file localMassBalance.cc. Referenced by pism::surface::PDDMassBalance::get_PDDs(). |