Fortran初心者です。サブルーチンを使ったコードを作りたいと思っています。まずサブルーチンの中身をmainプログラムとして書いてみたら、これはうまくいきました。下記がそれです。尚、「array.csv」は単純に1から100までが1次元に書いてあるファイルです。
program readcsv implicit none integer :: n = 100 real pressure integer i real, save :: n_array(100) real direchlet_pressure integer :: ntcabs = 0 do ntcabs = 1, n if(ntcabs==1)then open (17, file='array.csv', status='old') read (17, '()') ! ヘッダ行の読み飛ばし do i = 1, n read (17, *) pressure n_array(i) = pressure !配列のi番目にcsvファイルのi番目を格納 end do direchlet_pressure = n_array(ntcabs) print *,n_array(1) !出力したい配列要素 close (17) else direchlet_pressure = n_array(ntcabs) print *,direchlet_pressure !出力したい配列要素 endif end do end program readcsv
このファイルを次はサブルーチンとして書いたmainプログラムを書きたいとおもっています。具体的にはmainプログラムで「ntcabs」が変化して、その整数を引数としてサブルーチンで同様の計算をしたいのですが、出力が1と0ばかりが1000個ほど出力されてしまいます。私のほしい結果としては1から100まで順番に出力させたいと思っています。下記が現在書いているコードになります。何が問題でしょうか。
program main implicit none integer ntcabs integer :: nstep = 100 do ntcabs = 1, nstep call readcsv(ntcabs) end do end program main !Subroutine !====================================================================================== subroutine readcsv(nstep) implicit none integer nstep real pressure integer i real, save :: n_array(100) real direchlet_pressure ! integer :: n = 100 do nstep = 1, nstep if(nstep==1)then open (17, file='array.csv', status='old') read (17, '()') ! ヘッダ行の読み飛ばし do i = 1, nstep read (17, *) pressure n_array(i) = pressure !配列のi番目にcsvファイルのi番目を格納 end do direchlet_pressure = n_array(1) print *,direchlet_pressure !出力したい配列要素 close (17) else direchlet_pressure = n_array(nstep) print *,direchlet_pressure !出力したい配列要素 endif end do end subroutine readcsv
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。