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:       s=0.0d0
  c     t1 is the first elem in time array used as arbitrary epoch
        t1=t(1)
  c     These two loops below renormalize the magseries to 0.0
  c     by calculating the average of the magseries
  c     Also rescales the times to the tprime = time - t1
        do 103 i=1,n
           u(i)=t(i)-t1
           s=s+x(i)
  