水汽通量和散度剖面叠加图，error form HDIVG:invalid dimension environment

Syue

写好的gs，用grads运行了一下报错了，说是错误出现在27行？没看出来是什么情况，大家能不能帮帮忙看是怎么回事？以下是我写好的gs：
'open e:\20151105\qfdiv\q1.ctl'
'open e:\20151105\qfdiv\q2.ctl'
'open e:\20151105\fnl_20151104_00_00.ctl'
'set lon 101 120'
'set lat 23'
'set z 1 17'
'set t  1'
'set grid off'
'set grads off'

'set lon 101 120'
'set lat 23'
'set z 1 17'
'set t  1'


'define tt=TMPprs.3'
'define es=6.112*exp(17.67*(tt-273.15)/(tt-29.65))'
'define qs=0.62197*es/(lev-0.378*es)*1000'
'define q=qs*RHprs.3/100.0'
'define q1=q*ugrdprs.3/9.8'
'define q2=q*vgrdprs.3/9.8'
'define qfdiv=hdivg(q1,q2)*1e6'


* 'set lon 110'
*'set lat 11 40'
*'set lev 1000 100'
*'set clopts -1 -1 0.15'
* 'set xlopts 1 4 0.20'
*'set ylopts 1 4 0.20'
*'set rgb 16 228 228 228'
*'set rgb 17 216 216 216'
*'set rgb 18 200 200 200'
*'set rgb 19 175 175 175'
*set rgb 20 150 150 150'
*'set rgb 21 125 125 125'
*'set rgb 22 100 100 100'
*'set rgb 23 75 75 75'
*'set rgb 24 50 50 50'
*'set rgb 25 25 25 25'
*'set rgb 26 0 0 0'
*'set clevs  2  4  6  8  10'
*'set ccols 16  18   21 22 23  25 '


  *'set cint 2'
   *'d qflux'
   *'cbarn  0.5  0.5 '

*'set gxout barb'
*'d ugrdprs;vgrdprs'


'set gxout contour'
'set ccolor 0'
'set cstyle 1'
'set cthick 4'
'set csmooth on'
'set cint 4'
'd qfdiv'

'set lon 101 120'
'set lat 23'
'set z 1 17'
'set t  1'
'set gxout  vector'
'd q1.1;q2.2 '


'printim e:\20151105\qfdiv_qflux.gif white'
'disable print'
'reinit'

男紫汗

'set lon 101 120'
'set lat 23'计算的时候维度不要减，等画图的时候再减

river

'set lat 23'
你这都已经固定到一条线上去了，我问你怎么插值？

Syue

男紫汗 发表于 2017-6-20 08:49
'set lon 101 120'
'set lat 23'计算的时候维度不要减，等画图的时候再减

我把'set lat 23'改为'set lat 11 40'
然后后面倒数第三段加了'set lat 23'
                                 'set gxout contour'
                                 'set ccolor 0'
                                 'set cstyle 1'
                                 'set thick 4'
                                 'set csmooth on'
                                 'set cint 4'
                                 'd qfdiv'
运行后红字没了，但是错误还是在第27行，怎么回事呀？

Syue

river 发表于 2017-6-20 14:55
'set lat 23'
你这都已经固定到一条线上去了，我问你怎么插值？

'set lat 23'改成了'set lat 11 40'
倒数第三段加了'set lat 23'沿着北纬23画剖面图，运行后还是错误在27行。river说的插值在哪里啊？

男紫汗

Syue 发表于 2017-6-20 15:53
我把'set lat 23'改为'set lat 11 40'
然后后面倒数第三段加了'set lat 23'
                          ...

第27行是哪一行呢？我觉得你这个图有问题，散度是水平二维场，然后又要画纬向水汽通量剖面图，纵坐标不一样的，要么两个纵坐标分别画左边和右边。。。具体我也没画过

river

Syue 发表于 2017-6-20 15:58
'set lat 23'改成了'set lat 11 40'
倒数第三段加了'set lat 23'沿着北纬23画剖面图，运行后还是错误在2 ...

你都修改过了，我怎么知道第二十七行是哪一行？至于插值在哪不知道，请你去看hdivg的说明

yh65210226

谢谢分享，很有用。

惠子

谢谢分享，很有用。

carlye2002

谢谢分享，很有用。