f = addfile ("wrfout_d01_2003-07-15_00:00:00.nc", "r") tc = wrf_user_getvar(f,"tc",0) ; temperature in degC printVarSummary(tc) wks = gsn_open_wks("png" ,"panel1a") ; Open png file "panel1a.png" res = True res@gsnDraw = False ; Don't draw plots res@gsnFrame = False ; Don't advance frame res = wrf_map_resources(f,res) ; Add necessary resources for WRF map res@tfDoNDCOverlay = True ; Use the native WRF map projection res@gsnAddCyclic = False ; Don't add longitude cyclic point plot1 = gsn_csm_contour_map(wks,tc(0,:,:),res) ; Create contour plots plot2 = gsn_csm_contour_map(wks,tc(3,:,:),res) ; at different levels. plot3 = gsn_csm_contour_map(wks,tc(5,:,:),res) plot4 = gsn_csm_contour_map(wks,tc(7,:,:),res) ;---Create a panel of plots with 2 rows and 2 columns. gsn_panel(wks,(/plot1,plot2,plot3,plot4/),(/2,2/),False)