8#的程序有基础问题,可能没有考虑到5#补充的内容,注意到x(j,k)的由来和公式中x2()的修改,程序就好办多了。我尝试写了一条程序:
parameter(m=,n=,o=)!注意此处mno的定义和文章少有区别,m是年,n是月,o是个点数目
real x(m,n,o),x1(n,o),x2(m,o),xk(o)
do j=1,n
do k=1,o
do i=1,m
s=s+x(i,j,k)
enddo
x1(j,k)=s/m
enddo
enddo
do k=1,o
do j=1,n
s=s+x1(j,k)
enddo
xk(k)=s/n
enddo
do j=1,n
do k=1,o
x2(j,k)=x1(j,k)-xk(k)
enddo
enddo