anchor.diff:    x = xp*(c_length-s) + x
     z = zp*(c_length-s) + z
     s = c_length
  -  ! CRY ---------------------
  -
  -  ! CRY ---------------------
  +
     nabs = 0
     cry_proc(j) = iProc
     if(iProc == proc_AM) then
  
positive.code:       end if
  
        ! Transform back to particle coordinates with opening and offset
        x   = x + c_aperture/2 + mirror*c_offset
  
        ! Now mirror at the horizontal axis for negative X offset
        x   = mirror * x
        xp  = mirror * xp
  
        ! Last do rotation into collimator frame
        x_in(j)  =  x*cRRot +  z*sRRot
  
negative.code:     lint(j)   = -one
      indiv(j)  = -one
  
      x  = x_in(j)
      xp = xp_in(j)
      z  = y_in(j)
      zp = yp_in(j)
      p  = p_in(j)
  
      ! CRY ---------------------
      dpop   = (p-p0)/p0
  