grads t检验值的输出

天玖Y · 发表于 2014-12-30 16:19:31

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

强弱季风年差值场的t检验，想把T检验的值输出成.grd格式，但是作图发现有点问题，与原图中检验的值对应不起来~。不晓得哪里出错了，求指教呀~~

{:soso_e101:}
---GS
'reinit'
'open I:\OBSYYW\monsoon\air\obs-air.ctl'
'set gxout fwrite '
'set fwrite I:\OBSYYW\monsoon\air\t3.grd'
'set lon 0 360'
'set lat -90 90'
'set z 1'
*********计算强年的平均值，方差****************************
'define save=(air(t=5)+air(t=7)+air(t=10)+air(t=17))/4.0'
'define spp=(pow(air(t=5)-save,2)+pow(air(t=7)-save,2)+pow(air(t=10)-save,2)+pow(air(t=17)-save,2))/4.0'
*********计算ruo年的平均值，方差****************************
'define wave=(air(t=4)+air(t=13)+air(t=19))/3.0'
'define wpp=(pow(air(t=4)-wave,2)+pow(air(t=13)-wave,2)+pow(air(t=19)-wave,2))/3.0'***********计算EP与CP两者MJO的差值**************************
'define dif=save-wave'
********计算检验值test***************************************
'define test=(save-wave)/(sqrt(3*spp+2*wpp)/5)*sqrt(1/4+1/3))'
'define tt=abs(test)-2.57'
'd tt'
'disable fwrite'
;file:///C:\Users\hp\AppData\Roaming\Tencent\Users\364768116\QQ\WinTemp\RichOle\FKO$@WR{WNXT)TV1(`EJJ7D.png
QQ图片20141230161233.jpg

file:///C:\Users\hp\AppData\Roaming\Tencent\Users\364768116\QQ\WinTemp\RichOle\FKO$@WR{WNXT)TV1(`EJJ7D.png

来自群组: 南气院校友总会

zhouling · 发表于 2014-12-30 20:32:57

楼主可以参考这篇文献

youko · 发表于 2014-12-30 21:30:43

这位小哥，你试一下把set lon 0 360;set lat -90 90;改成set x 1 144;set y 1 73试试~~~~~

天玖Y · 发表于 2014-12-30 21:32:01

zhouling 发表于 2014-12-30 20:32
楼主可以参考这篇文献

好的谢谢学习一下

天玖Y · 发表于 2014-12-30 21:39:29

youko 发表于 2014-12-30 21:30
这位小哥，你试一下把set lon 0 360;set lat -90 90;改成set x 1 144;set y 1 73试试~~~~~

感谢这位小哥这回终于是对了

youko · 发表于 2014-12-30 21:52:49

大家互相学习，共同进步~

逍遥仙境 · 发表于 2015-3-11 15:41:47

明白了，谢谢！

小小的l · 发表于 2015-5-9 13:34:18

{:lxm_24:}

我是彤-Tina · 发表于 2015-11-21 21:28:39

搜到这个帖子正好解决了我的问题~

undisclosed · 发表于 2017-5-7 17:24:14

youko 发表于 2014-12-30 21:30
这位小哥，你试一下把set lon 0 360;set lat -90 90;改成set x 1 144;set y 1 73试试~~~~~

这位小哥，你试一下

		自动登录	找回密码
密码			立即注册

[脚本编辑] grads t检验 值的输出

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

[脚本编辑] grads t检验值的输出