anchor.diff:                 jn2   = j
                  rn3   = rn1
                  s3    = s
  +c
    2          continue
  +c
    1       continue
   c
            power = dsqrt(power)
  
positive.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
  
negative.code:             ph     = u(i)*f0
              ph     = ph-idint(ph)
  c
  c     This is binning things in phasebins
              j      = 1 + idint(nbins*ph)
              ibi(j) = ibi(j) + 1
              y(j) =   y(j) + v(i)
  