computeDayLengthDoy function

computeDayLengthDoy