C
character*20 dum1, dum2, dum3
real t1, t1a, t2, t2a
- real*8 dat1, dat2
- integer i, j, values(8)
+ real*4 dat1, dat2
+ integer i, j, values(8), k
t1 = secnds (0.0)
call date_and_time (dum1, dum2, dum3, values)
t1a = secnds (0.0)
- dat1 = 0.001*real (values(8)) + real (values(7)) +
- & 60.0*real (values(6)) + 3600.0* real (values(5))
- if (((dat1 - t1) < 0.) .or. ((dat1 - t1) > (t1a - t1))) call abort ()
+ dat1 = 0.001 * real(values(8)) + real(values(7)) +
+ & 60.0 * real(values(6)) + 3600.0 * real(values(5))
+ ! handle midnight shift
+ if ((t1a - t1) < -12.0*3600.0 ) t1 = t1 - 24.0*3600.0
+ if ((t1a - dat1) < -12.0*3600.0 ) dat1 = dat1 - 24.0*3600.0
+ if ((dat1 < nearest(t1, -1.)) .or. (dat1 > nearest(t1a, 1.)))
+ & call abort ()
t2a = secnds (t1a)
call date_and_time (dum1, dum2, dum3, values)
t2 = secnds (t1)
- dat2 = 0.001*real (values(8)) + real (values(7)) +
- & 60.0*real (values(6)) + 3600.0* real (values(5))
- if (((dat2 - dat1) < t2a) .or. ((dat2 - dat1) > t2)) call abort ()
+ dat2 = 0.001 * real(values(8)) + real(values(7)) +
+ & 60.0 * real(values(6)) + 3600.0 * real(values(5))
+ ! handle midnight shift
+ if ((dat2 - dat1) < -12.0*3600.0 ) dat1 = dat1 - 24.0*3600.0
+ if (((dat2 - dat1) < t2a - 0.008) .or.
+ & ((dat2 - dat1) > t2 + 0.008)) call abort ()
end