anchor.diff:            call root_EnergyDeposition(fluka_geo_index(fluka_ix), nnuc0-nnuc1, c1m3*(ien0-ien1))
           end if
   #endif
  -      ! hisix debugging:
  -      ! write out the particle distribution after the primary
  -        if(fluka_geo_index(fluka_ix).eq.11) then
  -          do j=1,napx
  -            write(unit210,*) naa(j), nzz(j), nucm(j),ejfv(j),mtc(j),dpsv(j)
  -          end do
  -        end if
         end if
   
   !     hisix: check which particle ids have not been sent back
  
positive.code:             pid_q = one
            end if
          end do
  
          if(pid_q.eq.zero.and.pids(j).ne.zero) then
            write(unit209,*) fluka_geo_index(fluka_ix), pids(j)
          end if
        end do
  
  !     empty places
  
negative.code:             pid_q = one
            end if
          end do
        end do
  
  !     empty places
  