;
; These four load statements are optional in NCL V6.4.0.

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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "./gsn_csm_pres_hgt_mod.ncl"

begin
 f      = addfile ("$NCARGTEST/nclscripts/cdf_files/uvt.nc" , "r")
 lat    = f->lat                                    ; get latitude
 ind14S = ind(lat.eq.-1.395307)                     ; subscript index of 1.4S
 u      = f->U(0,:,ind14S,:)                        ; get u at 1.4S

 wks   = gsn_open_wks("x11","pres_hgt")

;---Default plot
 res                 = True                        ; plot mods desired
 res@gsnMaximize     = True
 res@tiMainString    = "January 1988"              ; title
 res@cnLevelSpacingF = 5.0                         ; contour level spacing

 plot  = gsn_csm_pres_hgt(wks, u, res ) 

;---Change the right Y axis
 res@gsnPresHgtHeightLabel = "My own height label"
 res@tiMainString          = "Change the height label on the right axis"
 plot  = gsn_csm_pres_hgt(wks, u, res ) 

;---Change the left Y axis too
 res@tiYAxisString         = "My own pressure label"
 res@tiMainString          = "Change the pressure label on the left axis too"
 plot  = gsn_csm_pres_hgt(wks, u, res ) 

end