anchor.diff:  c     bpow = value of {p(i)} at the highest peak
   c     depth= depth of the transit at   *bper*
   c     qtran= fractional transit length  [ T_transit/bper ]
  -c     in1  = bin index at the start of the transit [ 0 < in1 < nb+1 ]
  -c     in2  = bin index at the end   of the transit [ 0 < in2 < nb+1 ]
  +c     in1  = bin index at the start of the transit [ 0 < in1 < nbins+1 ]
  +c     in2  = bin index at the end   of the transit [ 0 < in2 < nbins+1 ]
   c
   c
   c     Remarks:
   c     ~~~~~~~~
   c
   c     -- *fmin* MUST be greater than  *1/total time span*
  -c     -- *nb*   MUST be lower than  *nbmax*
  +c     -- *nbins*   MUST be lower than  *nbinsmax*
   c     -- Dimensions of arrays {y(i)} and {ibi(i)} MUST be greater than
  -c     or equal to  *nbmax*.
  +c     or equal to  *nbinsmax*.
   c     -- The lowest number of points allowed in a single bin is equal
  -c     to   MAX(minbin,qmi*N),  where   *qmi*  is the minimum transit
  +c     to   MAX(minbin,qmin*N),  where   *qmin*  is the minimum transit
   c     length/trial period,   *N*  is the total number of data points,
   c     *minbin*  is the preset minimum number of the data points per
   c     bin.
  
positive.code:             ph     = u(i)*f0
              ph     = ph-idint(ph)
  c
              j      = 1 + idint(nb*ph)
              ibi(j) = ibi(j) + 1
              y(j) =   y(j) + v(i)
   102     continue
  c
  c-----------------------------------------------
  c     Extend the arrays  ibi()  and  y() beyond
  c     nb   by  wrapping
  c
           do 104 j=nb1,nbkma
              jnb    = j-nb
              ibi(j) = ibi(jnb)
              y(j) =   y(jnb)
   104     continue
  
negative.code: c
  c
        subroutine eebls(n,t,x,u,v,nf,fmin,df,nbins,qmin,qmax,
       &     p,bper,bpow,depth,qtran,in1,in2)
  c
  c------------------------------------------------------------------------
  