anchor.diff:      fluka_clight    = clight
   
       call alloc(pids,               npart, 0, "pids")
  -! These are now "global" ID variables and should be handled in other modules
  -!    call alloc(partID,             npart, 0, 'partID')
  -!    call alloc(parentID,           npart, 0, 'parentID')
  -!    call alloc(fluka_weight,       npart, one, 'fluka_weight')
       call alloc(fluka_type,         nele, FLUKA_NONE, 'fluka_type')
       call alloc(fluka_geo_index,    nele, 0, 'fluka_geo_index')
       call alloc(fluka_synch_length, nele, zero, 'fluka_synch_length')
  -
  -! See above global comment
  -!    do j = 1, npart
  -!      partID(j) = j
  -!      parentID(j) = j
  -!    end do
  -
  -!    fluka_weight       = one
  -!    fluka_type         = FLUKA_NONE
  -!    fluka_geo_index    = 0
  -!    fluka_synch_length = zero
   
       if(unit208 == -1) then
         call f_requestUnit(fort208,unit208)
  
positive.code:   subroutine fluka_mod_end()
      implicit none
      call dealloc(pids,"pids")
  !    call dealloc(partID,'partID')
  !    call dealloc(parentID,'parentID')
  !    call dealloc(fluka_weight,'fluka_weight')
      call dealloc(fluka_type,'fluka_type')
      call dealloc(fluka_geo_index,'fluka_geo_index')
      call dealloc(fluka_synch_length,'fluka_synch_length')
  
negative.code:     integer :: npart_new, nele_new, j
  
      call alloc(pids,               npart_new, 0, "pids")
      call alloc(fluka_type,         nele_new, FLUKA_NONE, 'fluka_type')
      call alloc(fluka_geo_index,    nele_new, 0, 'fluka_geo_index')
      call alloc(fluka_synch_length, nele_new, zero, 'fluka_synch_length')
  
      fluka_max_npart = npart_new
  
  