一个求年代间差异的fortran程序让我变成死循环了 求指教

凉茶

总共54年 要求年代间差异 即后27年与前27年的差 下面是我的程序 求纠错 谢谢
program ndj
implicit none
integer,parameter::nt=54,nx=144,ny=73,nz=8
integer:: it,ix,iy,iz
real wwinave(nx,ny,nz) !存放各高度各网点的年平均值，这里只是3维数组
real winave(nx,ny,nz,nt) !实际测量数据
real xx(nx,ny,nz),yy(nx,ny,nz),bb(nx,ny,nz)  !分别存放前27年，后27年平均湿度和年代间差
open(11,file='d:\lunwen\yw\win1.grd',form='binary')
do it=1,nt
  do iz=1,nz
           read(11) ((winave(ix,iy,iz,it),ix=1,nx),iy=1,ny)
enddo
enddo
wwinave=0
do it=1,nt/2
do iz=1,nz
    do iy=1,ny
      do ix=1,nx
            wwinave(ix,iy,iz)=wwinave(ix,iy,iz)+winave(ix,iy,iz,it)!将各高度网格点按年份累加
                !yy=yy+it
         print*,winave(ix,iy,iz,it)
          end do
    end do
  end do
end do
do iz=1,nz
    do iy=1,ny
      do ix=1,nx
          xx(nx,ny,nz)=wwinave(ix,iy,iz)/27
          !print*,wwinave(ix,iy,iz)
           end do
    end do
end do
wwinave=0
do it=28,nt
do iz=1,nz
    do iy=1,ny
      do ix=1,nx
            wwinave(ix,iy,iz)=wwinave(ix,iy,iz)+winave(ix,iy,iz,it)!将各高度网格点按年份累加
                !yy=yy+it
          end do
    end do
  end do
end do
do iz=1,nz
    do iy=1,ny
      do ix=1,nx
          yy(nx,ny,nz)=wwinave(ix,iy,iz)/27
           end do
    end do
end do
do iz=1,nz
    do iy=1,ny
      do ix=1,nx
          bb(nx,ny,nz)=yy(nx,ny,nz)-xx(nx,ny,nz)
           end do
    end do
end do
open(12,file='d:\lunwen\yw2\ndj.grd',form='binary')
do iz=1,nz
       write (12) ((bb(ix,iy,iz),ix=1,nx),iy=1,ny)
enddo
close(12)
open(13,file='d:\lunwen\yw2\ndj.txt')
do iz=1,nz
       write (13,*) ((bb(ix,iy,iz),ix=1,nx),iy=1,ny)
enddo
close(13)
end

凉茶

出的文件数值都是0  怎么回事 啊 求大神指导

mofangbao

凉茶 发表于 2012-5-9 18:36

                               
登录/注册后可看大图

出的文件数值都是0  怎么回事 啊 求大神指导

太长了，自己先定位到出错的那一段代码再说吧~