- 积分
- 44
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-9-20
- 最后登录
- 1970-1-1

|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
写了个关于水汽通量散度剖面图的代码,但是总是出错。错误直接上图,求解答。
代码如下:
'reinit'
'set grid off'
'set timelab off'
'set grads off'
'sdfopen e:\201409\air.2014.nc'
'sdfopen e:\201409\rhum.2014.nc'
'sdfopen e:\201409\shum.2014.nc'
'sdfopen e:\201409\uwnd.2014.nc'
'sdfopen e:\201409\vwnd.2014.nc'
'set t 992 1044'
'set lev 1000 300'
'define p=lev'
'define e=(6.112*exp((17.67*(air.1-273.16))/(air.1-29.65)))'
'define qs=(0.622*e/(p-0.378*e))'
'define q=qs*rhum.2/100'
'define qdiv=hdivg(q*uwnd.4/9.8,q*vwnd.5/9.8)'
'set lon 108'
'set lat 34'
'd qdiv'
'enable print e:\201409\1\qdiv.gmf'
'set zlog on'
'set ylevs 1000 925 850 700 600 500 400 300 250 200 150 100'
'set gxout shaded'
'set clab on '
'set cint 0.3'
'set rbcols auto'
'set csmooth on'
'run cbar.gs'
'set ylpos 0 r'
'set gxout contour'
'set cint 0.3'
'set cthick 11'
'd qdiv'
'draw title \sep.6-18 qdiv'
'draw map'
'draw xlab \time'
'draw ylab \presure'
'print'
'disable print'
'printim e:\201409\1\qdiv1.png white'
'c'
|
|