load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/cd_string.ncl" ;---Read SLP variable off NetCDF file filename = "SLP.1979_2003.nc" a = addfile(filename,"r") slp = a->slp ; 300 x 73 x 144 newtime = cd_string(slp&time,"%c %D %Y") ; i.e. Mar 01 2002 ;---Open an X11 window wks = gsn_open_wks("x11","slp1b") ;---Set some plotting options res = True res@gsnDraw = False res@gsnFrame = False res@gsnLeftString = "" ; turn off long_name string on left res@gsnRightString = "" ; turn off units string on right res@cnLevelSelectionMode = "ExplicitLevels" res@cnLevels = ispan(976,1060,4) res@lbLabelBarOn = False 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 = 16 ; 16 plots 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 4 x 4 configuration pres = True pres@gsnPanelLabelBar = True ; turn on common labelbar pres@txString = slp@long_name + " (" + slp@units + ")" ; turn on a main title pres@gsnPanelFigureStrings = newtime(0:nplots-1) gsn_panel(wks,plot,(/4,4/),pres)