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

cathy

太感谢了！楼主这么好的分享！我用我的数据试了一下，有点问题，不知道什么原因，求大神指点~~~
这是我照着写的gs：
'reinit'
'sdfopen E:\temp.nc'
'enable print E:\oblique.gmf'

'set x 1'
'set y 1'
'set z 1 26'
'set t 121'

lon1 = -147.0
lon2 = -132.0
lat1 = 0.0
lat2 = 5.0
lon = lon1
'collect 1 free'
while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(temp_a,'lon','lat')'
  lon = lon + 1
endwhile

'set vpage 0.25 10.75 0.2 8.08'
'set parea 0.5 10.5 1.02 7.58'
'set grads off'
'set grid off'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'

'set x 213 232'
'set xaxis 'lon1' 'lon2
'set clab on'
'set csmooth on'
'set gxout shaded'
'd coll2gr(1,-u)-273.15'
'set gxout contour'
*'set cint 1'
'set cthick 11'
'd coll2gr(1,-u)-273.15'
'draw title pottmp along (0,147W) to (5S,132W)'
'print'
'disable print'
;
出现的错误提示为：

river

cathy 发表于 2013-11-22 20:12

                               
登录/注册后可看大图

太感谢了！楼主这么好的分享！我用我的数据试了一下，有点问题，不知道什么原因，求大神指点~~~
这是我照着 ...

gs看不出问题。你的nc资料有26层吗？

cathy

我的nc数据一共有28层，我只用了其中26层，不然我把我的nc资料ctl贴出来一下吧，我也很奇怪

river

cathy 发表于 2013-11-22 22:27

                               
登录/注册后可看大图

我的nc数据一共有28层，我只用了其中26层，不然我把我的nc资料ctl贴出来一下吧，我也很奇怪

你的资料经度是从121.5开始的，你怎么写的都是负的呢lon1 = -147.0  lon2 = -132.0

cathy

river 发表于 2013-11-23 06:23

                               
登录/注册后可看大图

你的资料经度是从121.5开始的，你怎么写的都是负的呢lon1 = -147.0  lon2 = -132.0

非常感谢哈，我怎么就没看出来呢，可是我改为lon1 = 207.5 lon2 = 228.5之后还是有错呢
还有就是原来的gs中'd coll2gr(1,-u)-273.15' 写错了，应该是'd coll2gr(1,-u)'。

river

cathy 发表于 2013-11-23 08:54

                               
登录/注册后可看大图

非常感谢哈，我怎么就没看出来呢，可是我改为lon1 = 207.5 lon2 = 228.5之后还是有错呢
还 ...

你去看一下coll2gr的官方文档，上面有详细用法和注意事项，我看你的资料好像不是等高面上的资料，后面-u的那个参数可能要不同的

cathy

river 发表于 2013-11-23 11:59

                               
登录/注册后可看大图

你去看一下coll2gr的官方文档，上面有详细用法和注意事项，我看你的资料好像不是等高面上的资料，后面-u的 ...

谢谢@river  ！问题解决了，只有等压面上的资料可以用-u，若等高面上的资料则需要改为所用等高面的层数，比如此处的26 。

river

cathy 发表于 2013-11-23 14:35

                               
登录/注册后可看大图

谢谢@river  ！问题解决了，只有等压面上的资料可以用-u，若等高面上的资料则需要改为所用等高面的层数， ...

解决问题就好，这也给其他人一个提示

jeffggg

为什么我做出来的风场矢量箭头很稀疏？本来分辨率挺高的资料 结果出来很不好，这是什么原因造成的呢？我完全按照上述方法做的。

陌上尘

大神啊，收藏之！！！！！！