anchor.diff:          ! first time particle hits collimator: entering jaw
           linside(j)=.true.
           if(dowrite_impact) then
  -           if ( tiltangle.gt.zero ) then
  -              x_Dump=(x+c_aperture/two+tiltangle*sp)*mirror+c_offset
  -           else
  -              x_Dump=(x+c_aperture/two+tiltangle*(sp-c_length))*mirror+c_offset
  -           end if
  -           xpDump=(xp+tiltangle)*mirror
  -           y_Dump=z
  -           ypDump=zp
  -           s_Dump=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,1
  +          if ( tiltangle.gt.zero ) then
  +            x_Dump=(x+c_aperture/two+tiltangle*sp)*mirror+c_offset
  +          else
  +            x_Dump=(x+c_aperture/two+tiltangle*(sp-c_length))*mirror+c_offset
  +          end if
  +          xpDump=(xp+tiltangle)*mirror
  +          y_Dump=z
  +          ypDump=zp
  +          s_Dump=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,1
           end if
         end if
   !JUNE2005
  
positive.code:     call f_close(coll_flukImpAllUnit)
      call f_close(coll_scatterUnit)
      call f_close(coll_fstImpactUnit)
    endif
    call f_close(profiling_jaw_unit)
  
    call f_requestUnit('amplitude.dat', amplitude_unit)
    call f_requestUnit('amplitude2.dat', amplitude2_unit)
  
negative.code:   integer, private, save :: coll_twissLikeUnit  = -1
    integer, private, save :: coll_sigmaSetUnit   = -1
    integer, private, save :: coll_settingsUnit   = -1
    integer, private, save :: coll_jawProfileUnit = -1
  
    integer, private, save :: distsec_unit, efficiency_unit, efficiency_dpop_unit
    integer, private, save :: coll_summary_unit, amplitude_unit, amplitude2_unit, betafunctions_unit, orbitchecking_unit
    integer, private, save :: efficiency_2d_unit
    integer, private, save :: outlun
  
  #ifdef HDF5
  