positive degree day for normally distributed temperatures

PDD fraction in a normally distributed climate

Calculates the probability weighted positive degree days for normally distributed temperatures:
fPDD=int(normpdf(T,Tmu,sigma)*T,'T',0,inf)

Usage:
fPDD=normalpdd(Tmu,sigma)

Example:
T_june = 3; %degC
sigma_june = 7; %degC =std(daily_average_temps_june)
PDD_june=30*normalpdd(T_june,sigma_june) % units: degC*days
PDD_factor = 8e-3 %m/day/degC <- different factors exist in litterature.
melt_june = PDD_factor*PDD_june % meters of melt.

reference: Braithwaite 1985.

Aslak Grinsted 2012