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('obox133973414.20.0')
box01 = read_data('obox133973414.20.1')
box02 = read_data('obox133973414.21.1')
box03 = read_data('obox133985718.10.1')
box04 = read_data('obox133985718.11.0')
box05 = read_data('obox133985718.12.0')
box06 = read_data('obox147849927.02.0')
box07 = read_data('obox147849927.02.1')
box08 = read_data('obox182093783.22.0')
box09 = read_data('obox182093783.22.1')
box10 = read_data('obox211225576.10.1')
box11 = read_data('obox211225576.12.0')
box12 = read_data('obox211266743.09.0')
box13 = read_data('obox211266743.10.1')
box14 = read_data('obox216595223.07.0')
box15 = read_data('obox216595223.07.1')
box16 = read_data('obox247354464.12.0')
box17 = read_data('obox247354464.12.1')

bat = read_data('bat.dat')

latvals = 33.7 + findgen(11) * 0.216666666666667
lonvals = -300.966331871482 + 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, 34.7, -299.75, /isotropic, $
	limit=[33.7, -300.966331871482, 35.7, -298.533668128518], $
	title='Sky Map Centered at ASM Error Box Intersection', color=0

map_grid, label=1, latalign=0.5, latlab=-298.658668128518, lonalign=0.5, $
        lonlab=33.8083333333333, 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, -bat(0,*), bat(1,*), color=40

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

end