相关系数问题

MARZOLINI

求中国某区域的降水与海温的相关系数。
我是想把中国某区域的每个格点都与海温场做相关，例如中国某区域的格点数是lat：25*lon：37，然后把每个格点都与海温场做时间相关，得到的时间相关系数相加然后除以格点数25*37（类似于我得到25*37个海温与降水的时间相关系数场然后求平均）。
下面是我的脚本：
我在将这些相关系数累加的过程中，发现此循环没有累加成功，输出的ccr_1995的值都是缺设值，请求各位大神帮我看看这是什么问题，或者是我求相关系数的思路就不对。
nerravg_1995和nsst_1995的第三维是时间维。
ccr_1995(:,:)=0.0
    do i=0,24,1
        do j=0,36,1
            ccr_1995(:,:)=ccr_1995(:,:)+escorc(nerravg_1995(i,j,:), nsst_1995(:,:,:))
        end do
    end do

   print(ccr_1995(:,:))

APTX4869

nsst_1995(:,:,:))这是三位的而nerravg_1995(i,j,:)是一维，你这个循环应该是需要都改成一维 的

MARZOLINI

应该不需要改为一维得吧，因为我如果指定i=0，j=0，不是循环，
则escorc(nerravg_1995(0,0,:), nsst_1995(:,:,:))是能计算出结果的。

MARZOLINI

APTX4869 发表于 2018-9-11 17:24
nsst_1995(:,:,:))这是三位的而nerravg_1995(i,j,:)是一维，你这个循环应该是需要都改成一维 的

应该不需要改为一维得吧，因为我如果指定i=0，j=0，不是循环，
则escorc(nerravg_1995(0,0,:), nsst_1995(:,:,:))是能计算出结果的。

anran

请问楼主解决了吗，我最近也遇到了海温相关结果输出全是缺测值

好久不见

循环里不是应该写成 ccr_1995(i,j)=escorc(nerravg_1995(i,j,:), nsst_1995(:,:,:)) 吗。。。。。。