anchor.diff:        (crnnumxv(j),  j=1,crnapxo),       &
         (crpartID(j),  j=1,crnapxo),       &
         (crparentID(j),j=1,crnapxo),       &
  +      (crpairID(:,j),j=1,crnapxo),       &
         (crpstop(j),   j=1,crnapxo),       &
         (crxv1(j),     j=1,crnapxo),       &
         (cryv1(j),     j=1,crnapxo),       &
  
positive.code:   integer,          allocatable, private, save :: crnnumxv(:)   ! (npart)
    integer,          allocatable, private, save :: crpartID(:)   ! (npart)
    integer,          allocatable, private, save :: crparentID(:) ! (npart)
  
    logical,          allocatable, private, save :: crpstop(:)    ! (npart)
    logical,          allocatable, private, save :: crllostp(:)   ! (npart)
  
negative.code: 
    ! Recompute the thick arrays
    if(ithick == 1) call synuthck
  
    ! Recompute the map of particle pairs
    call updatePairMap
  
    ! Aperture data
    aperv(1:napxo,1:2) = craperv(1:napxo,1:2)
  