anchor.diff:  
     partID(1:napxo)   = crpartID(1:napxo)
     parentID(1:napxo) = crparentID(1:napxo)
  +  pairID(:,1:napxo) = crpairID(:,1:napxo)
     pstop(1:napxo)    = crpstop(1:napxo)
     llostp(1:napxo)   = crllostp(1:napxo)
   
  
positive.code:       (crnnumxv(j),  j=1,crnapxo),       &
        (crpartID(j),  j=1,crnapxo),       &
        (crparentID(j),j=1,crnapxo),       &
        (crpstop(j),   j=1,crnapxo),       &
        (crxv1(j),     j=1,crnapxo),       &
        (cryv1(j),     j=1,crnapxo),       &
  
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)
  