相关系数的空间分布

wpd0518

'reinit'
'sdfopen f:\fore.nc'
'sdfopen f:\obs.nc'
'set gxout fwrite'
'set fwrite e:\corr.grd'
'set x 1 134'
'set y 1 61'
'set t 1'
i=1
while(i<=12)
'set z 'i''
'define avefore=ave(sst,t=1,t=1560)'
'define aveobs=ave(sst.2,t=1,t=1560)'
'define a=ave((avefore-sst)*(aveobs-sst.2),t=1,t=1560)'
'define b=sqrt(ave(pow(avefore-sst,2),t=1,t=1560))'
'define c=sqrt(ave(pow(aveobs-sst.2,2),t=1,t=1560))'
'define corr=a/(b*c)'
'd corr'
i=i+1
endwhile
'disable fwrite'
'reinit'

mnkingking

应该没有问题，你检查哈原始场数据看看~