一个有误的假相当位温的倾斜剖面图gs，求各位前辈指点！

南南偶卖糕Oo

各位前辈好，我是grads新手，最近在做一个关于假相当位温的倾斜剖面图，只需要2017年10月8日20时这一个时间的数据，但是我的图无法做出来，以下是我的gs文件，希望前辈们能帮我看看，指点一下！万分感谢！
'reinit'
'sdfopen G:\shuju\T_201710.nc'
'sdfopen G:\shuju\rhum_201710.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 lat 30 40'
'set lev 850 200'
'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'

'c'

river

你这个问题竟然还能回答。你去看一下发的斜剖面的帖子吧。你这个设置的顺序有问题，计算过程应该放在最前面，计算出来三维的数据，然后才能做剖面啊。你一上来就固定成一个点，就求了一条线上的数据，你让它怎么切出来一个面？

zaqbob

楼主最后解决了吗?也遇到相同的问题

南南偶卖糕Oo

zaqbob 发表于 2020-1-3 19:31
楼主最后解决了吗?也遇到相同的问题

没有呢，很抱歉帮不上你的忙，我最后换成其他图来说明问题了

river

几年以前的帖子突然诈尸了，吓死了

南南偶卖糕Oo

river 发表于 2022-4-20 12:55
几年以前的帖子突然诈尸了，吓死了

哈哈哈哈谢谢老师指点啦

茄子

南南偶卖糕Oo 发表于 2022-4-22 17:22
哈哈哈哈谢谢老师指点啦

请问怎么解决的呢

四叶草

river 发表于 2022-4-20 12:55
几年以前的帖子突然诈尸了，吓死了

哈哈，出来诈尸玩下

river

四叶草 发表于 2022-10-27 14:11
哈哈，出来诈尸玩下

今天把你诈出来了，哈哈。来，一起诈！