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" ,"panel1d") ; Open png file "panel1d.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@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 = 2 ; ; Loop across four levels and create a contour plot for each one. ; Change the main title each time. ; plots = new(4,graphic) do i=0,3 res@tiMainString = "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 2 columns. pres = True ; Set panel resources. pres@gsnPanelLabelBar = True ; Turn on panel labelbar. gsn_panel(wks,plots,(/2,2/),pres) ; Create a panel of plots with ; 2 rows and 2 columns.