anchor.diff:              iskc=iskc+1
               if(mod(iskc,iskip).ne.0) goto 460
               if((ia-nstart).lt.0) goto 460
  -            if(progrm.eq.'MAD') then !NON-STF only
  -#ifndef STF
  -              c=c*c1e3
  -              d=d*c1e3
  -              e=e*c1e3
  -              f=f*c1e3
  -              g=g*c1e3
  -              p=p*c1e3
  -
  -              if(ntwin.eq.2) then
  -                c1=c1*c1e3
  -                d1=d1*c1e3
  -                e1=e1*c1e3
  -                f1=f1*c1e3
  -                g1=g1*c1e3
  -                p1=p1*c1e3
  -              end if
  -#else
  -        write(lout,*) "ERROR in postpr: program=MAD not valid for STF."
  -        call prror
  -#endif
  -
  +            if(progrm.eq.'MAD') then
  +              write(lerr,"(a)") "POSTPR> ERROR Program=MAD no longer supported"
  +              call prror
               end if
   !--LYAPUNOV
               if(ntwin.eq.2) then
  
positive.code:          h1=real(h1_tmp,fPrec)
           p1=real(p1_tmp,fPrec)
        endif
  #endif
        if(ierro.gt.0) then
          write(lout,10320) nfile
          goto 550
        endif
  #ifdef CR
  #ifndef STF
  !--Increment crbinrecs by 1
        crbinrecs(91-nfile)=crbinrecs(91-nfile)+1
  #else
        crbinrecs(posi1)=crbinrecs(posi1)+1
  #endif
  #endif
        if(ifipa.lt.1) goto 240
        if(progrm.eq.'MAD') then
  #ifndef STF
          c=c*c1e3
          d=d*c1e3
          e=e*c1e3
          f=f*c1e3
          g=g*c1e3
          p=p*c1e3
          if(ntwin.eq.2) then
            c1=c1*c1e3
            d1=d1*c1e3
            e1=e1*c1e3
            f1=f1*c1e3
            g1=g1*c1e3
            p1=p1*c1e3
          endif
  #else
          write(lout,*) "ERROR in postpr: program=MAD not valid for STF."
          call prror
  #endif
        endif
  !--LYAPUNOV
        if(ntwin.eq.2) then
  
negative.code: !--Now check that we have correct number of binrecs
  !--We can do this only if we know binrecs (NOT post-processing only)
        if (binrec.ne.0) then
          if (binrecs(posi1).ne.crbinrecs(posi1)) then
            write(lout,"(a)") "SIXTRACR> ERROR POSTPR Wrong number of binary records"
            write(lout,"(a,i0,a,3(1x,i0))") "SIXTRACR> Particle ",posi1,", binrec/binrecs/crbinrecs ",&
  