grads 合成分析距平场t检验

18200798924

用下面这个gs文件做距平高度场的检验时发现没有过检验的区域，把临界值2.447减小到1.447有过检验的区域，但是距平中心都不过检验，老师说不正常，应该距平中心更容易过检验才对，但是我检查不出错误，希望大家帮我看看。
'reinit'
'sdfopen F:\LJL-Graduation-thesis\Data\Reanalysis\height\hgt.mon.mean.nc'
'sdfopen  F:\LJL-Graduation-thesis\Data\Reanalysis\wind\uwnd.mon.mean.nc'
'sdfopen  F:\LJL-Graduation-thesis\Data\Reanalysis\wind\vwnd.mon.mean.nc'
'set grid off'
'set grads off'
'set lon  0 180'
'set lat  0 80'
'set lev 500'

'set xlint 20'
'set ylint 20'
'set map 1 1 6'
'set mpdset cnworld'


'define meanhgt10=ave(hgt.1,time=oct1982,time=oct2019,12)'*气候平均


'define parthgt10=(hgt.1(time=oct1982)+hgt.1(time=oct1985)+hgt.1(time=oct1987)+hgt.1(time=oct1988)+hgt.1(time=oct1997)+hgt.1(time=oct1998)+hgt.1(time=oct1999))/7'*异常大年平均

'define diffhgt10=parthgt10-meanhgt10'*异常大年相对于气候平均距平

'diffhgt82= hgt.1(time=oct1982)-parthgt10'*用于下面计算标准差
'diffhgt85= hgt.1(time=oct1985)-parthgt10'
'diffhgt87= hgt.1(time=oct1987)-parthgt10'
'diffhgt88= hgt.1(time=oct1988)-parthgt10'
'diffhgt97= hgt.1(time=oct1997)-parthgt10'
'diffhgt98= hgt.1(time=oct1998)-parthgt10'
'diffhgt99= hgt.1(time=oct1999)-parthgt10'

's=pow((pow(diffhgt82,2)+pow(diffhgt85,2)+pow(diffhgt87,2)+pow(diffhgt88,2)+pow(diffhgt97,2)+pow(diffhgt98,2)+pow(diffhgt99,2))/7,1/2)'*标准差
'a=s/pow(7,1/2)'
't005=diffhgt10/a'*t值

*风场
'define meanuwnd10=ave(uwnd.2,time=oct1982,time=oct2019,12)'
'define partuwnd10=(uwnd.2(time=oct1982)+uwnd.2(time=oct1985)+uwnd.2(time=oct1987)+uwnd.2(time=oct1988)+uwnd.2(time=oct1997)+uwnd.2(time=oct1998)+uwnd.2(time=oct1999))/7'
'diffuwnd10=partuwnd10-meanuwnd10'

'define meanvwnd10=ave(vwnd.3,time=oct1982,time=oct2019,12)'
'define partvwnd10=(vwnd.3(time=oct1982)+vwnd.3(time=oct1985)+vwnd.3(time=oct1987)+vwnd.3(time=oct1988)+vwnd.3(time=oct1997)+vwnd.3(time=oct1998)+vwnd.3(time=oct1999))/7'
'diffvwnd10=partvwnd10-meanvwnd10'

define_colors    *检验，过检验区域填色灰色
'set gxout shaded'
'set clevs -2.447 0 2.447'
'set ccols 51 0 0 51'
'd t005'
'cbarn 1 0'

'set gxout vector'
'set arrscl 0.5 3'
'set arrowhead 0.08'
'set arrlab on'
'd skip(diffuwnd10,3,3);skip(diffvwnd10,3,3)'

'set gxout contour'
'set clevs -30 -25 -20 -15 -10 -5 0  5 10 15 20 25 30'
'set cthick 5'
'set ccolor 1'
'd diffhgt10'

'set line 2 1 15'
'draw rec 4.75 3.35 5.55 4.0'

'draw title 500hpa difference of zone4.1 Arctic icec high index in October'

'printim F:\LJL-Graduation-thesis\Data\advance-prove\height-wind\500hpa\zone4.1-high-index-500hpa-uv10-hgt10.png white'

river

我觉得是你的气候平均计算有误，
'define meanhgt10=ave(hgt.1,time=oct1982,time=oct2019,12)'*气候平均
红色部分你要弄懂它代表什么意思才行。

你的要计算的是气候平均，我的理解是那么多年的数据最终得到一个平均数，而不是每年一个平均数吧？