;   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_1")       ; Create a plot workstation


  ter = a->HGT_M(0,:,:)                      ; Read the variable to memory
  
  plot = gsn_contour(wks,ter,True)           ; Create plot
   
end