; XY line plot for RegCM3 output
;
; MA Snyder and NS Diffenbaugh - May 15,2004

; pathnames must be changed based on local installation of NCL
load "/usr/local/ncarg/scripts/gsn_code.ncl"
load "/usr/local/ncarg/scripts/gsn_csm.ncl"
load "/usr/local/ncarg/scripts/contributed.ncl"
load "/usr/local/ncarg/scripts/shea_util.ncl"

begin

; Set path to input file
diri = "/vol1/msnyder/RegCM3/oct6/RegCM/Cal_Ncep_opt/"
fili = "ATM1960AVG.nc"

; Set path to HEAD_OUT.nc
diri1 = "/vol1/msnyder/RegCM3/oct6/RegCM/Cal_Ncep_opt/"
fili1 = "HEAD_OUT.nc"

; Pointers to input files
f  = addfile(diri+fili,"r")
f1 = addfile(diri1+fili1,"r")


; Read in variables
var    = f->TK(:,0,:,:)      ; (time,lev,lat,lon)
var1   = f->time

; Calculate area average temperature over the entire domain
varave = wgt_areaave(var,1.0,1.0,0)

 wks  = gsn_open_wks ("ps", "xy")            ; open workstation
 
  res                        = True               ; plot mods desired
 
  res@xyLineColors           = (/"blue"/)
  res@xyLineThicknesses      = (/"4.0"/)

  res@tiMainString           = "RegCM3: xy_plot"

plot = gsn_csm_xy(wks,var1,varave,res)

end