anchor.diff:        fluka_ix = -1
         fluka_nturn = -1
         flush(unit208)
  -      flush(unit209)
         return
   end subroutine kernel_fluka_element
   
  
positive.code:       fluka_ix = -1
        fluka_nturn = -1
        flush(unit208)
        flush(unit209)
  
        return
  end subroutine kernel_fluka_exit
  
  
negative.code:           call root_EnergyDeposition(fluka_geo_index(fluka_ix), nnuc0-nnuc1, c1m3*(ien0-ien1))
          end if
  #endif
        end if
  
  !     hisix: check which particle ids have not been sent back
  