anchor.diff:    call resize(dpsvLast,  npart_new, zero, "dpsvLast")  ! (npart)
     call resize(naaLast,   npart_new, 0, "naaLast")      ! nuclear mass [] (npart)
     call resize(nzzLast,   npart_new, 0 , "nzzLast")     ! atomic number [] (npart)
  -
  -  do i=npart, npart_new
  -    plost(i) = 0
  -  end do
   
   end subroutine aperture_expand_arrays
   
  
positive.code:            pstop(nlostp(j))=.true.
             ! Record for postpr
             if(.not.limifound.or.kape(ix).eq.0) then
                aperv(nlostp(j),1) = aper(1)
                aperv(nlostp(j),2) = aper(2)
             else
                aperv(nlostp(j),1) = min(ape(1,ix),ape(3,ix))
                aperv(nlostp(j),2) = min(ape(2,ix),ape(4,ix))
             end if
             ixv(nlostp(j))     = ix
             xvl(1,nlostp(j))   = xlos(1)
  
negative.code:             end select
  
              ! - aperture at current step
              call interp_aperture( iLast, ixLast, i, ix, kapert, aprr, slos )
  
              ! Check aperture
              if( lapeofftlt(ix).or.lapeofftlt(ixLast) ) then
  