Recent content by melmacianalf

How about program main implicit none character*50 what read (*,'(a50)'), what write (*,'(a50)'), what end program Gordon Shumway

ferry2: May be you can post a sample 'Input_data.dat' ? Gordon Shumway

I see nothing wrong in it. It is very common for someone who consults a forum for the first time to post the same question in many forums. Gordon Shumway

In cases like this, I have been using 'sed' (thanks to mikrom) to delete the lines containing '.' (in real numbers). The command is sed '/\./d' old_file > new_file You can use it in F90 as follows program using_sed implicit none character(*), parameter :: sed_cmd = "sed '/\./d'...

mikrom: Many thanks for the link to the previous post. Gordon Shumway

Tony1984: Do you still want to do it with fortran? Gordon Shumway

You can try something like this program extract implicit none integer, parameter :: skip_initial_lines = 52 integer, parameter :: print_middle_lines = 471 integer :: i character(len=100) :: line character(len=10) :: c1, c2, c3 double precision :: r1, r2, r3...

Tony1984: I think you ahould use other languages like 'sed' to do this task. In linux try sed -ne '/^DATA 3 467/,/^ENDDATA 469/p' file.data where file.data is the file containing your data. Gordon Shumway

First of all, people working with large codes use some version control programs. I use darcs to update F90 codes that have millions of lines. When I update a module or any part of the code, I only pass the 'update' or 'patches' that contains new lines. Using darcs my colleague can selectively...

program min_xyz implicit none...

You are right, I missed it. This is most likely the source of error. Gordon Shumway

stenssok: It is still not clear where the error lies. can you also post the full input from "combin.fig" In your first post you said that the error lies in the following statement READ(3,*) ((XPOS(J,IL),YPOS(J,IL),J=1,NABUR(IL)),IL=1,NL) Did u arrive at this conclusion b'coz only at...

stenssok: From all the informations that you had given, I am able to read and print the data using f77 without any problem. I see the following options 1. Add the statement REWIND(3) just above the read statement and try again 2. Perhaps you can try to post the full code or more code...

what are the values of NABURMAX and NLMAX? Gordon Shumway

Did you declare XPOS and YPOS as real data? Here is what I did program main implicit none integer :: J, IL integer, parameter :: NABUR(4) = (/8, 8, 0, 0/) integer, parameter :: NL = 2 real :: XPOS(8,2), YPOS(8,2) open(unit=3...