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