load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
;load "./WRFUserARW.ncl"

begin
;
; The WRF ARW input file.
; This needs to have a ".nc" appended, so just do it.
  a = addfile("../full/wrfout_d01_2012-01-27_00:00:00.nc","r")

; We generate plots, but what kind do we prefer?
  type = "x11"
; type = "pdf"
; type = "ps"
; type = "ncgm"
  wks = gsn_open_wks(type,"plt_PW")
  gsn_define_colormap(wks,"gsdtol")

; Set some basic resources
  res = True
  res@MainTitle                   = "REAL-TIME WRF"
  res@Footer = False

  pltres = True
  mpres = True

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;


  pw = wrf_user_getvar(a,"pw",(/1,2,3,4/))    		; get only a couple of times
  times = wrf_user_getvar(a,"times",(/1,2,3,4/))  	; get the same time information
  ntimes = dimsizes(times)         			; get number of times extracted
  
  opts = res                          
  opts@cnFillOn = True  
  opts@lbLabelBarOn = True

  do it = 0,ntimes-1  					; loop over all times we got from file
    opts@TimeLabel = times(it)      			; Set Valid time to use on plots
    contour  = wrf_contour(a,wks,pw(it,:,:),opts)
    plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)
  end do  

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;


end