;---Read SLP variable off NetCDF file filename = "SLP.1979_2003.nc" a = addfile(filename,"r") slp = a->slp ; 300 x 73 x 144 ;---Open an X11 window wks = gsn_open_wks("x11","slp1b") ;---Set some plotting options res = True res@cnFillOn = True ; Turn on color fill res@cnLinesOn = False ; Turn off contour lines res@cnLineLabelsOn = False ; Turn off contour line labels ;---Loop across first four time steps and create a plot nplots = 4 plot = new(nplots,graphic) do i=0,nplots-1 plot(i) = gsn_csm_contour_map(wks,slp(i,:,:),res) end do ;---Panel all four plots in a 2 x 2 configuration pres = True gsn_panel(wks,plot,(/2,2/),pres)