computeSunPositionDoyHour function

computeSunPositionDoyHour