anchor.diff:  
     if(.not. do_coll) return
     ! Arrays that are only needed if Collimation is enabled
  -
  -  call alloc(db_elemMap, nele_new, -1, "db_elemMap")
   
     call alloc(tbetax,  nblz_new, zero, 'tbetax')  !(nblz)
     call alloc(tbetay,  nblz_new, zero, 'tbetay')  !(nblz)
  
positive.code:       if (firstrun) then
          write(all_impacts_unit,'(a)') '# 1=name 2=turn 3=s'
          write(all_absorptions_unit,'(a)') '# 1=name 2=turn 3=s'
          write(FirstImpacts_unit,*)                                                   &
          '%1=name,2=iturn, 3=icoll, 4=nabs, 5=s_imp[m], 6=s_out[m], ',&
          '7=x_in(b!)[m], 8=xp_in, 9=y_in, 10=yp_in, ',                &
          '11=x_out [m], 12=xp_out, 13=y_out, 14=yp_out'
          write(coll_scatter_unit,*) &
          "#1=icoll, 2=iturn, 3=np, 4=nabs (1:Nuclear-Inelastic,2:Nuclear-Elastic,3:pp-Elastic,4:Single-Diffractive,5:Coulomb)", &
          ", 5=dp, 6=dx', 7=dy'"
        end if ! if (firstrun) then
  #ifdef HDF5
      end if
  
negative.code:   use mod_units
  !  use mod_ranecu
    use mod_ranlux
  
  #ifdef HDF5
    use hdf5_output
  