anchor.diff:          lint(j) = zlm - drift_length
         end if
       end if
  +
       !++  Transform back to particle coordinates with opening and offset
       if(x < 99.0d-3) then
         ! Include collimator tilt
  
positive.code:   real(kind=fPrec) x_Dump,xpDump,y_Dump,ypDump,s_Dump
    real(kind=fPrec) cRot,sRot,cRRot,sRRot
    real(kind=fPrec) s_impact,xinn,xpinn,yinn,ypinn
    real(kind=fPrec) x_in0,xp_in0,s_cry
  
    ! Initilaisation
  
  
negative.code:     z      = y_in(j)
      zp     = yp_in(j)
      p      = p_in(j)
      sp     = zero
      dpop   = (p - p0)/p0
      x_flk  = zero
  