;   Example script - read and plot data from WRF time series output
;   November 2008

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin

  type = "x11"
; type = "pdf"
; type = "ps"
; type = "ncgm"
  wks = gsn_open_wks(type,"plt_ts")       ; Create a plot workstation


  data = readAsciiTable("pt_a.d01.TS", 19, "float", 1)   ; read the data file

  res = True
  res@tmXTOn         = False        
  res@xyLineThicknesses = 2        
  res@xyLineColor    =  "blue"    
  res@tiYAxisString = "2m Temperature"
  res@tiXAxisString = "Simulation Time"
  res@tiMainString = "12Z 24 Jan 2000"

  plot  = gsn_csm_xy(wks,data(:,1),data(:,5),res)    ; we know column 2 is time and column 6 is T2
                                                     ; read README.tslist in WRF/run for more details
   
end