3 ! The IO library has an algorithm that involved repeated multiplication by 10,
4 ! resulting in introducing large cumulative floating point errors.
8 s = "-.18774312893273 "
9 read(unit=s, fmt='(g20.14)') d
10 if (d + 0.18774312893273d0 .gt. 1d-13) call abort