function ratime, value
hour = long(abs(value)/15.0)
min = long((abs(value)/15.0 - hour)*60.0 + 0.2)
return, string(hour, min, format="(i2.2,'h',i2.2,'m')")
end

function decdeg, value
deg = long(value)
min = long(abs((value - deg)*60.0)+0.1)
return, string(deg, min, format="(i3.2,'!9%!X',i2.2)")+string("222B)
end


pro plotso

set_plot, 'z'

loadct, 39

!P.background = 255

SPAWN, 'wc names', wc_string
num_points = long(STRMID(wc_string(0),0,8))
num_cols   = long(STRMID(wc_string(0),8,8))/num_points
data = STRARR(num_cols,num_points)
OPENR, unit, 'names', /GET_LUN
READF, unit, data
CLOSE, unit
FREE_LUN, unit

cor = read_data('catalog')

box00 = read_data('obox211036039.07.0')
box01 = read_data('obox211036039.07.1')
box02 = read_data('obox211213704.09.1')
box03 = read_data('obox211213704.10.0')
box04 = read_data('obox211254615.01.1')
box05 = read_data('obox211254615.02.0')
box06 = read_data('obox211271063.08.0')
box07 = read_data('obox211271063.08.1')
box08 = read_data('obox211277063.06.0')
box09 = read_data('obox211277063.06.1')




latvals = 4.53333333333333 + findgen(11) * 0.216666666666667
lonvals = -281.004681545625 + findgen(11) * 0.25
 
latlbs = strarr(11);
lonlbs = strarr(11);
for i=0,10 do begin
latlbs[i] = decdeg(latvals[i])
if (lonvals[i] le 0.0) then $
lonlbs[i] = ratime(lonvals[i]) $
else lonlbs[i] = ratime(lonvals[i]-360.0)

end

map_set, /gnomic, 5.53333333333333, -280, /isotropic, $
	limit=[4.53333333333333, -281.004681545625, 6.53333333333333, -278.995318454375], $
	title='Sky Map Centered at ASM Error Box Intersection', color=0

map_grid, label=1, latalign=0.5, latlab=-279.120318454375, lonalign=0.5, $
        lonlab=4.64166666666667, lats=latvals, lons=lonvals, charsize=0.7, $
	lonnames=lonlbs, latnames=latlbs, color=0

for i=0,num_points-1 do begin
	xyouts, -cor(0,i)+0.5, cor(1,i), data(0,i), color=0
endfor

oplot, -cor(0,*), cor(1,*), psym=4, symsize=0.6, color=0
oplot, -box00(0,*), box00(1,*), color=220
oplot, -box01(0,*), box01(1,*), color=220
oplot, -box02(0,*), box02(1,*), color=220
oplot, -box03(0,*), box03(1,*), color=220
oplot, -box04(0,*), box04(1,*), color=220
oplot, -box05(0,*), box05(1,*), color=220
oplot, -box06(0,*), box06(1,*), color=220
oplot, -box07(0,*), box07(1,*), color=220
oplot, -box08(0,*), box08(1,*), color=220
oplot, -box09(0,*), box09(1,*), color=220




filename =  'map.gif'
write_gif, filename, tvrd()

end