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

river

老A 发表于 2022-7-29 00:17
谢谢楼主回复，我是用楼主的地形和代码画的，可能是确实精度不够，不过那个区域确实是大部分在海拔5000m ...

orog

你看一下，可能确实是精度不够，海拔超过5000米的地方就那一点点，基本都在 30-35，85-90 之间

老A

谢谢楼主，确实是精度不够，超过5000m的看样子是在青藏高原，青藏高原西北部那边的喀喇昆仑山区域海拔高于5000m的没有显示出来，非常感谢您！

茄子

'reinit'
'sdfopen d:\by\xt.nc'
'enable print d:\by\wwpm-6-27-20.gmf'
'set grads off'
'set grid off'
'set mpdset hires'
'set mpdset  cnworld'
'set lon 105 115'
'set lat 28 34'
'set lev 1000 150'
'set t 13'
'define prs=lev'
'define mm=r'
'define e1=17.27*(t-273.16)/(t-35.86)'
'define es=6.112*exp(e1)'
'define qs=mm*(0.62197*es/(prs-0.378*es))/100'
'define e=prs*q/(0.62197+q)+1e-10'
'define tk=55.0+2840.0/(3.5*log(t)-log(e)-4.805)'
'define pot=t*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=pot*exp(((3376./tk)-2.54)*q*(1.0+0.81*q))'
lon1 = 105.0
lon2 = 115.0
lat1 = 40.0
lat2 = 30.0
lon = lon1
'collect 1 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(eqt,'lon','lat')'
lon=lon+1
endwhile
'set lon 105 115'
'set lev 1000 150'
*'set xlabs 34N,108E|33N,109E|32N,110E|31N,111E|30N,112E|29N,113E|28N,114E|'
'set gxout contour'
'd coll2gr(1,-u)'
'print'
'disable print'
;

上面是求假相当位温，任意两点剖面，运行后有问题，请大家帮忙看看。

PIZI

这个绝对是经典的经典