52 integer :: yyyymmdd,yyyy,mm,dd,hh,mi,ss,hhmiss
53 integer :: julday,jy,jm,ja
54 integer,
parameter :: igreg=15+31*(10+12*1582)
58 mm=(yyyymmdd-10000*yyyy)/100
59 dd=yyyymmdd-10000*yyyy-100*mm
61 mi=(hhmiss-10000*hh)/100
62 ss=hhmiss-10000*hh-100*mi
65 print*,
'there is no year zero.'
68 if (yyyy.lt.0) yyyy=yyyy+1
76 julday=int(365.25*jy)+int(30.6001*jm)+dd+1720995
77 if (dd+31*(mm+12*yyyy).ge.igreg)
then
79 julday=julday+2-ja+int(0.25*ja)
82 juldate=
real(julday,kind=dp) +
real(hh,kind=dp)/24._dp + &
83 real(mi,kind=dp)/1440._dp +
real(ss,kind=dp)/86400._dp
real(kind=dp) function juldate(yyyymmdd, hhmiss)