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

djcj520

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

djcj520

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

djcj520

'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

上个图不是很好吗

hxyj

楼主辛苦了。

nzy718d

果断收藏，谢谢

狙击枪

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

风卷云

FNL的细网格在哪里下的？

djcj520

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

djcj520

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