anchor.diff:  
   ! Log output energy + nucleons as per the FLUKA coupling
       nnuc1       = nnuc1 + naa(j)                          ! outcoming nucleons
  -    ien1        = ien1  + rcp(j) * c1e3                   ! outcoming energy
  +    ien1        = ien1  + (rcp(j) * c1e3)                 ! outcoming energy
   
   ! Fix hits
   ! if(part_hit_pos(j) .eq.ie .and. part_hit_turn(j).eq.iturn)
  
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
  