anchor.diff:        if(drift_length.gt.c1m15) then
           linside(j)=.false.
           if(dowrite_impact) then
  -           if ( tiltangle.gt.zero ) then
  -              x_Dump=(x+c_aperture/two+tiltangle*(s+sp))*mirror+c_offset
  -           else
  -              x_Dump=(x+c_aperture/two+tiltangle*(s+sp-c_length))*mirror+c_offset
  -           end if
  -           xpDump=(xp+tiltangle)*mirror
  -           y_Dump=z
  -           ypDump=zp
  -           s_Dump=s+sp+real(j_slices-1,fPrec)*c_length
  -           write(profiling_jaw_unit,'(3(1x,i7),5(1x,e17.9),1x,i1)') &
  -                icoll,iturn,name(j),x_Dump,xpDump,y_Dump,ypDump,s_Dump,2
  +          if ( tiltangle.gt.zero ) then
  +            x_Dump=(x+c_aperture/two+tiltangle*(s+sp))*mirror+c_offset
  +          else
  +            x_Dump=(x+c_aperture/two+tiltangle*(s+sp-c_length))*mirror+c_offset
  +          end if
  +          xpDump=(xp+tiltangle)*mirror
  +          y_Dump=z
  +          ypDump=zp
  +          s_Dump=s+sp+real(j_slices-1,fPrec)*c_length
  +          write(coll_jawProfileUnit,"(3(1x,i7),5(1x,e17.9),1x,i1)") &
  +            icoll,iturn,name(j),x_Dump,xpDump,y_Dump,ypDump,s_Dump,2
           end if
           if(iexact) then
             zpj = sqrt(one-xp**2-zp**2)
  
positive.code:     call f_open(unit=coll_impactUnit,    file=coll_impactFile,    formatted=.true.,mode="w")
      call f_open(unit=coll_flukImpUnit,   file=coll_flukImpFile,   formatted=.true.,mode="w")
      call f_open(unit=coll_flukImpAllUnit,file=coll_flukImpAllFile,formatted=.true.,mode="w")
  
      write(coll_allImpactUnit,"(a)") "# 1=name 2=turn 3=s"
      write(coll_allAbsorbUnit,"(a)") "# 1=name 2=turn 3=s"
  
negative.code:     call f_requestUnit(coll_impactFile,    coll_impactUnit)
      call f_requestUnit(coll_flukImpFile,   coll_flukImpUnit)
      call f_requestUnit(coll_flukImpAllFile,coll_flukImpAllUnit)
      call f_requestUnit(coll_jawProfileFile,coll_jawProfileUnit)
  
      call f_open(unit=coll_allImpactUnit, file=coll_allImpactFile, formatted=.true.,mode="w")
      call f_open(unit=coll_allAbsorbUnit, file=coll_allAbsorbFile, formatted=.true.,mode="w")
  