anchor.diff:        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     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
  +c     Also rescales the times to the tprime = time - t1
         do 103 i=1,n
            u(i)=t(i)-t1
            s=s+x(i)
  
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:                jn2   = j
                 rn3   = rn1
                 s3    = s
  c
   2          continue
  c
   1       continue
  c
           power = dsqrt(power)
  