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('obox0.-2147483648.0')
box01 = read_data('obox287731536.17.1')

bat = read_data('bat.dat')

latvals = 2.26666666666667 + findgen(11) * 0.216666666666667
lonvals = -286.751627507423 + 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, 3.26666666666667, -285.75, /isotropic, $
	limit=[2.26666666666667, -286.751627507423, 4.26666666666667, -284.748372492577], $
	title='Sky Map Centered at ASM Error Box Intersection', color=0

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

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

end