; Generate a series of orbit_info structures
; input
;   t_start - time of beginning of first orbit to generate. Assumed to be at an ascending node
;   n_start - Orbit number of first orbit to generate
;   t_stop - generate orbits up to the one which overlaps this time
;   (all times are in GPS microseconds)
; return 
;   An array of orbit_info structures. Note that the preliminary flag is not set either way,
;   so it takes the default value, definitive.
;
; Commonly, if you already have orbits in the orbit table and you want more, call it like this:
;
;   last_orbit=get_last_orbit_info()
;   more_orbits=generate_orbit_info(last_orbit.stop_time,last_orbit.orbit_number+1,t_stop)

function generate_next_orbit_info,last_orbit_info
  n_orbit=last_orbit_info.orbit_number+1
  orbit_info=find_next_orbit(last_orbit_info.stop_time)
  orbit_info.orbit_number=n_orbit
  result=orbit_info
  print,"New orbit #"+string(orbit_info.orbit_number,format='(%"%d")')+", starting at "+usec2vms(orbit_info.start_time)
  return,result
end