データ数が未知のファイルを読み込む

program test_read

character(len=500)::infle

real,allocatable::x(:) 

character(len=2000)::strm

namelist /para/infle

read(*,nml=para)

open(11,file=infle,action="read")

n=0

count_vaild_data: do

  read(11,'(A)',iostat=ios)strm

  if(ios<0)exit

  if(strm(1:1) == "#")then

    cycle

  else

    n=n+1

  endif

enddo count_vaild_data

nt=n

rewind(11)

allocate(x(nt))

n=0

skip_comment: do

  read(11,'(A)',iostat=ios)strm

  if(ios<0)exit

  if(strm(1:1) == "#")then

    cycle

  else

    n=n+1

    read(strm,*)x(n)

  endif

enddo skip_comment

print *

do n=1,nt

  print *,'x(', n ,')=',x(n)

end do

print *

exe=test_read

infle="test_data.txt"

cat << EOF > $infle

# Comment 1

# Comennt 2

1

2

3

# Comment 4

4

5

EOF

echo "${infle}: "

cat $infle

echo

namelist=${exe}.namelist.txt

cat <<EOF > $namelist

&para

infle=$infle

&end

EOF

$exe < $namelist

exit 0