请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6603|回复: 13

[分享资料] 求助 关于时间相关系数的空间分布

[复制链接]

新浪微博达人勋

发表于 2015-6-15 10:18:17 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
最近在学习相关系数,根据论坛上牛人们给出的程序,自己也试了下,可是用fortran算的相关系数和直接用grads里tcorr做出来的图竟不一样,看了半天也没看明白到底是怎么回事,求牛人指导。
fortran:
parameter(nx=144,ny=73,nz=1,n=52,undef=-9.96921e+36)
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),t(nx,ny,nz)
open(1,file='D:\1957_2014\ecof1.txt')
read(1,*) (ty(it),it=1,n)
close(1)

open(2,file='D:\1957_2014\slp.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 j=1,ny
  do i=1,nx
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/n
wa(i,j,iz)=wa(i,j,iz)/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
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 j=1,ny
do i=1,nx
if(tw(i,j,iz,it)==undef)then
r(i,j,iz)=0.0
endif
enddo;enddo;enddo;enddo
open (3,file="D:\1957_2014\rslp.grd",form='binary')
do iz=1,nz
  do j=1,ny
do i=1,nx
write (3) r(i,j,iz)
enddo;enddo;enddo;
end


grads:
'reinit'
'open D:\1957_2014\ya.ctl'
'open D:\1957_2014\slp.ctl'
'set dfile 1'
'set x 1'
'set y 1'
'set t 1 52'
'define a=ya'
'set dfile 2'
'set x 1 144'
'set y 1 73'
'set t 1'
'define r=tcorr(ya,slp.2,t=1,t=52)'
'set gxout shaded'
'set black -0.2732 0.2732'
'd r'
'set gxout contour'
'd r'

'printim   D:\1957_2014\slp.png white'
'c'
'reinit'

fortran算的相关系数

fortran算的相关系数

用grads里tcorr算的相关系数

用grads里tcorr算的相关系数
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-6-15 11:53:29 | 显示全部楼层
我用过tcorr函数,这个函数是没有问题的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-6-15 12:11:25 | 显示全部楼层
刷牙 发表于 2015-6-15 11:53
我用过tcorr函数,这个函数是没有问题的

那就是fortran的问题了,我照着论坛里的帖子弄得,能帮我看看什么地方不对吗?我看了好几遍也没看出来。。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-6-15 12:34:01 | 显示全部楼层
zhaozhao 发表于 2015-6-15 12:11
那就是fortran的问题了,我照着论坛里的帖子弄得,能帮我看看什么地方不对吗?我看了好几遍也没看出来。 ...

Fortran不熟,我用matlab和tcorr结果对比过,虽然有一些差异但是问题不大,你就直接用tcorr就行了呗
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-6-16 08:19:02 | 显示全部楼层
采纳大家的建议,应该没错!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-6-28 08:31:19 | 显示全部楼层
wsjqxjy 发表于 2015-6-16 08:19
采纳大家的建议,应该没错!!

谢谢,已经解决了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-11-10 08:20:32 | 显示全部楼层
刷牙 发表于 2015-6-15 12:34
Fortran不熟,我用matlab和tcorr结果对比过,虽然有一些差异但是问题不大,你就直接用tcorr就行了呗

matlab程序能让我参考一下吗
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-11-10 08:20:54 | 显示全部楼层
{:5_231:}
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2015-12-23 09:22:01 | 显示全部楼层
谢谢各位的指导
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-5 14:20:21 | 显示全部楼层
请问我设置set black -0.3247 0.3247,但是画出来的图是-4和4之间不填色,这是为什么呢
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表