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