Fnl细网格Grib2资料用Grads2.1新版沿着任意某条线经纬度画剖面图

djcj520 · 发表于 2015-8-22 17:59:43

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

借鉴前人的研究本人研究了下用fnl细网格grib2资料，并用最新Grads2.1沿着任意经纬度画图，实现得到清晰的剖面图。

djcj520 · 发表于 2015-8-22 18:01:49

本帖最后由 djcj520 于 2015-8-22 18:05 编辑

不好意思，东西还没有贴上就发表了，这貌似不能加图和附件

djcj520 · 发表于 2015-8-22 18:03:06

本帖最后由 djcj520 于 2015-8-22 18:05 编辑

'reinit'
'open G:\fnl_20150601_00_00.ctl'
'set grads off'
'set grid off'
'set zlog on'
'set csmooth on'
'run define_colors.gs'
'set font 0'
'set xlopts 1 4 0.20'
'set ylopts 1 4 0.20'
'set x 1'
'set y 1'
'set t 5'
'set lev 1000 100'
lon1 = 100.0
lon2 = 120.0
lat1 = 40.0
lat2 = 20.0
lon = lon1
'collect 1 free'
while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(VVELprs*100,'lon','lat')'
  lon = lon + 1
endwhile

'set map 1 1 9'
'set x 10 20'
'set xaxis 'lon1' 'lon2
'set clab on'

'set gxout contour'
   'set cthick 5'
   'set clopts -1 4 0.15'
   'set cmin 0'
   'set cstyle 1'
   'set ccolor rainbow'
   'set cint 20'

   'd coll2gr(1,-u)'
   'set cmax 0'
   'set cthick 5'
   'set clopts -1 4 0.15'
   'set cstyle 3'
   'set ccolor rainbow'
   'set cint 100'
   'd coll2gr(1,-u)'

*'set gxout shaded'
*'set cmax 0'
*'d coll2gr(1,-u)'
*'set gxout contour'
*'set cthick 6'
*'d coll2gr(1,-u)'
'draw title Omega along (20N,40E) to (40N,120E)'
'gxprint g:\test.png white'
;

river · 发表于 2015-8-22 19:58:36

上个图不是很好吗

hxyj · 发表于 2015-8-22 21:22:31

楼主辛苦了。

nzy718d · 发表于 2015-8-23 08:22:23

果断收藏，谢谢

狙击枪 · 发表于 2015-8-23 08:41:46

作者是否可以详细介绍如何处理grib2，生成ctl等等？

风卷云 · 发表于 2015-8-23 08:48:43

FNL的细网格在哪里下的？

djcj520 · 发表于 2015-8-24 07:36:12

http://rda.ucar.edu/datasets/ds083.2/index.html#sfol-wl-/dsszone/ds083.2?g=1

djcj520 · 发表于 2015-8-24 07:37:04

自己注册一个账号，用163的应该可以

		自动登录	找回密码
密码			立即注册

[分享资料] Fnl细网格Grib2资料用Grads2.1新版沿着任意某条线经纬度画剖面图

登录后查看更多精彩内容~

相关帖子