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:       real*8, dimension(30000) :: y
        integer, dimension(30000) :: ibi
        integer :: minbin,nbmax,kmi,kma,kkmi,nb1,nbkma,i,jf,j,jnb,k,kk,
       &nb2,jn1,jn2
        real*8 :: rn,tot,s,t1,f0,p0,ph,power,rn1,pow,rn3,s3
  
  c
  
negative.code: c     t    = array {t(i)}, containing the time values of the time series
  c     x    = array {x(i)}, containing the data values of the time series
  c     u    = temporal/work/dummy array, must be dimensioned in the
  c     calling program in the same way as  {t(i)}
  c     v    = the same as  {u(i)}
  c     nf   = number of frequency points in which the spectrum is computed
  c     fmin = minimum frequency (MUST be > 0)
  