斜剖面不能画计算量大点的物理量吗？

qxcong201611

{:5_213:}{:5_213:}

毛毛虫106

顶{:eb502:}{:eb502:}{:eb502:}

毛毛虫106

金钱不够，评论来凑，帮楼主顶

a930013

感谢大大分享

ITCZ

太好了，终于找到了，谢谢楼主分享

南南偶卖糕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'