;---Read SLP variable off NetCDF file filename = "SLP.1979_2003.nc" a = addfile(filename,"r") slp = a->slp ; 300 x 73 x 144 slp@units = "mb" ; "millibars" ;---Open an X11 window wks = gsn_open_wks("x11","slp1f") ;---Set some plotting options res = True res@gsnDraw = False res@gsnFrame = False res@gsnLeftString = "" res@gsnRightString = "" res@cnFillOn = True ; Turn on color fill res@cnLinesOn = False ; Turn off contour lines res@cnLineLabelsOn = False ; Turn off contour line labels res@lbLabelBarOn = False ; Turn off labelbar res@cnLevelSelectionMode = "ExplicitLevels" res@cnLevels = ispan(976,1060,4) ;---Loop across nplots time steps and create a plot nplots = 7 plot = new(nplots,graphic) do i=0,nplots-1 plot(i) = gsn_csm_contour_map(wks,slp(i,:,:),res) end do ;---Panel all seven plots in a 2 x 3 x 2 configuration pres = True pres@gsnPanelLabelBar = True pres@txString = slp@long_name + " (" + slp@units + ")" pres@gsnPanelRowSpec = True gsn_panel(wks,plot,(/2,3,2/),pres)