水汽通量散度图，能运行，图上空白，求大神指点

玮ya

'reinit'
'sdfopen d:\aa\uwnd.mon.mean.nc'
'sdfopen d:\aa\vwnd.mon.mean.nc'
'sdfopen d:\aa\rhum.mon.mean.nc'
'sdfopen d:\aa\air.mon.mean.nc'
'enable print d:\aa\gmf\sd.gmf'
'set lon 70 140'
'set lat 15 55'
'set lev 850'
i=597
while(i<=599)
'set t 'i
'set grads off'
'define prs=lev'
'define g=9.8'
'define es=(6.112*exp((17.67*(air.4-273.15))/(air.4-29.65)))'
'define qs=(0.62197*es/(prs-0.378*es))'
'define q=qs*rhum.3/100'
'define x=q*uwnd.1'
'define y=q*vwnd.2'
'define sqdiv=hdivg(x,y)*1e9'
'set gxout shaded'
'd sqdiv'
'print'
'c'
i=i+1
endwhile
'disable print'
;

小傅之追梦

运行过程中有没有显示什么异常？

玮ya

小傅之追梦 发表于 2013-5-11 22:31

                               
登录/注册后可看大图

运行过程中有没有显示什么异常？

没有，源文件也打开过都是好的

river

看着没什么问题哎

逸儿~~

玮ya 发表于 2013-5-11 22:39

                               
登录/注册后可看大图

没有，源文件也打开过都是好的

我来凑个热闹~~傅兄的论文做好没？

kongfeng0824

报错信息是什么？

小傅之追梦

逸儿~~ 发表于 2013-5-12 10:31

                               
登录/注册后可看大图

我来凑个热闹~~傅兄的论文做好没？

还在等数据，急急急急急急死我了

玮ya

逸儿~~ 发表于 2013-5-12 10:31

                               
登录/注册后可看大图

我来凑个热闹~~傅兄的论文做好没？

我真是服了你了

玮ya

小傅之追梦 发表于 2013-5-12 13:41

                               
登录/注册后可看大图

还在等数据，急急急急急急死我了

导师催得紧，我也灰常急啊

xiaoxiao2013

'reinit'
'sdfopen f:\thesis\data\levels.nc'
'sdfopen f:\thesis\data\output.nc'

'set lon 90 115'
'set lat 25 45'
'set lev 850'

id=109
while(id<=118)
'set grads off'
'set parea 1 10.5 1 7.5'
'set xlopts 1 5 0.12'
'set ylopts 1 5 0.12'
'set t 'id''
'define prs=lev'
'define g=9.8'
'define es=(6.112*exp((17.67*(t.2-273.15))/(t.2-35.86)))'
'define qs=(0.62197*es/(prs-0.378*es))'
'define q=qs*r.1/100'
'define x=q*u.1'
'define y=q*v.1'
'define sqdiv=hdivg(x,y)*1e9'

'set gxout contour'
'set cthick 5'
'set ccolor 4'
'set mpdset cnriver'
'set map 2 1 7'
'draw map'
'd sqdiv'
'set mpdset mres'
'set mpdset cnworld'
'set map 1 1 6'
'draw map'

'q time'
res=subwrd(result,3)
res=substr(res,1,12)
'draw title 850hPa  divergence of moisture flux on 'res''

'enable print f:\thesis\data\850\'res'sqdivg.gmf'
'print'
'disable print'
pull dummy
'c'
id=id+1
endwhile
我的是这个样子的，能出图，你看看吧！