pro aod500

file=[‘./data3.prn’]
aot_av=fltarr(3,81)
a=fltarr(3)
dum=’ ‘

close,1
openr,1,file
readf,1,dum

for i =0,80 do begin
readf,1,a,format='(f9.3,f8.3,f9.6)’
aot_av(*,i)=a(*)
endfor
close,1
print, aot_av

lon=aot_av[0,*]
lat=aot_av[1,*]
aot=aot_av[2,*]

levels=indgen(70)*0.1
device, decomposed=0
loadct, 39
!p.background=255
!p.color=0

window, 0, xsize=720, ysize=500
map_set, /CYLINDRICAL, POSITION=[0.03, 0.05, 0.88, 0.95], title=’AOT 500′
contour, aot, lon, lat, /fill, /irregular, levels=levels, /overplot
plots, lon, lat, psym=4, symsize=2, thick=2, color=255
colorbar,lev, col, /vertical, /right, range=[0,6], $
position=[0.92, 0.15, 0.95, 0.86], Xminor=0, ncolors=255

map_continents, /coast, /countries
map_grid, box_axes=0.1, londel=15, latdel=10, charsize=0.65

end

contour와 plots를 이용하여 지도 위에 점을 표시, 그 주변 일정 반경안에 등고선처럼 같은 값일때 같은 색이 나타나도록 표현하고 싶습니다.
이것저것 시도해 봤는데 contour로는 한계가 있는것 같네요.
혹시 다른 프러시저를 써야 하나요??

http://xbigfile.naver.com/bigfileupload/download?fid=xqnm1zGOD6dwHqujaxbdaxUqaIYZKAuqKxUmaxumKAKmFIYwKq2lKxvqKxvlHqEXM4tlKoUdpxMZKA0SKqvZFqv9FAvXp6JvpApSFxvZ
위의 링크는 데이터 파일입니다.