anchor.diff:  c     Start period search     *
   c******************************
   c
  -      do 100 jf=1,nf
  +      do 100 jf=1,nfreqs
   c
            f0=fmin+df*dfloat(jf-1)
            p0=1.0d0/f0
  
positive.code: 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)
  c     df   = frequency step
  c     nbins   = number of bins in the folded time series at any test period
  
negative.code:       if(kmi.lt.1) kmi=1
        kma=idint(qmax*dfloat(nbins))+1
        kkmi=idint(rn*qmin)
        if(kkmi.lt.minbinpts) kkmi=minbinpts
        bpow=0.0d0
  c
  c     The following variables are defined for the extension
  