helicity at 1000m height
helicity at 3000m height
This shows helicity change every hour from August 17, 2018 at 12:00 to August 24, 2018 at 06:00. We referred to the NCL site’s code.
begin
idir = "/home/numodel2/SOULIK/WRFOUT/"
ifil = "wrfout_d01_2018-08-15_12:00:00"
a = addfile(idir + ifil, "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
Comments