anchor.diff:  c
   c
         subroutine eebls(n,t,x,u,v,nf,fmin,df,nb,qmi,qma,
  -     &p,bper,bpow,depth,qtran,in1,in2)
  +     &     p,bper,bpow,depth,qtran,in1,in2)
   c
   c------------------------------------------------------------------------
   c     >>>>>>>>>>>> This routine computes BLS spectrum <<<<<<<<<<<<<<
   c
  -c         [ see Kovacs, Zucker & Mazeh 2002, A&A, Vol. 391, 369 ]
  +c     [ see Kovacs, Zucker & Mazeh 2002, A&A, Vol. 391, 369 ]
   c
   c     This is the slightly modified version of the original BLS routine
   c     by considering Edge Effect (EE) as suggested by
  
positive.code:       s=0.0d0
        t1=t(1)
        do 103 i=1,n
        u(i)=t(i)-t1
        s=s+x(i)
   103  continue
        s=s/rn
        do 109 i=1,n
        v(i)=x(i)-s
   109  continue
  c
  c******************************
  
negative.code: c     ~~~~~~~~~~~~~~~~~~
  c
  c     p    = array {p(i)}, containing the values of the BLS spectrum
  c     at the i-th frequency value -- the frequency values are
  c     computed as  f = fmin + (i-1)*df
  c     bper = period at the highest peak in the frequency spectrum
  c     bpow = value of {p(i)} at the highest peak
  c     depth= depth of the transit at   *bper*
  