anchor.diff:  
   #ifdef CR
   !--   Initiate count of binary records
  -#ifndef STF
  -      crbinrecs(91-nfile)=1
  -#else
         crbinrecs(posi1)=1
   #endif
  -#endif
   
    210  ifipa=0
  -#ifndef STF
  -      if(ntwin.eq.1) then
  -         read(nfile,end=530,iostat=ierro) ia,ifipa,b_tmp,c_tmp,d_tmp,e_tmp,f_tmp,g_tmp,h_tmp,p_tmp
  -         b=real(b_tmp,fPrec)
  -         c=real(c_tmp,fPrec)
  -         d=real(d_tmp,fPrec)
  -         e=real(e_tmp,fPrec)
  -         f=real(f_tmp,fPrec)
  -         g=real(g_tmp,fPrec)
  -         h=real(h_tmp,fPrec)
  -         p=real(p_tmp,fPrec)
  -
  -      elseif(ntwin.eq.2) then
  -         read(nfile,end=530,iostat=ierro) ia,ifipa,b_tmp,c_tmp,d_tmp,e_tmp,f_tmp,g_tmp,h_tmp,p_tmp,ilapa,b_tmp,c1_tmp,d1_tmp,    &
  -     &e1_tmp,f1_tmp,g1_tmp,h1_tmp,p1_tmp
  -
  -         b=real(b_tmp,fPrec)
  -         c=real(c_tmp,fPrec)
  -         d=real(d_tmp,fPrec)
  -         e=real(e_tmp,fPrec)
  -         f=real(f_tmp,fPrec)
  -         g=real(g_tmp,fPrec)
  -         h=real(h_tmp,fPrec)
  -         p=real(p_tmp,fPrec)
  -
  -         c1=real(c1_tmp,fPrec)
  -         d1=real(d1_tmp,fPrec)
  -         e1=real(e1_tmp,fPrec)
  -         f1=real(f1_tmp,fPrec)
  -         g1=real(g1_tmp,fPrec)
  -         h1=real(h1_tmp,fPrec)
  -         p1=real(p1_tmp,fPrec)
  -      endif
  -#else
  -!     STF case: read tracking data until one reaches right particle.
  +      ! read tracking data until one reaches right particle.
         if(ntwin.eq.1) then
            read(nfile,end=530,iostat=ierro) ia_stf,ifipa_stf,b_tmp,c_tmp,d_tmp,e_tmp,f_tmp,g_tmp,h_tmp,p_tmp
   
  
positive.code: 
            rewind nfile
            !Skip headers
  #ifdef STF
            do j=1,itopa,2
  #endif
            read(nfile,iostat=ierro)
  #ifdef STF
            enddo
  #endif
            if(ierro.gt.0) then
              write(lout,10320) nfile
  #ifdef CR
  
negative.code:          h1=real(h1_tmp,fPrec)
           p1=real(p1_tmp,fPrec)
        endif
        if(ierro.gt.0) then
          write(lout,10320) nfile
          goto 550
        endif
  #ifdef CR
        crbinrecs(posi1)=crbinrecs(posi1)+1
  #endif
        if(ifipa.lt.1) goto 240
        if(progrm.eq.'MAD') then
          write(lerr,"(a)") "POSTPR> ERROR Program=MAD no longer supported"
          call prror
        endif
  !--LYAPUNOV
        if(ntwin.eq.2) then
  