dd_dbltri_half function

Compute half-day GDD using triangle functions. This is an internal function, called by dd_calc()