利用Fortran求相关系数

wangyijie

       大家好，这个求相关系数的程序我找了好久也没有发现到底是哪里出毛病了，恳求大佬们帮助！

  !求个例1982夏季RMM与冬季太平洋地区OLR相关系数（西经90-东经120，南北纬15）
    program  main
    parameter(nx=85,ny=13,nt=123)
    integer it,i,j
    real rmm(nt),olr(nx,ny,nt)
    real ya,wa(nx,ny),syw(nx,ny),sy,sw(nx,ny),r(nx,ny)
    open(1,file='d:\biye\xiaRMM.grd',form='binary')
    open(3,file='d:\biye\olr.grd',form='binary')
    !读取olr数据
    do it=1,nt
        do j=1,ny
            do i=1,nx
                read(3) olr(i,j,it)
            enddo
        enddo
    enddo
    close(3)
    do it =1,nt
        read(1) rmm(it)
    enddo
    close(1)
    !进行相关系数的计算
    do i=1,nx
        do j=1,ny
            ya=0.0
            wa(i,j)=0.0
            do it=1,nt
                ya=ya+rmm(it)
                wa(i,j)=wa(i,j)+olr(i,j,it)
            enddo
            ya=ya/real(nt)
            wa(i,j)=wa(i,j)/real(nt)
            syw(i,j)=0.0
            sy=0.0
            sw(i,j)=0.0
            do it=1,nt
                syw(i,j)=syw(i,j)+(rmm(it)-ya)*(olr(i,j,it)-wa(i,j))
                sy=sy+(rmm(it)-ya)**2
                sw(i,j)=sw(i,j)+(olr(i,j,it)-wa(i,j))**2
            enddo
            r(i,j)=0.0
            r(i,j)=syw(i,j)/sqrt(sy*sw(i,j))
        enddo
    enddo
    open(4,file='d:\biye\r.grd',form='binary')
    do j=1,ny
        do x=1,nx
            read(4) r(i,j)
        enddo
    enddo
    end

wangyijie

问题已经解决了，谢谢大家！