想看看台风，用FNL资料画风场图。

morrocco

最近北美台风新闻很多
所以试着画画
ncl官方教程没找到画fnl的资料，所以用的era的模板
想把第一层的U V计算个合成风速大小画出来
未果，求指点

;************************************************
; era40_1.ncl
;
; Concepts illustrated:
;   - Reading GRIB data
;   - Plotting ERA40 data
;   - Spanning the full color map for contour fill
;   - Drawing color-filled contours over a cylindrical equidistant map
;   - Selecting a different color map
;
;*************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"  
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"  
;************************************************
begin
;************************************************
; create pointer to file and read in data
;************************************************
  in = addfile("fnl_20170906_00_00.grib2","r")
  U = in->UGRD_P0_L100_GLL0(0,:,:)
  V = in->VGRD_P0_L100_GLL0(0,:,:)
;************************************************
; create default plot
;************************************************
  spd=sqrt(U*U+V*V)

  wks = gsn_open_wks("png","fnl")               ; open a ps file
  gsn_define_colormap(wks,"BlAqGrYeOrRe")        ; choose colormap

  res                       = True               ; plot mods desired
  res@cnFillOn              = True               ; turn on color fill
  res@cnLinesOn             = False              ; turn off contour lines

;---This resource not needed in V6.1.0
  res@gsnSpreadColors       = True               ; use full range of color map

  plot = gsn_csm_contour_map_ce(wks,spd(:,:),res)     
end