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

river

晓文闯天下 发表于 2016-4-5 11:21
您好，我看了你整理的关于任意两点间斜剖图做法的总结，有一点疑问，想问问您。
   您的做法是将某一个点 ...

一个格点上的数据插值到斜线上十个点上面是什么原理？这个帖子里讲的东西和你说想法是一样的。只是利用矢量分解，把这个点上的值投影到斜线上的一个点，然后所有的点都投影到斜线上之后，再利用插值函数把这些点化成等值线······

晓文闯天下

river 发表于 2016-4-5 21:33
一个格点上的数据插值到斜线上十个点上面是什么原理？这个帖子里讲的东西和你说想法是一样的。只是利用矢 ...

脚本中，
'set x 1'
'set y 1'
'set lev 1000 100'
再加上默认的时间为第一个时刻，不是应该已经把格点确定为 1000-100 hPa上的一个（1,1）点了吗？
说的那个斜线上的十个点只是假设说插值过程中用到的斜线上的点。

晓文闯天下

river 发表于 2016-4-5 21:33
一个格点上的数据插值到斜线上十个点上面是什么原理？这个帖子里讲的东西和你说想法是一样的。只是利用矢 ...

脚本中，
'set x 1'
'set y 1'
'set lev 1000 100'
再加上默认的时间为第一个时刻，不是应该已经把格点确定为 1000-100 hPa上的一个（1,1）点了吗？
说的那个斜线上的十个点只是假设说插值过程中用到的斜线上的点。

river

晓文闯天下 发表于 2016-4-5 21:55
脚本中，
'set x 1'
'set y 1'

你没看见后面的循环？
while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(omega,'lon','lat')'
  lon = lon + 1

设置'set x 1'  'set y 1'在大部分情况下确实是你说的固定一个点，但是后面还有循环的，这时候就变成了表示一次插值一个点的意思了。就像你如果使用到ave函数求某个变量的时间平均，如果要出图的话，就一定要在使用ave函数的前面加一个set t 常数，那个常数是随便的，但是图是不会随着常数不同而变化的。你做这个也可以试试改变 x y的值，这个时候图应该是没有变化的

晓文闯天下

river 发表于 2016-4-6 07:41
你没看见后面的循环？
while (lon

受教了，麻烦楼主了

ymm9106

任意两点斜剖面，值得共享

ymm9106

thank you 楼主辛苦啦

river

晓文闯天下 发表于 2016-4-6 08:51
受教了，麻烦楼主了

不客气，帖子对你有用就是最好的

yuyuxuan

river 发表于 2016-4-6 16:41
不客气，帖子对你有用就是最好的

楼主  假如我想用FNL资料求假相当位温，那应该怎么弄呢，那个gr2stn函数应该如何使用呢 谢谢楼主

yuyuxuan

*'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))'
楼主你看这个是我求假相当位温的gs，我如果想画假相当位温的任意剖面图，我是不是要把这个里面每一个变量都用collect 1 gr2stn() 这么写一下，还是只需要插值最后算出的eqt然后画图呢？谢谢楼主