19 jday=
iw3jdn(idat(1),idat(2),idat(3))
20 call w3fs26(jday,jy,jm,jd,jdow,jdoy)
function iw3jdn(IYEAR, MONTH, IDAY)
Computes julian day number from year (4 digits), month, and day.
subroutine w3doxdat(idat, jdow, jdoy, jday)
Program history log:
subroutine w3fs26(JLDAYN, IYEAR, MONTH, IDAY, IDAYWK, IDAYYR)
Computes year (4 digits), month, day, day of week, day of year from julian day number.