anchor.diff:        xv2(j)   * sqrt(emitY)*tas(1,3,3) + &
         yv2(j)   * sqrt(emitY)*tas(1,3,4) + &
         sigmv(j) * sqrt(emitZ)*tas(1,3,5) + &
  -      ejv(j)   * sqrt(emitZ)*tas(1,3,6)*c1m3
  +      ejv(j)   * sqrt(emitZ)*tas(1,3,6) * c1m3
   
       tmpYP = &
         xv1(j)   * sqrt(emitX)*tas(1,4,1) + &
  
positive.code:       xv2(j)   * sqrt(emitY)*tas(1,4,3) + &
        yv2(j)   * sqrt(emitY)*tas(1,4,4) + &
        sigmv(j) * sqrt(emitZ)*tas(1,4,5) + &
        ejv(j)   * sqrt(emitZ)*tas(1,4,6)*c1m3
  
      tmpS = &
        xv1(j)   * sqrt(emitX)*tas(1,5,1) + &
  
negative.code:   real(kind=fPrec) :: emitX, emitY, emitZ
    integer          :: j
  
    if(iclo6 == 0) then
      write(lout,"(a)") "COLLDIST> ERROR The 6D closed orbit (iclo6 > 0) is required for format 6."
      call prror
    endif
  
    emitX = cdist_emitX
    emitY = cdist_emitY
    emitZ = cdist_bunchLen*c1m3 * cdist_spreadE
  
    do j=1,napx
  
      tmpX = &
        xv1(j)   * sqrt(emitX)*tas(1,1,1) + &
  