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

yuyuxuan

river 发表于 2016-4-14 07:53
函数本身是grads自带的，不需要添加

又来麻烦你了，我画FNL1°×1°自带的数据，比如温度，垂直速度我已经可以画出来了，但是假相当位温还是有问题，运行完gs没有错误 但是图是这样的F:/G/TU.png
'reinit'
'open F:/G/FNL/test1/fnl_20121101_00_00.ctl'
*t=1
*while(t<=20)
*'clear'
*'set t 't
*'q time'
*m_time=subwrd(result,3)
'set t 2'
'set grads off'
'set grid off'
'set zlog on'
'set csmooth on'
'set ylint 100'
'set xlint 2'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set x 1'
'set y 1'
'set lev 1000 100'

lon1=116
lon2=125
lat1=37
lat2=42
lon=lon1
'collect 1 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)

'define tmp=TMPprs'
'define rh=RHprs'
'define prs=lev'
'define es=(6.11*exp(17.27*(tmp-273.16)/(tmp-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(tmp)-log(e)-4.805)'
'define theta=tmp*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'

'collect 1 gr2stn(eqt,'lon','lat')'
lon=lon+1
endwhile
'set lon 116 125'
*'set xaxis 116 125 2'
'set clab on'
'set gxout contour'
'd coll2gr(1,-u)'
*'draw title fnl,'m_time
*'gxprint F:/G/FNL/test1/jxdww/'t'.png white'
*t=t+1
*endwhile

这个是我的gs

yuyuxuan

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

file:///F:/G/TU.png

river

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

先计算出假相当位温，然后再插值

yuyuxuan

river 发表于 2016-4-14 16:13
先计算出假相当位温，然后再插值

好的谢谢楼主，还有就是，楼主下面的例子画剖面风场，楼主的剖面是东北-西南走向，U、V正向和剖面都在第一象限，而如果剖面走向是西北-东南走向的话，是不是要在下面这个式子，U前面加一个负号呢？谢谢楼主
'collect 2 gr2stn(Uwnd.2*cos(alfa)+Vwnd.3*sin(alfa),'lon','lat')'

yuyuxuan

river 发表于 2016-4-14 16:13
先计算出假相当位温，然后再插值

我想明白了楼主，不用改的，直接那个公式就也可以求风，再次感谢

river

yuyuxuan 发表于 2016-4-14 17:05
我想明白了楼主，不用改的，直接那个公式就也可以求风，再次感谢

想明白就好了

ChrisLee

好东西学习一下

haddad

好帖，楼主好人的

熊猫二货

楼主为啥我用grads打不开你给的地形数据

river

熊猫二货 发表于 2016-4-25 15:32
楼主为啥我用grads打不开你给的地形数据

怎么打不开了，什么错误提示，问题要说明白