求助 求位温

heyufei

大家能不能给看看 我这个求116经度下剖面的位温GS有什么问题 为什么只能出现500HPA 到 100hpa的

'reinit'
'open E:\aa\fnl\fnl_20130705_12_00.ctl'
'set lon 100 125'
'set lat 25 45'
'set lev 1000 100'
'set mpdset cnworld'
'set grid off'
'set display color white'
'c'
'set grads off'

'set csmooth on'

'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)))'


'set lat 25 45'
'set lon 116'
'set lev 1000 100'
'set clevs 330 335 340 345 350 355 360 365 370 375 380'
'set ccols 14 9 4 11 5 10 7 12 8 2 6'


'd theta'
'cbarn'
'draw title 23-18-116E-2theta'

'printim E:\aa\fnl\23-18-116E-2theta.jpg'
'print'

传说中的谁

这不是计算的问题，是绘图的问题。
请问楼主知道set clevs的作用吗？如果知道，那您应该能弄明白问题出在哪；如果不知道，请弄清楚再去使用它。

river

你把这两句去掉'set clevs 330 335 340 345 350 355 360 365 370 375 380'
'set ccols 14 9 4 11 5 10 7 12 8 2 6'

heyufei

谢谢 这样就做出来了 E:/aa/fnl/theta/06-06-115-theta

heyufei

还有一个问题 想问一下  我想怎么把这张图的 间距变小一点呢 现在每条等温线的差值是5 想变成 2
'reinit'
'open E:\aa\fnl\fnl_20130706_12_00.ctl'
'set lon 100 125'
'set lat 20 45'
'set lev 1000 300'
'set mpdset cnworld'
'set grid off'
'set display color white'
'c'
'set grads off'

'set csmooth on'

'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.'
'define e=prs*q/(0.62197+q)+1e-10'
'define theta=t*pow((1000/prs),(0.2854*(1.0-0.28*q)))'


'set lat 31'
'set lon 100 120'
'set lev 1000 100'


'd theta'
'cbarn'
'draw title 06-12-31N-theta'

'printim E:\aa\fnl\theta\06-12-31N-theta.jpg'
'print'

xry_2005

set cint 2

19921219

是的，可以set cint 2

19921219

请问这个是位温还是假相当位温？

19921219

请问是位温还是假相当位温？

krysunder

这句“'define e=prs*q/(0.62197+q)+1e-10' ” 有什么用？