; *************************************************************************
; veceff_1.ncl
;
; Concepts illustrated:
;   - Drawing vector plots
;   - Changing the length of the vectors
;   - Changing the length of the smallest vector as a fraction of the reference vector
;   - Coloring vectors based on magnitude
;   - Spanning the full color map for colored vectors
;   - Creating a main title
;   - Transposing an array
;   - Naming dimensions of an array
;   - Generating dummy data
;
; *************************************************************************
;
; These files are loaded by default in NCL V6.2.0 and newer
; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin
;
; Generate dummy vector data arrays
;
  M  = 30
  N  = 25
  PI = 3.14159

  V = 10.0 * cos(onedtond((2.0 * PI / M) * ispan(0,M-1,1),(/N,M/)))

;---Use transpose function to exchange dimensions of U to avoid manual dimension reordering
  U  = transpose(10.0 * cos(onedtond((2.0 * PI / N) * ispan(0,N-1,1),(/M,N/))))
 
  wks   = gsn_open_wks ("png", "veceff")       ; send graphics to PNG file

  res                      = True

  res@gsnMaximize          = True    ; Maximize plot in frame

  res@tiMainString         = "No resources set"

  vector = gsn_csm_vector(wks,U,V,res)

  res@vcMonoLineArrowColor = False  ; color arrows based on magnitude
  res@vcRefLengthF         = 0.03313608
  res@vcMinFracLengthF     = 0.3
  res@vcLevelPalette       = "temp1" ; set color map

  res@tiMainString         = "Line-Drawn Vectors (colored by magnitude)"

  vector = gsn_csm_vector(wks,U,V,res)

  res@vcLineArrowThicknessF= 3. 
  res@tiMainString         = "Thickened Vectors"

  vector = gsn_csm_vector(wks,U,V,res)
end