fortran编程

凌晨24点

integer::i,ii,j,jj,k,kk,r,aa
integer(kind=4),parameter::n=1321,mm=7320
integer,parameter::m=100,nn=1000
integer::id(n),year(n),month(n),day(n),clock(n),cp(n),yearyear(mm),monthmonth(mm),dayday(mm),llll
real::lat(n),long(n),ws(n),fs(n),latitude(m),longitude(m),pre(mm),rain(nn),wind(nn),winds(mm)
character::station(m)*5,stn(nn)*5,ystn(nn)*5,tmp1*1,tmp11*1,tmp2*2,tmp22*2,yr*4,sta(mm)*5,tmp4*6
real::yyr(nn),mmon(nn),dday(nn),da(nn),ye(nn),mo(nn)
open(10,file='E:\munich\lujing.txt')
read(10,*)
do i=1,n
read(10,*)id(i),year(i),month(i),day(i),clock(i),lat(i),long(i)
enddo
close(10)
open(11,file='E:\munich\stations.txt')
do j=1,m
read(11,*) station(j),latitude(j),longitude(j)
enddo
close(11)

open(13,file='E:\munich\pre.txt')
do r=1,mm
read(13,*)sta(r),yearyear(r),monthmonth(r),dayday(r),pre(r),winds(r)
enddo
close(13)

open(501,file='D:\tt\tt.txt')
do i=1,n
k=1
do j=1,m
if(sqrt((((latitude(j)-lat(i))*111)**2)+(((longitude(j)-long(i))*101)**2))<=200)then
stn(k)=station(j)
k=k+1
endif
enddo
!降水对应站点
do jj=1,k-1
rain(jj)=0.0
wind(aa)=0.0
do r=1,mm
if(yearyear(r)==year(i).and.monthmonth(r)==month(i).and.dayday(r)==day(i).and.sta(r)==stn(jj))then
rain(jj)=pre(r)
wind(aa)=winds(r)
endif
enddo
enddo
! write(day(i),'(i3,i4)')k-1,clock(i)
do jj=1,k-1
write(501,'(i4,i5,i3,i3,i4,i3,a6,f8.1,f8.1)')id(i),year(i),month(i),day(i),clock(i),k-1,stn(jj),rain(jj),wind(aa)
!write(501,'(i4,i5,i3,i3,i3,i4)')event(i),year(i),month(i),day(i),k-1,clock(i)
!do jj=1,k-1
! write(day(i),"(a5,f8.1)")event(i), stn(jj),rain(jj)
!write(501,"(i4,a6,f8.1)")event(i),stn(jj),rain(jj)
enddo
enddo
close(501)
end



这个程序运行到“open(11,file='E:\munich\stations.txt')"这句时提示array bounds exceeds,数组超出边界，但是我定义m=100是远远多于实际的行数的，数据的实际行数只有24行，我改成24行出错，改成25也出错，求高手指点！

mofangbao

程序有点多，楼主先定位一下然后再自己筛选一下试试啦

凌晨24点

mofangbao 发表于 2011-12-18 15:09

                               
登录/注册后可看大图

程序有点多，楼主先定位一下然后再自己筛选一下试试啦

谢谢啦，现在搞定了，多谢兄台提醒