- 积分
- 7756
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-10-28
- 最后登录
- 1970-1-1
|
发表于 2012-7-6 20:57:41
|
显示全部楼层
本帖最后由 njzqxt 于 2012-7-6 21:01 编辑
请XD们帮着看看,我作的图有什么问题。谢谢了。
'reinit'
'sdfopen e:\nc\day\uwnd.2012.nc'
'sdfopen e:\nc\day\vwnd.2012.nc'
'sdfopen e:\nc\day\rhum.2012.nc'
'sdfopen e:\nc\day\air.2012.nc'
'set lon 90 110'
'set lat 18 33'
'set lev 850'
'set t 151'
'set grads off'
'set grid off'
'set xlopts 1 1 0.15'
'set ylopts 1 1 0.15'
'set clopts 1 1 0.15'
'set xlint 3'
'set ylint 3'
'define tt=air.4+273.16'
'define es=6.11*exp(17.26*(air.4-273.16)/(tt-35.86))'
'define qs=0.622*es/(lev-0.378*es)'
'define qq=qs*rhum.3/100'
'define g=9.8'
'define xflux=qq*uwnd.1/g'
'define yflux=qq*vwnd.2/g'
'define sqdiv=hdivg(xflux,yflux)*1e8'
'd sqdiv'
'run d:\grads19\map\yn.gs'
'printim e:\201205yjks\sqtl.gif'
'print e:\201205yjks\sqtl.gmf'
'print'
'disable print'
; |
|