fortran程序应该是在批量输出txt文件时出错，出不来结果，求助！

不断长大

求助各位大神帮我看看我这个程序。。。是想求站点观测资料的每一年的月平均值，一共有50年，因为有700多个测站，所以需要进行批处理，以及批量的输出文件。下面的程序是给出三个测站资料作为测试。。。此程序没有报错，但是出不来结果。。估计是批量输出的编写有问题，希望给位大神帮忙看下。。非常感谢 50434.txt (1.1 MB, 下载次数: 14)

2013-11-29 16:16 上传

点击文件名下载附件


50353.txt (1.1 MB, 下载次数: 19)

2013-11-29 16:16 上传

点击文件名下载附件


50136.txt (1.1 MB, 下载次数: 10)

2013-11-29 16:16 上传

点击文件名下载附件


ave.f90 (3.05 KB, 下载次数: 7)

2013-11-29 16:16 上传

点击文件名下载附件


program aveimplicit none
   character(79) aaaa
   character(79) buffer

   integer station   ! 站点数
   integer parameter :: n=50   ! 资料序列总年数
   interger  y, m, d, p, year, runnian, yy(31,12,n), mm(31,12,n), dd(31,12,n)   !年、月、日等变量
   real tave(31,12,n), tmax(31,12,n), tmin(31,12,n),  tavemonth(12,n), tmaxmonth(12,n), tminmonth(12,n), tavemonth_ave(12,n), tmaxmonth_ave(12,n), tminmonth_ave(12,n)   

   open(10, file="dir.txt", access="sequential", status="old")
   do station =1, 786
      read(10, format="(A79)") aaaa
!     write(*,*) aaaa
!  end  do
      open(11, file=aaaa, access="sequential", status="old")
!      read(11. format="(A79,24X)")buffer               ! 批量处理txt文件
       open(12, file=aaaa_ave, access="sequential", status="new")  !计算结果保存的目标文件
!  write(*,*)buffer


      do y=1, n                                             ! 年循环
           year=1962+y    !year是时间起始年为1963年
           call sub(year, runnian)  !判断是否闰年，输入year，输出闰年值
        do m=1, 12                                       ！月循环
          select case(m)
             case(1,3,5,7,8,10,12)
               p=31

             case(2)
               if (runnian==0) then
                         p=28
               else
                         p=29
                       end if
             case(4,6,9,11)
               p=30
          end select

          tavemonth(m,y)=0
          tmaxmonth(m,y)=0
          tminmonth(m,y)=0
           do d=1,p                                        !日循环
             read(11, format="(A79,24X)") yy(d,m,y), mm(d,m,y), dd(d,m,y), tave(d,m,y),tmax(d,m,y), tmin(d,m,y)
             tavemonth(m,y)=tavemonth(m,y)+tave(d,m,y)
             tmaxmonth(m,y)=tmaxmonth(m,y)+tmax(d,m,y)
             tminmonth(m,y)=tminmonth(m,y)+tmin(d,m,y)
           end do                                            !日循环结束


        tavemonth_ave(m,y)=tavemonth(m,y)/real(p)
        tmaxmonth_ave(m,y)=tmaxmonth(m,y)/real(p)
        tminmonth_ave(m,y)=tminmonth(m,y)/real(p)      
        write(12,*)   tavemonth_ave(m,y), tmaxmonth_ave(m,y),tminmonth_ave(m,y)
        end do                                              !月循环结束

      end do                                                 !年循环结束



    end do                                                   !批处理结束

    close(12)
    close(11)
    close(10)

    subroutine sub(year, runnian)            !判断闰年的小程序
      implicit none
        integer year, runnian
        runnian=0
        if(mod(year, 100)==0)  then
                if(mod(year, 400)==0) then
                        runnian=1
                end if
        else
                if(mod(year, 4)==0) then
                        runnian=1
                end if
        end if
        if (runnian==1) then
                write(*,*) year,"年是闰年"
        else
                write(*,*) year,"年不是闰年"
        end if
     return
    end subroutine

andrewsoong

程序没看出什么问题，看来我功力不够~~~

不断长大

已解决，就是输出的问题

yeah...

不断长大 发表于 2013-11-30 17:18

                               
登录/注册后可看大图

已解决，就是输出的问题

输出的问题？请问一下你的数据是怎么读的呀？

不断长大

yeah... 发表于 2013-11-30 22:11

                               
登录/注册后可看大图

输出的问题？请问一下你的数据是怎么读的呀？

就是按行读，每一个数据都放在数组中

yeah...

不断长大 发表于 2013-12-1 09:40

                               
登录/注册后可看大图

就是按行读，每一个数据都放在数组中

我做的和你的差不多 不过我只需要做一个站就可以，和你的程序编出来差不多，可是我发现用excel验证时出现了问题哎 请问一下你这个程序编出来的和用excel验证的相同么

不断长大

yeah... 发表于 2013-12-2 21:30

                               
登录/注册后可看大图

我做的和你的差不多 不过我只需要做一个站就可以，和你的程序编出来差不多，可是我发现用excel验证时出现 ...

也不同。。。但是找到原因啦，因为数据问题，数据有些缺测