NCL画风场图

Windy1998

begin
f = addfile("总数据.nc","r")
u1=short2flt(f->u10)
v1=short2flt(f->v10)
u=dim_avg_n_Wrap(u1,0);
v=dim_avg_n_Wrap(v1,0);
;*************************************************
; calculate divergence: Use Wrap to include meta data
;*************************************************
  div = uv2dvG_Wrap(u,v)                ; u,v ==> divergence
;*************************************************
; calculate divergent wind components
;*************************************************   
  ud    = new ( dimsizes(u), typeof(u), u@_FillValue )
  vd    = new ( dimsizes(v), typeof(v), v@_FillValue )
  dv2uvg(div,ud,vd) ; div  ==> divergent  wind components
  copy_VarCoords(u, ud )
  copy_VarCoords(u, vd )
  ud@long_name  = "Zonal Divergent Wind"
  ud@units      = u@units
  vd@long_name  = "Meridional Divergent Wind"
  vd@units      = v@units
;*************************************************
; plot results
;*************************************************   
  wks  = gsn_open_wks("png","wind")           ; send graphics to PNG file
                                             
  res                 = True
  res@vcRefMagnitudeF = 3.                    ; make vectors larger
  res@vcRefLengthF    = 0.050                 ; reference vector length
  res@vcGlyphStyle    = "CurlyVector"         ; turn on curly vectors
  res@vcMinDistanceF  = 0.012                 ; thin the vectors
  res@gsnLeftString   = "Divergent Wind"
                                              ; plot 1st time step
  plot= gsn_csm_vector_map(wks,ud,vd,res)

end