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

river

yuyuxuan 发表于 2016-4-13 16:08
*'define es=(6.11*exp(17.27*(tmp-273.16)/(tmp-35.86)))'
*'define q=rh*(0.62197*es/(prs-es))/100.'
...

那还用说啊，肯定是只插值你要的最后的那个变量啊

yuyuxuan

river 发表于 2016-4-13 21:24
那还用说啊，肯定是只插值你要的最后的那个变量啊

谢谢谢谢您，我这就去看看，还有还有，这两句话是啥意思嘞，
'set x 14 16'
'set xaxis 'lon1' 'lon2
为啥您下面的那个例子 就没有设置这两句话嘞，麻烦您了

yuyuxuan

river 发表于 2016-4-13 21:24
那还用说啊，肯定是只插值你要的最后的那个变量啊

还有就是 楼主  调用那几个函数  需不需要额外安装什么补丁呀，就像兰溪发的那个处理FNL资料要安装g2ctl.exe那些东西一样。谢谢楼主。

river

yuyuxuan 发表于 2016-4-13 21:47
谢谢谢谢您，我这就去看看，还有还有，这两句话是啥意思嘞，
'set x 14 16'
'set xaxis 'lon1' 'lon2

这不是麻烦不麻烦的问题，这是基础问题，建议你去看grads实用手册，论坛就能下载到

river

yuyuxuan 发表于 2016-4-13 21:59
还有就是 楼主  调用那几个函数  需不需要额外安装什么补丁呀，就像兰溪发的那个处理FNL资料要安装g2ctl. ...

函数本身是grads自带的，不需要添加

yuyuxuan

cathy 发表于 2013-11-22 20:12
太感谢了！楼主这么好的分享！我用我的数据试了一下，有点问题，不知道什么原因，求大神指点~~~
这是我照 ...

您好   我想问一下  'set x 213 232'
这句话是什么意思嘞

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')'