f = addfile ("wrfout_d01_2003-07-15_00:00:00.nc", "r") tc = wrf_user_getvar(f,"tc",-1) ; temperature in degC tc@units = "degC" ; better units wks = gsn_open_wks("png" ,"panel1i") ; Open png file "panel1i.png" res = True ; Plot options desired res@gsnDraw = False ; Don't draw plots res@gsnFrame = False ; Don't advance frames 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 res@cnFillOn = True ; Turn on color res@cnLineLabelsOn = False ; Turn off line labels res@cnLinesOn = False ; Turn off contour lines res@cnInfoLabelOn = False ; Turn off info label res@lbLabelBarOn = False ; Turn off individual labelbars res@cnLevelSelectionMode = "ManualLevels" ; Set contour levels res@cnMinLevelValF = 8 res@cnMaxLevelValF = 38 res@cnLevelSpacingF = 1 res@gsnLeftString = "" ; Turn off left subtitle res@gsnRightString = "" ; Turn off right subtitle ;---Loop across six levels and create a contour plot for each one. plots = new(6,graphic) do i=0,dimsizes(plots)-1 res@gsnCenterString = "bottom_top index " + i plots(i) = gsn_csm_contour_map(wks,tc(0,i,:,:),res) end do ;---Create a panel of plots with 2 rows and 3 columns. pres = True ; Set panel resources. pres@gsnPanelRowSpec = True ; Specify plots per row pres@gsnMaximize = True ; Maximize plots in panel. pres@gsnPanelLabelBar = True ; Turn on panel labelbar. pres@pmLabelBarWidthF = 0.8 ; Change labelbar width. pres@gsnPanelMainString = tc@description + " (" + tc@units + ")" gsn_panel(wks,plots,(/2,1,3/),pres)