;*************************************************
; gsn_xy_2.ncl
;*************************************************
load          "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
;*************************************************
begin
  data  = asciiread("./data.asc",(/500,6/),"float")
  x     = data(:,1)
  y     = data(:,4)
;**************************************************
; create plot
;**************************************************
  wks = gsn_open_wks("ps","gsn_xy")        ; Open a workstation.

  res               = True                   ; plot mods desired
  res@tiMainString  = "An xy plot Example"   ; title
  res@tiYAxisString = "Dust (ppm)"           ; y axis title
  res@tiXAxisString = "Time"                 ; x axis title
  plot = gsn_xy(wks,x,y,res)               ; Draw an XY plot with 1 curve.
  
end