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('obox139593926.21.0')
box01 = read_data('obox139593926.21.1')
box02 = read_data('obox139593926.22.0')
box03 = read_data('obox139593926.22.1')
box04 = read_data('obox139599974.25.1')
box05 = read_data('obox139599974.28.0')
box06 = read_data('obox139669526.00.2')
box07 = read_data('obox139669526.01.2')
box08 = read_data('obox139669526.02.2')
box09 = read_data('obox139721477.18.1')
box10 = read_data('obox139721477.20.0')
box11 = read_data('obox139796119.02.2')
box12 = read_data('obox139796119.03.2')
box13 = read_data('obox139796119.04.2')
box14 = read_data('obox139796119.06.2')
box15 = read_data('obox139882279.00.2')
box16 = read_data('obox139882279.03.2')
box17 = read_data('obox139882279.04.2')
box18 = read_data('obox139951399.04.2')
box19 = read_data('obox139951399.05.2')
box20 = read_data('obox291847296.03.0')
box21 = read_data('obox291847296.03.1')
box22 = read_data('obox292454608.09.0')
box23 = read_data('obox292454608.10.1')

bat = read_data('bat.dat')

latvals = -34.2 + findgen(11) * 0.216666666666667
lonvals = -268.44507964664 + 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, -33.2, -267.25, /isotropic, $
	limit=[-34.2, -268.44507964664, -32.2, -266.05492035336], $
	title='Sky Map Centered at ASM Error Box Intersection', color=0

map_grid, label=1, latalign=0.5, latlab=-266.17992035336, lonalign=0.5, $
        lonlab=-32.3083333333333, 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
oplot, -box10(0,*), box10(1,*), color=220
oplot, -box11(0,*), box11(1,*), color=220
oplot, -box12(0,*), box12(1,*), color=220
oplot, -box13(0,*), box13(1,*), color=220
oplot, -box14(0,*), box14(1,*), color=220
oplot, -box15(0,*), box15(1,*), color=220
oplot, -box16(0,*), box16(1,*), color=220
oplot, -box17(0,*), box17(1,*), color=220
oplot, -box18(0,*), box18(1,*), color=220
oplot, -box19(0,*), box19(1,*), color=220
oplot, -box20(0,*), box20(1,*), color=220
oplot, -box21(0,*), box21(1,*), color=220
oplot, -box22(0,*), box22(1,*), color=220
oplot, -box23(0,*), box23(1,*), color=220

oplot, -bat(0,*), bat(1,*), color=40

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

end