Time of sunrise/sunset calculation in sunGeometry.f90 breaks in high and low latitude locations. Fix implements specific checks for the edge cases where -TAN(LAT1)TAN(D) < -1 or > 1, which causes the follow-up equation TD=ACOS(-TAN(LAT1)TAN(D)) to break.
Time of sunrise/sunset calculation in sunGeometry.f90 breaks in high and low latitude locations. Fix implements specific checks for the edge cases where -TAN(LAT1)TAN(D) < -1 or > 1, which causes the follow-up equation TD=ACOS(-TAN(LAT1)TAN(D)) to break.