關於算式不成立程序無法運行

li913

[Fortran] 纯文本查看 复制代码

program pwv_mon_avg
implicit none
integer :: stat, stat2, i, year, mon, mon_num(12)
character :: filename*12, filename_new*16
real(8) :: pwv_mm, pwv_mon_mean(12)

open (11,file='list.txt')
close (11,status='delete')
call system ('for %f in (*_mon.pwv) do echo %f >> list.txt')
open (11,file='list.txt')
do 
    read (11,*,iostat=stat) filename
    if (stat/=0) exit
    filename_new=filename(1:4)//'_mon_avg.pwv'
    open (33,file=filename_new)
    open (22,file=filename)
    mon_num(:) = 0
    pwv_mon_mean(:) = 0.0d0
    do 
      read (22,*,iostat=stat2) year, mon, pwv_mm
      if (stat2/=0) exit
      pwv_mon_mean(mon) = pwv_mon_mean(mon) + pwv_mm
      mon_num(mon) = mon_num(mon) + 1
    end do
    do i = 1, 12
        if( mon_num(i)/=0 ) pwv_mon_mean(i) = pwv_mon_mean(i) / mon_num(i)
        write (33,'(i2.2,3x,f12.7)') i, pwv_mon_mean(i)
    end do
    close(22)
    close(33)
end do
close (11,status='delete')
write(*,*) "completed..."
read(*,*)
end program