anchor.diff:      nnumxv(j)=crnnumxv(j)
       nlostp(j)=crnlostp(j)
       pstop(j)=crpstop(j)
  +    llostp(j)=crllostp(j)
       xv1(j)=crxv(1,j)
       yv1(j)=cryv(1,j)
       xv2(j)=crxv(2,j)
  
positive.code:   integer,             allocatable, public, save :: crnlostp(:)  ! (npart)
  
    logical,             allocatable, public, save :: crpstop(:)   ! (npart)
  
    integer,                          public, save :: crnpart_old = -1
  
  
negative.code:   call alloc(crnnumxv,     npart_new,      0,       "crnnumxv")
    call alloc(crnlostp,     npart_new,      0,       "crnlostp")
    call alloc(crpstop,      npart_new,      .false., "crpstop")
    call alloc(crllostp,     npart_new,      .false., "crllostp")
  
    crnpart_old = npart_new
  
  