Fortran main関数をサブルーチンに書き換えたい

Question

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

Answer 1

もし下記が間違っていたら教えてください。

program main
  implicit none

  integer ntcabs
  integer :: total_nstep = 100

  do ntcabs = 1, total_nstep
  	call readcsv(ntcabs,total_nstep)
  end do
end program main


!Subroutine
!======================================================================================
subroutine readcsv(current_nstep,total_nstep)
  implicit none
  integer i, current_nstep
  integer :: total_nstep

  real pressure, direchlet_pressure
  real, save :: n_array(100)



  do current_nstep = 1, total_nstep
	if(current_nstep==1)then
		open (17, file='array.csv', status='old')
		read (17, '()')       ! ヘッダ行の読み飛ばし
		do i = 1, total_nstep
			read (17, *) pressure
			n_array(i) = pressure   !配列のi番目にcsvファイルのi番目を格納
		end do
		close (17)
		direchlet_pressure = n_array(1)
		print *,direchlet_pressure  !出力したい配列要素
  	
	else
		direchlet_pressure = n_array(current_nstep)
		print *,direchlet_pressure  !出力したい配列要素
	endif
  end do
end subroutine readcsv