幂函数的问题

443753891

用GRADS算幂函数不报错，但画出来的图是空白，没值。请高手帮看看，幂函数写得对不对，哪里错了。谢谢
'reinit'
'open F:\20090728\wrf\wrfout\outd02.ctl'
'enable print F:\20090728\wrf\amk\amk.gmf'
i=25
while(i<=31)
'set mpdset cnworld'
'set lat 15 50'
'set lon 105 135'
'set grid off'
'set grads off'
'set xlint -1'
*'set ylint -1'
'set lev 850'
'set t 'i
'define l=2500'
'define cp=1004'
'define g=9.8'
'define es=6.1078*exp(17.2693882*(tk-273.16)/(tk-35.86))'
'define qs=622*es/(lev-0.378*es)'
'define q=rh*qs/100'
'define te=tk+l*q/cp'
'define thete=te*pow(1000/lev,0.2857)'
'define pr=-16*thete+1805'
'define amk=cp*(1-pow(pr/lev,0.2857))*(tk+l*q/cp)/9.8'
'set lon 115 125'
'set lat 25 35'
'set lev 850'
'd amk''
'print'
'c'
i=i+1
endwhile
'disable print'
;

兰溪之水

'define thete=te*pow(1000./lev,0.2857)'
'define pr=-16*thete+1805'
'define amk=cp*(1-pow(1.*pr/lev,0.2857))*(tk+l*q/cp)/9.8'

试试

443753891

兰溪之水 发表于 2012-12-18 09:52

                               
登录/注册后可看大图

'define thete=te*pow(1000./lev,0.2857)'
'define pr=-16*thete+1805'
'define amk=cp*(1-pow(1.*pr/l ...

还是不行，前一个POW可以算出来，但加上第二个POW就不行，通过分步算，好像是与第二个POW中的pr/lev有关

兰溪之水

443753891 发表于 2012-12-18 10:03

                               
登录/注册后可看大图

还是不行，前一个POW可以算出来，但加上第二个POW就不行，通过分步算，好像是与第二个POW中的pr/lev有关

检查数和数据格式之类的，自个试试吧~！

river

好像没什么问题，我用NC资料可以画出来
不过没太研究算出来的值对不对

443753891

river 发表于 2012-12-18 12:12

                               
登录/注册后可看大图

好像没什么问题，我用NC资料可以画出来
不过没太研究算出来的值对不对

噢，我用的WRF模拟ARW后处理资料。

river

443753891 发表于 2012-12-18 12:26

                               
登录/注册后可看大图

噢，我用的WRF模拟ARW后处理资料。

就是说我用nc可以画出来，那就说明是你的资料的问题吧，不是gs的问题

kongfeng0824

兰溪之水 发表于 2012-12-18 09:52

                               
登录/注册后可看大图

'define thete=te*pow(1000./lev,0.2857)'
'define pr=-16*thete+1805'
'define amk=cp*(1-pow(1.*pr/l ...

我用兰溪的办法试了一下是可以的。楼主可以再仔细检查一下