Compute half-day GDD using triangle functions. This is an internal function, called by dd_calc()
dd_dbltri_half(case, tmin, tmax, thresh_low, thresh_up)