fortran循环读入ERA月资料

cycycycycy

很奇怪，用了之前的方法，读不出来
大神帮忙看看吧！谢谢！
program EX1
integer,parameter::in=240,jn=121,zn=27,mon=12,tt=38
integer it,imo,iz,k,i,j
character timenames(tt)*4
real w(in,jn,zn,mon,tt),hgt(in,jn,zn,mon,tt),air(in,jn,zn,mon,tt),q(in,jn,zn,mon,tt),rh(in,jn,zn,mon,tt)
real u(in,jn,zn,mon,tt),v(in,jn,zn,mon,tt)
DATA timenames/'1979','1980','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','2011','2012','2013','2014','2015','2016'/

!!!需要100HPA以下的data
u=0.0
v=0.0
w=0.0
hgt=0.0
air=0.0
q=0.0
rh=0.0
do it=1979,2016
open(9,file='I:\data\ChenYLdata\'//trim(timenames(it-1978))//'.grd',form='binary',status='old')
  do imo=1,mon
    read(9)(((u(i,j,iz,imo,it-1978),i=1,in),j=1,jn),iz=1,zn)
    read(9)(((v(i,j,iz,imo,it-1978),i=1,in),j=1,jn),iz=1,zn)
    read(9)(((w(i,j,iz,imo,it-1978),i=1,in),j=1,jn),iz=1,zn)
    read(9)(((hgt(i,j,iz,imo,it-1978),i=1,in),j=1,jn),iz=1,zn)
read(9)(((air(i,j,iz,imo,it-1978),i=1,in),j=1,jn),iz=1,zn)
    read(9)(((q(i,j,iz,imo,it-1978),i=1,in),j=1,jn),iz=1,zn)
    read(9)(((rh(i,j,iz,imo,it-1978),i=1,in),j=1,jn),iz=1,zn)
  enddo
close(9)
enddo
print*,u(12,12,1,1,1)


end