想画单站水汽通量散度随时间和高度的变化图

筱筱筱筱筱雅痴

以下是我的脚本，希望好心人指导一下，谢谢
'reinit'
'open J:\lw-data\grads\chuli\fnl_20161230_12_00.ctl'
*'set lon 100 130'
*'set lat 30 50'
'set lev 1000 500'
'set t 3 10'
'define prs=lev'
'define g=9.8'
'define es=(6.112*exp((17.67*(TMPprs-273.15))/(TMPprs-29.65)))'
'define qs=(621.97*es/(prs-0.378*es))'
'define q=qs*rhprs/100'
'define Sqtl=mag(q*UGRDprs/g,q*VGRDprs/g)'
'define SqtlSD=hdivg(q*UGRDprs/g,q*VGRDprs/g)/100'
'set mpdset cnworld'
'set grads off'
'set grid off'
'set csmooth on'
'set gxout shaded'
'd SqtlSD*1e7(lon=117.12,lat=39.04)'
'cbarn 0.7 0'
'set gxout contour'
'd SqtlSD*1e7(lon=117.12,lat=39.04)'
'printim J:\lw-data\grads\gs\totalsqsd.png white'
;