;   Example script - plot terrain (using only basic NCL functions)
;   Note no map background info is available
;   November 2008

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

begin

  a = addfile("./geo_em.d01.nc","r")   ; Open a file


  type = "x11"
; type = "pdf"
; type = "ps"
; type = "ncgm"
  wks = gsn_open_wks(type,"plt_geo_2")     ; Create a plot workstation


  ter = a->HGT_M(0,:,:)                    ; Read the variable to memory

  res = True                               ; Create some plot resources
  res@cnFillOn = True                      ; Create a color fill plot
  res@gsnSpreadColors = True
  
  plot = gsn_contour(wks,ter,res)          ; Create plot
   
end