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

begin
eos_file = addfile("AIRS_SWath_test.he2","r")
;print(eos_file)

var = eos_file->topog_L2_Standard
lon = eos_file->Longitude_L2_Standard
lat = eos_file->Latitude_L2_Standard
var@lon2d = lon
var@lat2d = lat
;var@_FillValue = -1
printVarSummary(var)


xwks = gsn_open_wks("x11","airs_2002_topog")

resources = True
resources@cnFillOn = True

;define color map
cmap = (/"(/1.00, 1.00, 1.00/)", "(/0.00, 0.00, 0.00/)", \
         "(/.560, .500, .700/)", "(/.300, .300, .700/)", \
         "(/.100, .100, .700/)", "(/.000, .100, .700/)", \
         "(/.000, .300, .700/)", "(/.000, .500, .500/)", \
         "(/.000, .700, .100/)", "(/.060, .680, .000/)", \
         "(/.550, .550, .000/)", "(/.570, .420, .000/)", \
         "(/.700, .285, .000/)", "(/.700, .180, .000/)", \
         "(/.870, .050, .000/)", "(/1.00, .000, .000/)", \
         "CadetBlue", "Ivory", "LimeGreen", "DarkSalmon"/)

;rgbcmap = hsvrgb(cmap)
gsn_define_colormap(xwks,cmap)

plot = gsn_csm_contour_map_ce(xwks,var,resources)

delete(plot)
delete(var)
delete(resources)
end
~