anchor.diff:              ph     = u(i)*f0
               ph     = ph-idint(ph)
   c
  -c     I think this is binning things in phasebins
  +c     This is binning things in phasebins
               j      = 1 + idint(nbins*ph)
               ibi(j) = ibi(j) + 1
               y(j) =   y(j) + v(i)
  
positive.code:          power=0.0d0
  c
           do 1 i=1,nbins
              s     = 0.0d0
              k     = 0
              kk    = 0
              nb2   = i+kma
              do 2 j=i,nb2
                 k     = k+1
                 kk    = kk+ibi(j)
                 s     = s+y(j)
  
negative.code:       if(nbins.gt.nbinsmax) write(*,*) ' NBINS > NBINSMAX !!'
        if(nbins.gt.nbinsmax) stop
  c
  c     Check if the time-series length is larger than the minimum test freq
        tot=t(n)-t(1)
        if(fmin.lt.1.0d0/tot) write(*,*) ' fmin < 1/T !!'
        if(fmin.lt.1.0d0/tot) stop
  