テキストファイルからデータを読み取り、パーツ/変数を別のテキストファイルに保存/書き込むにはどうすればよいですか？

Question

私はFortranをかなり新しくしているので、これは素朴な質問かもしれません。私は行の＃540001と巨大な.txtファイルを読んでみたいと思います。データは、5つの高さでu、v、w、Tとして速度および温度の測定を伴う音響風速計からのものです。

"2011-07-10 09:30:00",9838,1.132,2.30225,-0.5635,29.18585,0.30275,0.689,-0.01125,29.67004,0.2165,-0.25475,0.12725,29.8923,0.51425,3.0405,-0.58375,29.5242,-0.0085,3.6235,-0.65175,29.61972,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 
"2011-07-10 09:30:00.05",9839,-0.21325,3.22775,-0.17,29.10953,0.33925,0.6867501,-0.0015,29.67874,0.1715,-0.196,0.1235,29.8923,0.035,2.6915,-0.3845,29.82806,-0.102,3.5505,-0.15825,29.61795,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 
"2011-07-10 09:30:00.1",9840,0.403,3.1195,-0.37175,29.22574,0.06550001,0.6655,0.1275,29.76208,0.196,-0.2,0.1,29.901,0.16225,2.31525,-0.5975,29.69263,0.24175,3.11925,-0.3725,29.57977,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 

など、私は保存/行列U（5,540001）、V（5,540001）を書き、（Wたい

...：ここでは、ファイルの最初の数行があります5,540001）、およびT（5,540001）を使用して、計算とプロットを行うことができます。ここで私は使用していますものです：

PROGRAM READ_MAIZE 
    IMPLICIT NONE 

    REAL,DIMENSION(:,:),Allocatable :: u_r, v_r, w_r, T_r 

    CHARACTER(len=*) :: fname 
    fname='FINALDATA.txt' 
    open(unit=1,file=fname,status='old',action='read') 

    do i=1,540001 
    READ(1,*)timestamp(i),count(i),u_r(5,i),v_r(5,i),w_r(5,i), & 
    T_r(5,i),u_r(2,i),v_r(2,i),w_r(2,i),T_r(2,i), & 
u_r(1,i),v_r(1,i),w_r(1,i),T_r(1,i), & 
    u_r(3,i),v_r(3,i),w_r(3,i),T_r(3,i), & 
u_r(4,i),v_r(4,i),w_r(4,i),T_r(4,i),flags(1:20) 
    end do 
    close(1) 

    WRITE(U_maize,'(A,I7.7,A,I7.7,A)'), & 
     '.txt' 
    open(11,file=U_maize,status='unknown',action='write') 
    write(11,'(F20.14)')(u_r) 
    end 

がu_r順番を気にしない（5、I）u_rが続く（2、i）は...（彼らは秩序の外にある異なる高さに相当します）。これは動作していません。

Answer 1

あなたのコードにはかなりのことがあり、最初にやろうとしていることを理解することが難しくなります。私はあなたのコードに以下の注釈をつけ、コンパイルして出力を生成するものに変換しました。多分それが助けになるでしょう。

PROGRAM READ_MAIZE 
    IMPLICIT NONE ! This means that every variable has to be declared 
        ! explicitly. You don't do that. So I did it for you 

    REAL,DIMENSION(:,:),Allocatable :: u_r, v_r, w_r, T_r 
    integer, dimension(:), allocatable :: data_count ! using fortran keywords 
      ! (such as "count") as variables is dangerous and should be avoided 
    character(len=len("2011-07-10 09:30:00.05")), allocatable :: timestamp(:) 

    CHARACTER(len=*), parameter :: fname='FINALDATA.txt' 
    character(len=len("U_XXXXXXX_XXXXXXX.txt")) :: U_maize 
    integer :: in_unit, out_unit ! Use integer variables for the unit. 
            ! together with newunit keyword, this is 
            ! safer 
    integer, parameter :: num_records = 3 ! you need to up this number to 
              ! 540001 again 
    integer :: i 

    ! If you have allocatable arrays, you need to allocate them before you 
    ! can use them 
    allocate(u_r(5, num_records)) 
    allocate(v_r(5, num_records)) 
    allocate(w_r(5, num_records)) 
    allocate(T_r(5, num_records)) 
    allocate(data_count(num_records)) 
    allocate(timestamp(num_records)) 

    ! the "newunit" keyword is a safe way to create a unique unit 
    ! identifier. You should really use this. 
    open(newunit=in_unit,file=fname,status='old',action='read') 

    do i=1,num_records 
     READ(in_unit,*) timestamp(i), data_count(i), & 
      u_r(5,i),v_r(5,i),w_r(5,i),T_r(5,i),  & 
      u_r(2,i),v_r(2,i),w_r(2,i),T_r(2,i),  & 
      u_r(1,i),v_r(1,i),w_r(1,i),T_r(1,i),  & 
      u_r(3,i),v_r(3,i),w_r(3,i),T_r(3,i),  & 
      u_r(4,i),v_r(4,i),w_r(4,i),T_r(4,i)  ! I left out the flags 
                ! since I didn't know what 
                ! that was. 
    end do 
    close(in_unit) 

    ! I don't know how the file name should be constructed, except 
    ! that it should end in a .txt, and the format. So I made something up. 
    write(U_maize, '(A, I7.7, A, I7.7, A)') 'U_', 35, '_', 6, '.txt' 

    open(newunit=out_unit,file=U_maize,status='unknown',action='write') 

    ! To make it more readable, I tell write to write 5 numbers per row, 
    ! Not sure whether this is what you want. 
    write(out_unit,'(5(X, F20.14))') u_r 
    close(out_unit) ! I know it isn't technically needed, but please always 
     ! close files when finished, even if the program terminates anyway. 
end program READ_MAIZE ! tell the compiler what you want to end here.