anchor.diff:             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)
  +             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))
  +             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)
  
positive.code:   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
  
  
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
  