对水汽通量图进行t检验，对显著性区域填色但填色图没有颜色

WOT

我用的是ERA数据，时间范围：1961-2020年7、8月，经纬度20S-60N，20E-120E，变量有z、q、u、v四个变量。
我现在的问题是可以画出水汽通量的u,v方向的分量，但是无法用填色图画出显著性区域，我用的是t检验的方法，不知道哪里出错了，请各位老师帮我看一下，谢谢！
*************1992-2020与1961-1985七、八月整层水汽通量差**************
'reinit'
'sdfopen E:\GrADS\testwww\ERA\shuiqitongliang.nc'




'set parea 1.5 10.0 1.5 7.0'
'set mproj scaled'
'set map 1 1 3'
'set mpdset cnworld cnriver'


'set xlopts -10'
'set ylopts -10'
'set xlint 10'
'set ylint 10'


*****************1992-2020年7、8月整层水汽通量***************
************七月*************************
'set z 1 16'
'define aqu7=u*q'
'define aqv7=v*q'
'set t 63'
'define aqud7=ave(aqu7,t+0,t=119,1)'
'define aqvd7=ave(aqv7,t+0,t=119,1)'
**************八月***********************
'set z 1 16'
'define aqu8=u*q'
'define aqv8=v*q'
'set t 64'
'define aqud8=ave(aqu8,t+0,t=120,1)'
'define aqvd8=ave(aqv8,t+0,t=120,1)'
************七、八月平均*********
'set z 1 16'
'define aqud=(aqud7+aqud8)/2'
'define aqvd=(aqvd7+aqvd8)/2'
******************水汽通量***************
'set z 16'
'define aqus=vint(z(z=16)/98,skip(aqud,7,7),500)'
'define aqvs=vint(z(z=16)/98,skip(aqvd,7,7),500)'
'define ashuiqi=mag(aqus,aqvs)'


*****************1961-1985年7、8月整层水汽通量***************
************七月*************************
'define bqu7=u*q'
'define bqv7=v*q'
'set t 1'
'define bqud7=ave(bqu7,t+0,t=49,1)'
'define bqvd7=ave(bqv7,t+0,t=49,1)'
**************八月***********************
'define bqu8=u*q'
'define bqv8=v*q'
'set t 2'
'define bqud8=ave(bqu8,t+0,t=50,1)'
'define bqvd8=ave(bqv8,t+0,t=50,1)'
************七、八月平均*********
'set z 1 16'
'define bqud=(bqud7+bqud8)/2'
'define bqvd=(bqvd7+bqvd8)/2'
******************水汽通量***************
'set z 16'
'define bqus=vint(z(z=16)/98,skip(bqud,7,7),500)'
'define bqvs=vint(z(z=16)/98,skip(bqvd,7,7),500)'
'define bshuiqi=mag(bqus,bqvs)'
************************差值***************************************
'define shuiqi=ashuiqi-bshuiqi'
'define qus=aqus-bqus'
'define qvs=aqvs-bqvs'
**************************水汽通量方差*********************
**********************1992-2020年水汽通量方差******
i=63
'set z 1 16 '
'define aaqu7=u*q'
'define aaqv7=v*q'
'define assum=aaqu7-aaqu7'
while (i<=119)
'set t ' i
'set z 1 16 '
'define aaqu7=u*q'
'define aaqv7=v*q'


'set t ' i+1
'set z 1 16 '
'define aaqu8=u*q'
'define aaqv8=v*q'


'set z 1 16 '
'define aaqud=(aaqu7+aaqu8)/2'
'define aaqvd=(aaqv7+aaqv8)/2'
'set z 16'
'define aaqus=vint(z(z=16)/98,skip(aaqud,7,7),500)'
'define aaqvs=vint(z(z=16)/98,skip(aaqvd,7,7),500)'


'define aashuiqi=mag(aaqus,aaqvs)'
'define assum=pow(aashuiqi-ashuiqi,2)+assum'
i=i+2
endwhile
'define avari=assum'


*******************1961-1985年水汽通量方差**********
i=1
'set z 1 16 '
'define bbqu7=u*q'
'define bbqv7=v*q'
'define bssum=bbqu7-bbqu7'
while (i<=49)
'set t ' i
'set z 1 16 '
'define bbqu7=u*q'
'define bbqv7=v*q'


'set t ' i+1
'set z 1 16 '
'define bbqu8=u*q'
'define bbqv8=v*q'


'set z 1 16 '
'define bbqud=(bbqu7+bbqu8)/2'
'define bbqvd=(bbqv7+bbqv8)/2'
'set z 16'
'define bbqus=vint(z(z=16)/98,skip(bbqud,7,7),500)'
'define bbqvs=vint(z(z=16)/98,skip(bbqvd,7,7),500)'
'define bbshuiqi=mag(bbqus,bbqvs)'


'define bssum=pow(bbshuiqi-bshuiqi,2)+bssum'


i=i+2
endwhile
'define bvari=bssum'
*************************总体方差**************
'define s=(avari+bvari)/(25+29-2)'


***************************t值求算**********************
'define stutest=shuiqi/sqrt(s/29+s/25)'


***********************************************************
'set lon 20 120'
'set lat -20 60'




'set gxout shaded'
'd abs(stutest)''run  cbarn.gs'




'set gxout vector'
'set arrlab on'
'set arrowhead -0.2'
'set arrscl 0.7 20'
'd qus;qvs'
'set grid off'






'run C:\OpenGrADS\Classic\scripts\axis.gs -type b -interval 20 -position o -suffix `3.'
'run C:\OpenGrADS\Classic\scripts\axis.gs -type L -interval 10 -position o -suffix `3.'


'draw title 1000-500hpa vapor flux difference'


'printim E:\GrADS\test\erashui.png  white x2500 y2000'
'disable print'

river

难道不应该画的是通过检验的水汽通量的值吗？直接画 t 的值是啥意思？

WOT

river 发表于 2023-3-17 12:42
难道不应该画的是通过检验的水汽通量的值吗？直接画 t 的值是啥意思？

是的，我也想看看t值的分布

river

WOT 发表于 2023-3-17 16:32
是的，我也想看看t值的分布

看看这个帖子http://bbs.06climate.com/forum.p ... p;extra=&page=1

WOT

river 发表于 2023-3-18 12:34
看看这个帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=31696&extra=&page=1

谢谢！但是为什么没颜色呢？