anchor.diff:        call root_EnergyDeposition(icoll, nnuc0-nnuc1,c1m3*(ien0-ien1))
       end if
   #endif
  -    write(unit208,"(2(i5,1x),e24.16)") icoll, (nnuc0-nnuc1), c1m3*(ien0-ien1)
  +    write(unit208,"(2(i6,1x),e24.16)") icoll, (nnuc0-nnuc1), c1m3*(ien0-ien1)
       flush(unit208)
     end if
   
  
positive.code:   call k2coll_merlinInit
  #endif
  
  #ifdef G4COLLIMATION
  
    if(n_slices /= 0) then
      write(lerr,"(a)") "COLL> ERROR Cannot use jaw fit in G4COLLIMATION version of SixTrack"
      call prror
    end if
  
    ! Open the edep file
    if(unit208 == -1) then
      call f_requestUnit(fort208,unit208)
      call f_open(unit=unit208,file=fort208,formatted=.true.,mode="w",status="replace")
    end if
  
    ! This function lives in the G4Interface.cpp file in the g4collimat folder
  
negative.code:     call f_requestUnit(coll_tracksFile,coll_tracksUnit)
      call f_open(unit=coll_tracksUnit,file=coll_tracksFile,formatted=.true.,mode="w",status="replace")
      write(coll_tracksUnit,"(a)") "# name turn s x xp y yp DE/E type"
    end if
  
    call f_requestUnit(coll_pencilFile,coll_pencilUnit)
    call f_open(unit=coll_pencilUnit, file=coll_pencilFile,formatted=.true.,mode="w",status="replace")
    write(coll_pencilUnit,"(a)") "# x xp y yp"
  
    ! Crystal Files
    if(coll_hasCrystal .and. dowrite_crycoord) then
  