2 ! Contributed by Dominique Dhumieres <dominiq@lps.ens.fr>
4 character(15) :: str="+ .339 567+2"
5 real, parameter :: should_be = .339567e2
6 real, parameter :: eps = 10 * epsilon (should_be)
9 read(str,'(BN,F15.6)') x
14 if (abs (x - should_be) > eps .or. abs (y - should_be) > eps) then