Fortran Fortran reading scientific form

[Axios]

Hi,

I have been struggling with the Fortran IO for a while. Here is the question, any comments will be much appreciated.

I output some double precision numbers with the format:

write (file_id,'(6(ES17.10E3,2X))'), a,b,c,d,e,f ! say the filename is 'data.dat'

It works fine.

But when I open this file with the command

open (unit = file_id_2, file = "data.dat", status="old")
read (file_id_2,'(6(ES17.10E3,2X))'), a,b,c,d,e,f

I cannot read back what I previously write to the file correctly. Is there any obvious mistake that I made?

Any comments will be appreciated. Thanks!

 

[jtbell]

Can you give some specific examples of numbers that you wrote to the file and what they looked like when you read them in again? (That is, how are they "wrong"?)

 

[Axios]

Can you give some specific examples of numbers that you wrote to the file and what they looked like when you read them in again? (That is, how are they "wrong"?)

The input file looks like this: (with "B" representing the blank, and there are in total 6 number per line)

0.1234567000E+000BB1.0378216950E+000BB1.0378779900E-009BB...

These numbers are output previously by using the command:

write (file_id,'(6(ES17.10E3,2X))'), a,b,c,d,e,f

The problem is, when I use the command

read (file_id_temp,'(6(ES17.10E3,2X))'), a,b,c,d,e,f

(file_id and file_id_temp refers to the same file, I checked this.)

what the code reads are (I just show the first three number corresponding to the one shown above in the input file)

1.0399802540E+009BB1.0656704870E+009BB8.1465474200E+008BB...

It seems that the code reads something, but just not reading correctly...

Any comments will be appreciated, thanks.

--Axios

 

[jtbell]

I'm just guessing wildly here:

1. If the file is being written and read by separate programs, are you sure you're opening the correct file for reading?

2. If the the file is being written and read in the same program: after you finish writing to the file, you need to close it and re-open it so it reads from the beginning again. I also seem to remember a REWIND statement that resets the file to its beginning, but it's been a long time since I used it.

3. When the input data is separated by blanks, you normally don't even need to specify an input format, but can use the default format instead. Try:

read (file_id_temp,*), a,b,c,d,e,f

 

[Axios]

I'm just guessing wildly here:

1. If the file is being written and read by separate programs, are you sure you're opening the correct file for reading?

2. If the the file is being written and read in the same program: after you finish writing to the file, you need to close it and re-open it so it reads from the beginning again. I also seem to remember a REWIND statement that resets the file to its beginning, but it's been a long time since I used it.

3. When the input data is separated by blanks, you normally don't even need to specify an input format, but can use the default format instead. Try:

read (file_id_temp,*), a,b,c,d,e,f

Thanks for the reply.

1. The reason I'm sure I'm reading the correct file is that when I changed the number in the file I meant to read, the output is changed, but changing in a wield way. For example, when I changed the first three number in the first line of input to:

1.1234567000E+000BB1.1378216950E+000BB1.0378779900E-009BB...

the output becomes:

1.0663888460E+009  1.0665093480E+009  8.1465474200E+008

while the original input is (also shown in previous post):

0.1234567000E+000BB1.0378216950E+000BB1.0378779900E-009BB...

and output is:

1.0399802540E+009BB1.0656704870E+009BB8.1465474200E+008BB...

We can see that when the input number changes, the output also changes, in a strange way though...

2. The file to read and the file to write are different, although the file to read was previously generated using the same output format of the same code, but in a different run.

3. Thanks for the comment on using * as the format, I tried this method, and the input and output are:

open (unit = file_id_temp, file = particle_filename, status="old")
read (file_id_temp,*), a,b,c,d,e,f

input:

1.1234567000E+000BB1.1378216950E+000BB1.0378779900E-009BB...

output:

1.0000000000E+000BB1.0000000000E+000BB0.0000000000E+000BB...

So...quite wield...

Thanks for the comments!

--Axios

 

[Axios]

I'm using OSX 10.6, and ifort 11.1, if that helps...thanks.