GrADS画固定经纬度某一时次位温垂直分布图

ミ青池の-ō

用ncep2.5*2.5资料画固定经纬度某一时次位温垂直分布图，
用‘d (air/10)*pow((1000/lev),0.286)’公式对吗？
从论坛里找到下面.gs的位温公式：
'define t=tmpprs'
'define rh=rhprs'
'define prs=lev'
'define es=(6.112*exp(17.67*(t-273.15)/(t-29.65)))'
'define q=rh*(0.62197*es/(prs-es))/100.0'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(t)-log(e)-4.805)'
'define theta=t*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
到'define e=prs*q/(0.62197+q)+1e-10' 出现错误，“1e”是什么东西？
应该用哪个公式是正确的？
我想画的是这个图，GrADS新手纠结了好久搞不定，求各位大神指教！

temperature

谢谢楼主 先看看

river

假相当位温和位温应该不一样吧，你查下书吧

oydl1985

楼主应该把错误提示放出来的呀~