请问不同长度文件用一维数组怎么顺序读取？

whmooncool

我这样大概100个站点，每个站点开始的年月日不同，数据长度不同，比如有的站点是从51年开始记录，有的从60年记录。我想顺序读下去这些个文件，就编写了如下程序，程序链接时没有错误，但每次读取的时候总是显示数组溢出，可能是红色的问题。请问我的程序改怎么改？
module typedef
  type rheight
integer :: station, lat,lon,height,year!  站点号，纬度，经度，高度
  end type
type day   
integer :: year,day
  end type
end module
program main
  use typedef
  implicit none
  integer, parameter :: fileid=9
integer, parameter :: days=16802  ！选取最大文件的长度
character(len=80)  :: tempstr
  type(rheight) :: r(days) !原数据
  type(day) :: d(days)       ! 存储提取出数据
  integer i,j,k,n
  integer a,b,c

  character stationfile(7)*9
  data stationfile/'52378.txt','52495.txt','52674.txt','52679.txt','52681.txt','52787.txt','52797.txt'/
   n=1
  do  j=1,58

  open(n, file='f:\50-07/'//trim(stationfile(j))//'')
  


  ! 用循环读入每日站点数据
  do i=1,days1
read(n,*) r(i)%station, r(i)%lat,r(i)%lon,r(i)%height,r(i)%   
if ((r(i)%year)>1960) then
a=i-3652
d(a)%year=r(i)%year
d(a)%day=r(i)%day
end if                                         
  end do
  close(n)
   k=n+60
  ! 在文件中输出日降水量
  
  open(k, file='f:\d'//trim(stationfile(j))//'')
  
  do b=1,days2
  write(k,"(I4,I2,I2,f8.1)") d(b)%year, d(b)%month,d(b)%day,d(b)%dr
  end do
  close(k)
   n=n+1
  end do

  stop
end program

mofangbao

首先确定你组合成的文件名是没有错误的
然后像这种行数不一致的，建议read的时候加上 end=的标签，这样读取到文件末尾的时候就会自动跳转到指定行去执行
还有行数不定时的循环建议使用 continue的方式来循环，而不是do一个固定的数值，具体方法我的很多帖子里面都用过，你自己找找看

topmad

可以先判断一下文件的行数，然后用动态数组

whmooncool

topmad 发表于 2012-3-21 13:37

                               
登录/注册后可看大图

可以先判断一下文件的行数，然后用动态数组

可以详细说下嘛？谢谢

topmad

whmooncool 发表于 2012-3-22 11:12

                               
登录/注册后可看大图

可以详细说下嘛？谢谢

详细请见fortran版求取文件行数的帖子