code to calculate degree days with modified sine method