【已解决】回归系数的程序

muggle

我打开了一个原始场，一个平均场的grd文件，分别放在一个四维数组和一个三维数组里，现在要计算每个格点31年的回归系数，并且存放在一个三位数组里面，调适后程序没有错，但是运行结果却显示这样，说不存在unit11？这是为什么，路径明明是对的。以下是错误截图和我的程序，求指点



PROGRAM main
implicit none
integer,parameter::n=31

real yave

integer i,m,sum
PARAMETER ix=29,iy=17,iz=17,it=62
real x,y,z,t
integer uv(ix,iy,iz,it),year(it),pjuv(ix,iy,iz),k(ix,iy,iz)


m=1979
sum=0
do i=1,n
year(i)=m
m=m+1
sum=sum+year(i)
end do
yave=sum/n



open (11,file='f:/yq/zl/uvwinter-nc.grd',form='binary')
do t=1,it
  do z=1,iz
    do y=1,iy
          do x=1,ix
        read(11) uv(x,y,z,t)
      end do
        end do
  end do
end do
close(11)

open (22,file='f:/yq/pj/pjuvwinter-nc.grd',form='binary')
  do z=1,iz
    do y=1,iy
          do x=1,ix
        read(11) pjuv(x,y,z)
      end do
        end do
  end do
close(22)


  do z=1,iz
    do y=1,iy
          do x=1,ix
            do t=32,it
k(x,y,z)=(uv(x,y,z,t)*yave-n*pjuv(x,y,z)*yave)/(uv(x,y,z,t)*uv(x,y,z,t)-n*pjuv(x,y,z)*pjuv(x,y,z))
            end do
          end do
         end do   
  end do


open (33,file='f:/yq/hgxs/hgxs-nc.grd',form='binary')
do z=1,iz
    do y=1,iy
          do x=1,ix
     read(33) k(x,y,z)
      end do
        end do
  end do
close(33)

end program




[img]file:///C:\Users\Administrator\AppData\Roaming\Tencent\Users\369956658\QQ\WinTemp\RichOle\RZ1$DR8[KI@`(PRX9~7]2{E.png[/img]

Fred2016

多回复几个就可以到5了，不是吗？

muggle

我还以为哪个大神出现了

言深深

写错东西了，22里面的read你写成了11

muggle

言深深 发表于 2016-4-20 12:29
写错东西了，22里面的read你写成了11

oh...no
摸索了半天原来是这里出了错，出现这种低级失误真的不应该
调整之后仍然出现了forrtl severe（24）我去摸索一下，如果还是不明白的，希望前辈你多多指教

梦轩蔷薇

不懂啊，路过的初学者一枚！

卢姁luxu

33是不是应该为write

muggle

卢姁luxu 发表于 2017-2-11 10:40
33是不是应该为write

是的！就是这个问题