;   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_3")       ; 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
  res@cnLevelSelectionMode = "ManualLevels"  ; Set levels manually
  res@cnMinLevelValF = 0.
  res@cnMaxLevelValF = 1000.
  res@cnLevelSpacingF = 50.
  
  plot = gsn_contour(wks,ter,res)            ; Create plot
   
end