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

begin

  a = addfile("../wrfout_d01_2000-01-24_12:00:00.nc","r")

  type = "x11"
; type = "pdf"
; type = "ps"
; type="png"

  wks = gsn_open_wks(type,"plt_helicity")
  gsn_define_colormap(wks, "BlAqGrYeOrRe")

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

  pltres = True
  mpres = True


  times  = wrf_user_list_times(a)   ; get times in the file
  ntimes = dimsizes(times)          ; number of times in the file

  do it =0,ntimes-1

    sreh = wrf_user_getvar(a,"helicity",it)    			; here a default of 3km is used
    ;sreh = wrf_user_getvar(a,(/"helicity","3000"/),it)    	; here 3km is specifically set - same as above
    ;sreh = wrf_user_getvar(a,(/"helicity","1000"/),it)    	; here a height of 1km is used
  
    cnres                      = res
    cnres@cnFillOn             = True
    cnres@cnSmoothingOn        = True
    cnres@cnSmoothingDistanceF = .005
  ; cnres@ContourParameters    = (/ -20., 90., 5./)
  ; cnres@gsnSpreadColorEnd    = -3  ; End 3rd from the last color in color map

    contour = wrf_contour(a,wks,sreh,cnres)
    plot   = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)

  end do

end