anchor.diff:            "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,"(a)") "# 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'"
  -        write(profiling_jaw_unit,"(a)") "# 1=icoll, 2=iturn, 3=np, 4=x, 5=xp, 6=y, 7=yp, 8=s, 9=[1:in,2:out,3:back]"
  +        write(profiling_jaw_unit,'("#",3(a7,1x),5(a17,1x),a12)') &
  +             "icoll", "iturn", "np", "x[m]", "xp[]", "y[m]", "yp[]", "s[m]", "[1:in,2:out]"
         end if ! if (firstrun) then
   #ifdef HDF5
       end if
  
positive.code:       call f_requestUnit('all_absorptions.dat', all_absorptions_unit)
        call f_requestUnit('Coll_Scatter.dat', coll_scatter_unit)
        call f_requestUnit('FirstImpacts.dat', FirstImpacts_unit)
        call f_requestUnit('JawProfiles.dat', profiling_jaw_unit)
  
        open(unit=all_impacts_unit, file='all_impacts.dat') !was 46
        open(unit=all_absorptions_unit, file='all_absorptions.dat') !was 47
        open(unit=coll_scatter_unit, file='Coll_Scatter.dat') !was 3998
        open(unit=FirstImpacts_unit, file='FirstImpacts.dat') !was 39
        open(unit=profiling_jaw_unit, file='JawProfiles.dat') !new
  
        if (firstrun) then
          write(all_impacts_unit,'(a)') '# 1=name 2=turn 3=s'
  
negative.code:   integer, allocatable :: nabs_type(:) !(npart)
  !MAY2005
  
    logical linside(napx)
    real(kind=fPrec), allocatable :: x_in(:) !(npart)
    real(kind=fPrec), allocatable :: xp_in(:) !(npart)
    real(kind=fPrec), allocatable :: y_in(:) !(npart)
  