The array elements include zeros in the first few rows and standard scientific notation real elements in Fortran format in other rows. The result of the code I wrote does not give me the exact input data. Negative signs are omitted and the values of the exponents are also ignored.
Find below the code written and the attached file containing the array elements in .txt format to serve as the input data
program test_processing
implicit none
integer :: j
integer :: i
integer, parameter :: n=800
integer, parameter :: m=800
real, dimension(m,n) :: g
real, parameter :: out_unit=6
open ( unit=15,file="recvx_4.txt", status='old', &
access='sequential', form='formatted', action='read', RECL= 1200)
call random_number(g)
do i=1, m
do j = 1, n
read (15, 110) g(i,j)
110 format (800(18X,E14.5E3))
write (6, 120) g
120 format (800(18X, E14.5E3))
open (6,file="matrix.txt",action="WRITE",status="replace", RECL = 1200)
write (6,120) g
PRINT *, g
close (6)
end do
end do
end program test_processing
There are a few areas where I can see that your code might cause a problem ... you have two WRITES and a PRINT ?? How much output are you trying to generate? Did you notice that your WRITE format statement skips 18 columns (18x), and does this 800 times? If you really have 800 wide array, your RECL is not correct (1200). Are you sure you have the correct syntax for reading scientific notation?
Some basic text on FORTRAN input and output will help. If you are an expert and just having some small problems,the areas I mention above might help.
I am not really an expert in FORTRAN programing, but only applied it in my M.Sc programme. However, I am trying to develop myself by applying same for seismic data processing in my current project. Thank you for your observation and suggestions.
- Badges
- Science topic
- Similar topics
- Mathematical Sciences
- Graphs