f = addfile ("wrfout_d01_2003-07-15_00:00:00.nc", "r") hgt = wrf_user_getvar(f,"HGT",0) wks = gsn_open_wks("png" ,"wrfgeo_wrf") gsn_define_colormap(wks,"BlAqGrYeOrReVi200") res = True res@cnFillOn = True res@ContourParameters = (/ 250., 3250., 100. /) contour = wrf_contour(f,wks,hgt,res) pltres = True mpres = True plot = wrf_map_overlays(f,wks,contour,pltres,mpres)