Hello,
I am trying to read a formatted 2D array from a file on disk into a variable. I have the write operation, which is rather simple, but am stuck on reading the same file. Could someone point me to a sample/writeup on how to do this? The net seems saturated, but i can't find a useful article.
Btw, the reason for the formatted file is to keep it human readable as it contains configuration options.
Thank you in advance.
I've actually found that the physical documentation that comes with the compiler is generally the most readable and most informative for Fortran compilers. Of course, that's not an option if you're using g95 or something like that.
Here's a pretty good page describing most of the technical specs of the read statement. Particularly, see the section on "Format Edit Descriptors" - very handy.
On a side note, if you have the exact write format string, you can usually drop that into a read format string, but if you're writing with WRITE(*,*) or something like that, you probably won't have a valid write format statement to use.
Finally, if you're dumping this out to ASCII so people can read it, and you don't have to worry about backward compatibility, consider dumping everything out as fixed-length fields, as they are by far the easiest things to read back in.
Sorry I can't think of better online resources, but Fortran is woefully underdocumented on the web. I remember once checking to see if g95 had Fortran reference docs, but they mostly only have docs on their specific compiler settings. Good luck, though!
Could the OP post the code finally used? I'm struggling with the same problem and cannot find an answer anywhere.
In particular I wonder if one can read line by line using something like
do line=1,nlines read(fid,*) myarray(line,:) enddo
or if one has to loop over all the numbers on each line too, using non-advancing reads.
Any help is appreciated.
TD
PS: Sorry for posting a new "answer". I'm new here and couldn't find out how to add to the comments to the answer by Mike.