pro find_single_node,time_ofs,t_lo,t_hi,a_node,a_node_count,d_node,d_node_count
  t_node=rtbis(t_lo,t_hi,0,"find_node_f",0.01,state=state,time_ofs=time_ofs)
  ;Is it an ascending or descending node?
  if find_node_f(t_node+1,time_ofs=time_ofs) gt 0 then begin
    ;It's an ascending node
    if a_node_count ge n_elements(a_node) then a_node=[a_node,dblarr(16)]
    a_node[a_node_count]=t_node+time_ofs
    a_node_count++
  end else begin
    ;It's a descending node
    if d_node_count ge n_elements(d_node) then d_node=[d_node,dblarr(16)]
    d_node[d_node_count]=t_node+time_ofs
    d_node_count++
  end
  ;Set up for the next node search
  t_lo=t_node+1d6
  t_hi=t_lo+1d6
  f_lo=find_node_f(t_lo,time_ofs=time_ofs)
  f_hi=find_node_f(t_hi,time_ofs=time_ofs)
  while(f_hi*f_lo gt 0) do begin
    t_hi=t_lo+(t_hi-t_lo)*2
    f_hi=find_node_f(t_hi,time_ofs=time_ofs)
  end
end