Skip reading N real numbers while reading data from a ascii file

AI Thread Summary
The discussion focuses on upgrading legacy Fortran 77 code to Fortran 90 standards, specifically regarding the implementation of dynamic arrays for easier parallelization. The main challenge is reading data from a text file while needing to skip a specific number of real numbers that are spread across multiple lines, where the number of real numbers per line is inconsistent. The user seeks an efficient method to handle this, as their current solution is not optimal. A proposed solution involves using a dummy variable to read the skipped numbers into a single variable, allowing the program to ignore them while still advancing through the file. This approach utilizes an implied DO loop, which can also be beneficial for output purposes. The discussion highlights the importance of efficient data handling in the context of upgrading legacy code for modern applications.
prashanth.hsn
Messages
2
Reaction score
0
Hi all,

I'm upgrading a legacy f77 code to f90 standards. I'm trying to make some of the global arrays as dynamic arrays so that at a later stage it would be easier for me to parallelize whole application.
Coming to my problem, I need to read some numbers from a text file so that I can allocate (just enough) memory to the multidimensional arrays. However, while reading the input text file I need to skip say 'n' real numbers and then resume reading. These real numbers span several records (=lines) in the input text file and the number of real numbers per record is not constant sad.
In summary, I know how many real numbers I've to skip, but I do not know in how many records these numbers span.

Ex:
Input.data
Code:
3 3 3
1.23445 0.2345235 0.345234 0.123425235 0.12343245
0.31415927 1.23445 0.2345235 0.345234 0.123425235 0.12343245 
0.31415927 1.23445 0.2345235 
0.31415927 1.23445 0.2345235 0.2345235 
0.123 0.3245 0.74565355 4.34623431 56.2352354
0.31415927 1.23445 0.2345235 0.2345235
3
33
3
33
12
6

I have to read first line and then skip 27 real numbers and then again start reading from line 7.
I currently have a code which works but it is not efficient,
My current code looks like this

Code:
Integer:: i, j, k, ii, jj, kk, num
Real, Allocatable, Dimension(:,:,:):: x 
!Start reading the data from file, UNIT=16
[B]read[/B](16,*) i, j, k
!I want to Skip reading data to dummy x(:,:,:) array
[B]read [/B](16,*) (((x(ii,jj,kk),ii=1,i),jj=1,j),kk=1,k)
! Resume reading
[B]read[/B](16,*) num
 
Technology news on Phys.org
I would declare a dummy variable (e.g. SKIP) and then do
READ ... I , J , K , (SKIP, IJK = 1 , I*J*K ) , ...

In other words, read 27 numbers into the same variable, and then ignore it. That will handle the fact that you don't know how many numbers are on each line of the input file.

Usually, the variables in an implied DO loop are array elements with subscripts, but they don't have to be.

This trick can be useful for output as well, if you want to write N copies of the same value for some reason.
 
dear AlephZero, thanks for your solution. I appreciate your help.
 
Dear Peeps I have posted a few questions about programing on this sectio of the PF forum. I want to ask you veterans how you folks learn program in assembly and about computer architecture for the x86 family. In addition to finish learning C, I am also reading the book From bits to Gates to C and Beyond. In the book, it uses the mini LC3 assembly language. I also have books on assembly programming and computer architecture. The few famous ones i have are Computer Organization and...
I have a quick questions. I am going through a book on C programming on my own. Afterwards, I plan to go through something call data structures and algorithms on my own also in C. I also need to learn C++, Matlab and for personal interest Haskell. For the two topic of data structures and algorithms, I understand there are standard ones across all programming languages. After learning it through C, what would be the biggest issue when trying to implement the same data...

Similar threads

Replies
12
Views
3K
Replies
4
Views
1K
Replies
5
Views
2K
Replies
5
Views
5K
Replies
18
Views
2K
Replies
1
Views
3K
Replies
2
Views
1K
Back
Top