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

f = addfile("./data/sst8292a.nc", "r")
T = f->SSTA

wks = gsn_open_wks("pdf", "output")

res = True
res@cnFillOn = True
res@gsnSpreadColors = True
res@lbOrientation = "Vertical"

res@gsnDraw = False
res@gsnFrame = False

;
; These viewport values I got from the gsn_panel routine
; by setting "panelres@gsnPanelDebug = True".
;
scale = 1.03     ; Make the width/height slightly bigger than what
                 ; gsn_panel is using.

res@vpXF      = 0.032272
res@vpYF      = 0.593734
res@vpWidthF  = 0.376281*scale
res@vpHeightF = 0.188141*scale

plot1 = gsn_csm_contour_map(wks, T(1,:,:),res)

res@vpXF      = 0.53           ; originally 0.55
res@vpYF      = 0.593734
res@vpWidthF  = 0.376281*scale
res@vpHeightF = 0.188141*scale

plot2 = gsn_csm_contour_map(wks, T(2,:,:),res)

a4_height   = 29.7
a4_width    = 23
cm_per_inch = 2.54

psres                = True
psres@gsnMaximize    = True
psres@gsnPaperWidth  = a4_width/cm_per_inch
psres@gsnPaperHeight = a4_height/cm_per_inch

maximize_output(wks,psres)      ; Special procedure for maximizing plots
                                ; in a frame.
end