%%HP: T(3)A(D)F(.);
@ Julian Date calculation routine
@ input : date : in mm.ddyyyy
@ output: julian date

\<< 0 0  
  \-> m d  \<< 
    DUP IP 'm' STO FP 100 * DUP IP 'd' STO FP 10000 * 
    IF m 3 < THEN 1 - END 365.25 * FLOOR m 1 + 
    IF m 2 \<= THEN 12 + END 30.6001 * FLOOR + d + 1720982 + 
  \>> 

\>>
