gd4 function

Growing degree days