Hoang Trong Minh Tuan wrote:
> Can anyone help me with this,
> I have a file "data.dat" containing real data with high precision
> like this (2 rows).
>
> 0.596598016919984 -0.382312557470785 -0.528482344762526
> -8.771967130768960E-002 -0.708964641303998 0.120784051948388
>
> I want to read it into an array. The variable is declared as follows
>
> REAL (KIND = SELECTED_REAL_KIND(R=16), DIMENSION(1:2, 1:3) :: INPUTDATA
Here is your problem. You have asked for a REAL type that can represent
numbers between 10**16 and 10**(-16). That is, "R" is for RANGE. You
want PRECISION, which is P=.
Cheers,
--
...........................Malcolm Cohen, Nihon NAG, Tokyo.
|