; 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