求助，海温与一维序列做相关后，陆地上出现数值

Crystal_PP

如题，我做了相关后，陆地上就出现数值，请问是怎么回事？

mofangbao

怎么做的，能说下过程吗？

天目神眉

支持一下庞Y，但是我不懂

冰枫

这个可能是GRADS边界和海温等资料的边界差异造成的，范围不大应该不要紧的，你用basemap覆盖就好了

Crystal_PP

天目神眉 发表于 2012-5-8 22:19

                               
登录/注册后可看大图

支持一下庞Y，但是我不懂

咦，认识？

Crystal_PP

mofangbao 发表于 2012-5-8 22:19

                               
登录/注册后可看大图

怎么做的，能说下过程吗？

就是简单的先提取多年的海温资料，再把它与一个一维的数值做了个相关，求相关系数的分布。提取的海温资料没有出现陆地上有值的情况，但是一做相关陆地上海岸线附近就出现数值，要怎么修改啊？

mofangbao

Crystal_PP 发表于 2012-5-9 08:44

                               
登录/注册后可看大图

就是简单的先提取多年的海温资料，再把它与一个一维的数值做了个相关，求相关系数的分布。提取的海温资料 ...

你把海岸附近的数值输出看看

Crystal_PP

mofangbao 发表于 2012-5-9 09:20

                               
登录/注册后可看大图

你把海岸附近的数值输出看看

我用 fwrite 命令截取海温nc资料
'set x 1 180'
'set y 1 89'
...

得到的海温，仔细一看海岸线附近的值还是没有，是我截取的海温资料有问题吗？

                               
登录/注册后可看大图

mofangbao

Crystal_PP 发表于 2012-5-9 09:37

                               
登录/注册后可看大图

我用 fwrite 命令截取海温nc资料
'set x 1 180'
'set y 1 89'

能把你求相关的脚本贴出来吗，我这里也有海温资料，可以试试~

Crystal_PP

mofangbao 发表于 2012-5-9 09:39

                               
登录/注册后可看大图

能把你求相关的脚本贴出来吗，我这里也有海温资料，可以试试~

parameter(nx=181,ny=89,nz=1,n=52,undef=32767)
dimension ty(n),tw(nx,ny,nz,n)
real ya,wa(nx,ny,nz),syw(nx,ny,nz),sy,sw(nx,ny,nz),r(nx,ny,nz)


open(1,file='i.txt')
read(1,*) (ty(it),it=1,n)
close(1)

open(2,file='sst.6.grd',form='binary')
read(2) ((((tw(i,j,iz,it),i=1,nx),j=1,ny),iz=1,nz),it=1,n)
close(2)

do iz=1,nz
do i=1,nx
do j=1,ny
ya=0.0
wa(i,j,iz)=0.0
do it=1,n
ya=ya+ty(it)
wa(i,j,iz)=wa(i,j,iz)+tw(i,j,iz,it)
end do
ya=ya/real(n)
wa(i,j,iz)=wa(i,j,iz)/real(n)

syw(i,j,iz)=0.0
sy=0.0
sw(i,j,iz)=0.0
do it=1,n
syw(i,j,iz)=syw(i,j,iz)+(ty(it)-ya)*(tw(i,j,iz,it)-wa(i,j,iz))
sy=sy+(ty(it)-ya)*(ty(it)-ya)
sw(i,j,iz)=sw(i,j,iz)+(tw(i,j,iz,it)-wa(i,j,iz))*(tw(i,j,iz,it)-wa(i,j,iz))
end do
syw(i,j,iz)=syw(i,j,iz)/real(n)
sy=sy/real(n)
sw(i,j,iz)=sw(i,j,iz)/real(n)
r(i,j,iz)=0.0
r(i,j,iz)=syw(i,j,iz)/sqrt(sy*sw(i,j,iz))
enddo;enddo;enddo

do it=1,n
do iz=1,nz
do i=1,nx
do j=1,ny
if(tw(i,j,iz,it)==undef)then
r(i,j,iz)=0.0
endif
enddo;enddo;enddo;enddo