任意两点间斜剖面做法总结【利用GRADS自带函数】

提拉米苏lei

折腾不同分辨率的地形剖面，楼主的学习贴需要好好学习，感谢分享！可惜还没到级别无法下载（自动回复：请不要使用迅雷等下载工具，点我查看下载帮助）……

提拉米苏lei

另外，想请教下，如果用ncep数据中的HGTsfc和HGTprs两组数据求差值，正值的地方是不是也代表地形呢？

提拉米苏lei

扣了金钱却没有下载下来怎么回事……

南南偶卖糕Oo

最近在学习画这种倾斜的剖面图，想楼主帮我看看我的程序，我的图一直出不出来，但是我不知道我哪里错了！！！万分感谢！！！
以下是我的gs:
'reinit'
'sdfopen G:\shuju\T_201710(duo).nc'
'sdfopen G:\shuju\rhum_201710(duo).nc'
'set display white'
'c'

'set grads off'
'set grid off'
'set zlog on'
'set x 1'
'set y 1'
'set lev 850 200'
'set time 12Z08OCT2017'

'define t=T.1'

'define rh=r.2'

'define prs=lev'

'define es=(6.1078*exp(17.27*(t-273.15)/(t-35.86)))'

'define q=rh*(0.62197*es/(prs-es))/100.'

'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 eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'

lon1 = 102.0
lon2 = 107.0
lat1 = 40.0
lat2 = 30.0
lon = lon1
'collect 1 free'
'collect 2 free'
while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1)/(lon2-lon1)
  'collect 1 gr2stn(eqt,'lon','lat')'
  lon = lon + 1
endwhile

'set map 15 1 9'
'set csmooth on'
'set lon 102 107'
'set lev 850 200'
'set time 12Z08OCT2017'
'set ylevs 850 800 700 600 500 400 300 200'
'set xlabs 40N,102E|38N,103E|36N,104E|34N,105E|32N,106E|30N,107E|'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set gxout contour'

'set csmooth on'
'set cthick 5'
'set cint 20'
'set ccolor 1'
'set clopts 0 6 0.15'
'd coll2gr(1,eqt)'
'printim G:\shuju\Z_new.png white'

'c'

南南偶卖糕Oo

yuyuxuan 发表于 2016-4-14 11:15
又来麻烦你了，我画FNL1°×1°自带的数据，比如温度，垂直速度我已经可以画出来了，但是假相当位温还是 ...

楼主的假相当位温gs倾斜剖面解决了吗？可以帮我看看我的gs吗？万分感谢你！
'reinit'
'sdfopen G:\shuju\T_201710(duo).nc'
'sdfopen G:\shuju\rhum_201710(duo).nc'
'set display white'
'c'

'set grads off'
'set grid off'
'set zlog on'
'set x 1'
'set y 1'
'set lev 850 200'
'set time 12Z08OCT2017'

lon1 = 102.0
lon2 = 107.0
lat1 = 40.0
lat2 = 30.0
lon = lon1
'collect 1 free'
'collect 2 free'

'define t=T.1'

'define rh=r.2'

'define prs=lev'

'define es=6.1078*exp(17.27*(t-273.15)/(t-35.86)))'

'define q=rh*(0.62197*es/(prs-es))/100.'

'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 eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'

while (lon <= lon2)
lat = lat1 + (lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(eqt,'lon','lat')'
lon = lon + 1
endwhile

'set map 15 1 9'
'set csmooth on'
'set lon 102 107'
'set lev 850 200'

'set time 12Z08OCT2017'
'set ylevs 850 800 700 600 500 400 300 200'
'set xlabs 40N,102E|38N,103E|36N,104E|34N,105E|32N,106E|30N,107E|'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set gxout contour'

'set csmooth on'
'set cthick 5'
'set cint 10'
'set ccolor 1'
'set clopts 0 6 0.15'
'd coll2gr(1,-u)'
'printim G:\shuju\Z_new.png white'

'c'

yuyuxuan

南南偶卖糕Oo 发表于 2018-11-5 11:07
楼主的假相当位温gs倾斜剖面解决了吗？可以帮我看看我的gs吗？万分感谢你！
'reinit'
'sdfopen G:\shuj ...

对不起啊，后来我没用grads画，用的MATLAB，抱歉了，还是两年前的呢

南南偶卖糕Oo

yuyuxuan 发表于 2018-11-5 19:01
对不起啊，后来我没用grads画，用的MATLAB，抱歉了，还是两年前的呢

哦哦，无论如何感谢你回复我哈~

yuyuxuan

南南偶卖糕Oo 发表于 2018-11-5 21:43
哦哦，无论如何感谢你回复我哈~

当时我是画出来了，但图不对，也不报错，就放弃了

最后一个匈奴

感谢楼主分享！！！

zyr

把贴子所有的楼层都翻了过来，终于找到自己的原因，并出图了。感谢楼主！