求助 求位温时出现的问题

heyufei

我弄了一个 求位温的GS  画出来的图象是这样子的
但是我还是有一个问题 能不能把等温线的间隔弄得小一点 这个元贝的间距是5度 想要弄成2度的 不然没法说明那个区域 位温异常啊

附上 Gs 求指教啊
'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'

传说中的谁

cint或者clevs

heyufei

传说中的谁 发表于 2014-5-8 12:37
cint或者clevs

亲 能不能详细地说一下

heyufei

传说中的谁 发表于 2014-5-8 12:37
cint或者clevs

就是说 在什么什么句子和什么什么句子之间 插一句什么什么的语句就去

一片落叶惊秋风

有两个办法，一个而是用‘set cint ’
一个是用'set zlog on'
'set ylevs 1000 925 850 700 600 500 400 300 200 100'

want_data

你的gs中的theta怎么和人家的假相当位温一样？
位温的公式不是这样的吧

19921219

这个是假相当位温？？

桔梗花开

这个位温的代码是经典的算法吗？我看挺多人是这么算的，有一篇挺老的文献是这样写的公式，但是里面是相对湿度不是比湿

周兰810626

我也觉得你的gs中的theta怎么和人家的假相当位温一样？确定是位温吗？正需要位温的