【求助】ERA数据画水汽通量填色图没颜色

WOT

这个脚本是画1992-2020年7、8月水汽通量的，想用填色图画一下水汽通量，但是运行脚本后没颜色，这是怎么回事呢？请各位老师帮忙看一下，大恩不言谢！
'reinit'
'sdfopen E:\GrADS\testwww\ERA\Geopotential.nc'
'sdfopen E:\GrADS\testwww\ERA\Specifichumidity.nc'
'sdfopen E:\GrADS\testwww\ERA\Uwnd.nc'
'sdfopen E:\GrADS\testwww\ERA\Vwnd.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'


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




*****************1992-2020年7、8月整层水汽通量***************
************七月*************************


'set lev 1000 500'


'define aqu7=ave(u.3,t=63,t=119)*ave(q.2,t=63,t=119)'
'define aqv7=ave(v.4,t=63,t=119)*ave(q.2,t=63,t=119)'


**************八月***********************
'set lev 1000 500'


'define aqu8=ave(u.3,t=64,t=120)*ave(q.2,t=64,t=120)'
'define aqv8=ave(v.4,t=64,t=120)*ave(q.2,t=64,t=120)'


************七、八月平均*********
'set lev 1000 500'
'define aqud=(aqu7+aqu8)/2'
'define aqvd=(aqv7+aqv8)/2'
******************水汽通量***************
'set lev 1000 '
'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)'
**********************************************
'set lon 20 120'
'set lat -20 60'
'set lev 1000 '
'set gxout shaded'
'd ashuiqi'
'set ccolor rainbow'
'run  cbarn.gs'


'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.'

黄鸿惠

试下先set ccolor ，再d 变量