;************************************
; CSM_Graphics: panel_6.ncl
;************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;************************************
begin

 in = addfile("h_avg_Y0191_D000.00.nc","r")
 T=in->T(:,:,:,{0:360})               ; get rid of cyclic points
 S=in->S(:,:,:,{0:360})  
;************************************
;creat plot
;************************************
 wks = gsn_open_wks("ps","panel")              ; open a ps  file
 gsn_define_colormap(wks,"gui_default")        ; choose colormap
 plot = new(4,graphic)                         ; create graphic array

 res                      = True
 res@cnFillOn             = True         ; turn on color fill
 res@gsnDraw              = False        ; do not draw picture
 res@gsnFrame             = False        ; do not advance frame
 res@gsnPolar             = "NH"         ; select northern hemisphere
 res@lbOrientation        = "Vertical"   ; vertical label bar

 plot(0)=gsn_csm_contour_map_polar(wks,T(0,0,:,:),res)
 plot(1)=gsn_csm_contour_map_polar(wks,T(0,5,:,:),res)
 plot(2)=gsn_csm_contour_map_polar(wks,S(0,0,:,:),res)
 plot(3)=gsn_csm_contour_map_polar(wks,S(0,3,:,:),res)

; draw panel without white space
 gsn_panel(wks,plot,(/2,2/),False)

; draw panel with white space added
 resP                 = True
 resP@gsnPanelYWhiteSpacePercent = 5
 resP@gsnPanelXWhiteSpacePercent = 5
 gsn_panel(wks,plot,(/2,2/),resP)

 
end