;************************************************
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 "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
;************************************************

begin
  a = addfile("try.nc","r")
  hgt = wrf_user_getvar(a,"HGT",0)
  hgt@lat2d = wrf_user_getvar(a,"XLAT",0)   ; latitude/longitude
  hgt@lon2d = wrf_user_getvar(a,"XLONG",0)  ; required for plotting

  wks = gsn_open_wks("png","hgt_gsn")

  res               = True

  res@gsnMaximize   = True   ; maximize plot in frame

  res@cnFillOn      = True  
  res@cnFillPalette = "OceanLakeLandSnow"
  res@cnLinesOn     = False

  res@mpProjection  = "CylindricalEquidistant"    ; The default
  res@mpDataBaseVersion = "MediumRes"

  res@gsnAddCyclic      = False

;---Zoom in on plot
  res@mpMinLatF     = min(hgt@lat2d)
  res@mpMaxLatF     = max(hgt@lat2d)
  res@mpMinLonF     = min(hgt@lon2d)
  res@mpMaxLonF     = max(hgt@lon2d)

  contour = gsn_csm_contour_map(wks,hgt,res)

end