2 ! Contributed by Dominique Dhumieres <dominiq@lps.ens.fr>
4 character(100) :: str1
= &
5 "123.00456.88 0.123E+01 +0.987+1 -0.2345+02 -0.6879E+2+0.7E+03 0.4E+03"
6 character(100), parameter :: should_be
= &
7 "123.00456.88 0.123E+01 0.987E+01-0.2345E+02-0.6879E+02 0.7E+03 0.4E+03"
8 character(100) :: output
9 complex :: c1
, c2
, c3
, c4
11 100 format ( 2F6
.2
, 2E10
.3
, 2E11
.4
, 2E8
.1
)
12 read (str1
,100) c1
, c2
, c3
, c4
13 write (output
, 100) c1
, c2
, c3
, c4
16 if (output
/= should_be
) then