; Example script - plot terrain (using some basic WRF scripts) ; Overwrite basic map background settings ; November 2008 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl" begin a = addfile("./geo_em.d01.nc","r") ; Open a file type = "x11" ; type = "pdf" ; type = "ps" ; type = "ncgm" wks = gsn_open_wks(type,"plt_geo_6") ; Create a plot workstation opts = True ; Set some Basic Plot options opts@MainTitle = "GEOGRID FIELDS" opts@InitTime = False ; Do not plot time or footers opts@Footer = False ter = wrf_user_getvar(a,"HGT_M",0) ; Read the variable to memory res = opts ; Use basic options for this field res@cnFillOn = True ; Create a color fill plot res@ContourParameters = (/ 0., 1000., 50. /) ; Set the levels contour = wrf_contour(a,wks,ter,res) pltres = True ; Set plot options mpres = True ; Set map options mpres@mpGeophysicalLineColor = "Black" ; Overwrite basic map settings mpres@mpGridLineColor = "Black" mpres@mpLimbLineColor = "Black" mpres@mpNationalLineColor = "Black" mpres@mpPerimLineColor = "Black" mpres@mpUSStateLineColor = "Black" plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres) ; Plot field over map background end