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