; Not quite the original
; My code is something like this:

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin
   file1   = addfile("/project/cas/shea/t_jja_mean.nc","r")
   tera_in = short2flt(file1->t(0,12:22,:,:))
  ;tera_in = tera_in(:,::-1,:)                      ; reverse lat order
   printVarSummary( tera_in )
  ;tera    = tera_in(longitude|:,latitude|:,levelist|:)
   tera    = tera_in(latitude|:,longitude|:,levelist|:)
   printVarSummary( tera )

   wks = gsn_open_wks("x11" ,"plot1")
   gsn_define_colormap(wks,"testcmap")
   res = True
   res@cnFillOn= True
   res@gsnSpreadColors= True
   res@mpCenterLonF   = 180
   printVarSummary(tera(:,:,0))
   
   printVarSummary(tera(:,:,0))                     ; dimension reduction [ 2D ]
   plot = gsn_csm_contour_map_ce(wks, tera(:,:,0), res )
end