求助！fnl资料做水汽通量散度随时间变化剖面图的定义问题！！！

cat1106123

新学grads不久，需要做水汽通量散度随时间变化的剖面图，只有fnl的资料，在论坛里搜了一些类似剖面图的gs,大多是nc的资料做的，按照fnl的变量改动了一下，发现出不了图，显示如下：


这是什么意思呢？急盼大神解答啊！！！

gs如下：
function b( args  )
dat = subwrd( args,1)
'open k:\grib\1\fnl6.ctl'
'set t 1 48'
'set lev 1000 300'
'define prs=lev'
'set ccolor 1'
'set grid off'
'define  l1=TMPprs'
if(l1>263.0);
'define a1=0.622*6.11*EXP(17.26*(TMPprs-273.16)/(TMPprs-35.86))'
'define b1=prs-0.278*EXP(17.26*(TMPprs-273.16)/(TMPprs-35.86))'
'define  qs1=a1/b1'
endif
if(l1<=263.0);
'define a1=0.622*6.11*EXP(21.87*(TMPprs-273.16)/(TMPprs-7.66))'
'define b1=lev-0.278*EXP(21.87*(TMPprs-273.16)/(TMPprs-7.66))'
'define  qs1=a1/b1'
endif
'define  q1=qs1*RHprs'
'define qdiv1 = q1 * hdivg(Ugrdprs,Vgrdprs)/9.8'
'set grads off'
'set grid off'
'set csmooth on'
'set map 15 1 9'
'set parea 1 10.2 1.5 8'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'
'set lon 126'
'set lat 43'
'set zlog on'
'set gxout shaded'
'd qdiv1'
'run cbar.gs'
'set ylpos 0 r'
'set gxout contour'
'set cthick 4'
*'set clskip 1'
'd qdiv1'
'draw title qdiv change'
'printim k:\grib\qdiv-change.gif gif white'
*'reinit'
;

cat1106123

gs瞎凑的比较乱，实在是不太懂，希望懂的帮帮我啊~

传说中的谁

function b( args  )
dat = subwrd( args,1)
楼主可以告诉我这两句有什么作用吗？看上去跟你后面的内容没有任何关联性。

cat1106123

传说中的谁 发表于 2013-10-29 09:16

                               
登录/注册后可看大图

function b( args  )
dat = subwrd( args,1)
楼主可以告诉我这两句有什么作用吗？看上去跟你后面的内容没 ...

这个......是从论坛里的一个讲制作“假相当位温剖面”的gs文件里copy下来的，我也不清楚呃，呵呵，因为发现那个gs里有比湿的计算，还有剖面的设置，所以就直接把计算部分改动了一下......