;---Read SLP variable off NetCDF file filename = "SLP.1979_2003.nc" a = addfile(filename,"r") slp = a->slp ; time x lat x lon (300 x 73 x 144) ;---Open an X11 window wks = gsn_open_wks("x11","slp1a") ;---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 ;---Create plots of first four time steps plot0 = gsn_csm_contour_map(wks,slp(0,:,:),res) plot1 = gsn_csm_contour_map(wks,slp(1,:,:),res) plot2 = gsn_csm_contour_map(wks,slp(2,:,:),res) plot3 = gsn_csm_contour_map(wks,slp(3,:,:),res) ;---Panel all four plots in a 2 x 2 configuration pres = True gsn_panel(wks,(/plot0,plot1,plot2,plot3/),(/2,2/),pres)