程序出不来、、求大神指点、、谢谢

appy!!!

这个是求30年850hPa,500hpa,200hpa三层的气候态，资料是三层30年的候资料
parameter(nx=144,ny=73,nz=3)
        parameter(ns=nx*ny*nz,np=73,mt=30)
            character*4 ch(mt)

        data ch/'1981','1982','1983','1984','1985','1986','1987','1988','1989','1990','1991','1992','1993','1994','1995','1996','1997','1998','1999','2000','2001','2002','2003','2004','2005','2006','2007','2008','2009','2010'/

        real th(ns,np),aveth(ns,np)
         

        open(30,file='f:/heigh/hh/avehgt.grd',form='binary')


        do 2013 kk=1,mt          !1981-2010

        open(22,file='f:/heigh/hh/'//ch(kk)//'.grd',form='binary')

        write(*,*)ch(kk)

        read(22) th
        close(22)

        do i=1,ns
        do k=1,np
        aveth(i,k)=aveth(i,k)+th(i,k)/float(mt)
        enddo
        enddo

2013        continue

        write(30)aveth

        end                求助啊、、

xo2001

提示这个grd文件已经读完了，你还要继续读取，就出错了