新学grads不会用，帮忙看一下，画不出图

天之边

'reinit'
'sdfopen e:\air.2002.nc'
'sdfopen e:\rhum.2002.nc'
'set X 1'
'set Y 1'
'set t 1 1459'
'set grads off'
'set grid on'
'set lev 850'
'define T850=air.1-273.16'
'define Td850=T850-((14.55+0.114*T850)*(1-0.01*rhum.2) + pow((2.5+0.007*T850)*(1-0.01*rhum.2),3) + (15.9+0.117*T850)*pow((1-0.01*rhum.2),14))'
'set lev 700'
'define T700=air.1-273.16'
'define Td700=T700-((14.55+0.114*T700)*(1-0.01*rhum.2) + pow((2.5+0.007*T700)*(1-0.01*rhum.2),3) + (15.9+0.117*T700)*pow((1-0.01*rhum.2),14))'
'set lev 500'
'define T500=air.1-273.16'
'define K=T850-T500+Td850-(T700-Td700)'
'Display  tloop(aave('K',lon=110,lon=117,lat=35,lat=40))'

lqouc

感觉是空间维有问题，你把报错截图上来。

天之边

lqouc 发表于 2013-11-25 10:44

                               
登录/注册后可看大图

感觉是空间维有问题，你把报错截图上来。

开始的时候是 'set lon 110 117' 'set lat 35 40',没有报错，也没有图，现在报错是

天之边

天之边 发表于 2013-11-25 10:48

                               
登录/注册后可看大图

开始的时候是 'set lon 110 117' 'set lat 35 40',没有报错，也没有图，现在报错是

cannot plot data

zzl

你先看一下'd k' 能不能出图

mofangbao

新学就拿别人做过的东西先做一遍  这种代码不太适合新学的时候做 最起码要对GrADS有一定了解 要不然即使你算出结果  你也不知道过程是否正确

lqouc

天之边 发表于 2013-11-25 10:48

                               
登录/注册后可看大图

cannot plot data

就这一句报错？话说你的k为什么要加引号？？？这个不对吧。