'reinit'
'open F:\grads\2012022100.ctl'
'set lon 60 100'
'set lat 30 60'
'set lev 850'
*'set grads off'
*'set grid off'
'define tt=tmpprs' 昨天这里出错了,变量没写对,引号也错啦
'define es=6.11*exp(17.26*(tt-273.16)/(tt-35.86))'
'define qs=0.622*es/(lev-0.378*es)'
'define qq=rhprs/100.*qs'
'define xflux=qq*ugrdprs'
'define yflux=qq*vgrdprs'
'define sqdiv=hdivg(xflux,yflux)*1e8'
'd sqdiv'