Bug: Bougamont gives incorrect albedo

[gampnico]

method_Bougamont divides by (t_star * dt_days) to give the wrong albedo. Spotted by @simsegoet.

Expected behaviour

The albedo in Bougamont et al. (2005) should be: $\alpha_{snow} - \frac{(\alpha_{snow} - \alpha_{firn})}{t^{\star}}\Delta t$

[MarcusGastaldello]

Hi,

I'm not sure this is a bug, unless I'm mistaken. I've attached a screenshot from a Jupyter Notebook showing the two albedo methods being equivalent - the original COSIPY one for the Oerlemans & Knapp and the one I coded for Bougamont. Because the decay rate changes you cannot use the first method and have to re-write the equation so that you continually update the albedo value for each timestep rather than recalculating every time with respect to the last snowfall event. This means the formula changes a bit.

As you've commented, the dilemma with this is it slightly ruins the fresh snow system.