anchor.diff:                 y_Dump=rcy (j)*cos_mb(c_rotation)-sin_mb(c_rotation)*rcx (j)
                  ypDump=rcyp(j)*cos_mb(c_rotation)-sin_mb(c_rotation)*rcxp(j)
                  s_Dump=c_length
  -               write(profiling_jaw_unit,'(i2,1x,i7,1x,i5,5(1x,e17.9),1X,i1)') &
  +               write(profiling_jaw_unit,'(3(1x,i7),5(1x,e17.9),1x,i1)') &
                       icoll,iturn,flukaname(j),x_Dump,xpDump,y_Dump,ypDump,s_Dump,2
               end if
            end do
         end if
       end if !if(cdb_cNameUC(icoll)(1:4).eq.'COLM') then
  -    call dealloc(linside, "linside")
     end if !if (found) then
   end subroutine collimate_do_collimator
   
  
positive.code:           "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]"
        end if ! if (firstrun) then
  #ifdef HDF5
      end if
  
negative.code:            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
          end if
        end if
  