Description: Gregorian Calendar to Julian Day Number.
The source code for the Fortran routine CAL2JD is available: