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

river

BLUEANGEL2012 发表于 2016-3-14 10:46
想问楼主的地形文件层数是和数据文件一样层数么？  还有每层的数据都一样的吗  还是从第二层开始可以赋值 ...

我帖子里的数据和地形数据都是17层的，具体资料你可以直接用grads打开地形ctl 画图看看，从500百帕以上基本都是负值了，也就是再往上没有地形了

wwwaaasss

{:5_214:}{:5_214:}

晓文闯天下

感谢楼主

晓文闯天下

您好，我看了你整理的关于任意两点间斜剖图做法的总结，有一点疑问，想问问您。
   您的做法是将某一个点的数值差值到我所做的斜剖面的直线上，（间隔一度或者是是其他），假设这个线上有10个相对应的经纬度，那也就是你所选格点上那一个值插值到斜线上的十个点上面。而我认为沿斜线的剖面图是这条斜线上的十个经纬度对应的附近的十个格点的数值插值到其分别对应的经纬度上面得到的数据，而不是一个格点差值到十个经纬度上所得到的数据。
   这是我的一点疑问，不知道理解的对不对，还请楼主帮忙分析一下。麻烦了，谢谢。

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

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