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