pro burster_position, bp, qstop = qstop
rad = 180./!pi
bp = {burst, day : 1, time : 1l, sep : 1., solar_n : 1.}
bp = replicate(bp,60)
b_ra = 17. + 44./60 + 28./3600
b_ra = b_ra*15.
b_dec = -28. + 45./60
day = indgen(60) + 1
month = intarr(60) + 12
hour = intarr(60)
year = intarr(60)+95
p_angle = fltarr(60)
for i = 0, 59 do begin
  if i gt 31 then begin
    day(i) = day(i) - 31
    year(i) = year(i) + 1
    month(i) = month(i) -11
  endif
endfor
ra = fltarr(60)
dec = fltarr(60)
for i = 0,59 do begin
  sun,year(i),month(i),day(i),hour(i),true_ra = true_ra, true_dec = true_dec,$
    pa = pa
  p_angle(i) = pa
  ra(i) = true_ra 
  dec(i) = true_dec 
endfor
ra = ra * 15.
ang = fltarr(60)
dlong = ra - b_ra
dlat = dec - b_dec
b_angle = atan(dlat,dlong)
ang = sphere_dist(dec, ra, b_dec, b_ra)
bp.time = lonarr(60)
bp.day = 6178 + day
for i = 0,59 do if i gt 31 then bp(i).day = bp(i).day + 31
bp.sep = ang
bp.solar_n = dlat * cos((p_angle + b_angle) / rad)
if n_elements(qstop) ne 0 then stop
end