anchor.diff:        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)
  +      call f_requestUnit('jaw_profiles.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
  +      open(unit=profiling_jaw_unit, file='jaw_profiles.dat') !new
   
         if (firstrun) then
           write(all_impacts_unit,'(a)') '# 1=name 2=turn 3=s'
  
positive.code:                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)') &
                      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
  
  
negative.code: !                  endif
  !CB
  
        linside(1:napx)=.false.
        
        if(n_slices.gt.1 .and. totals.gt.smin_slices .and. totals.lt.smax_slices .and. &
   &      (cdb_cNameUC(icoll)(1:4).eq.'TCSG' .or. cdb_cNameUC(icoll)(1:3).eq.'TCP' .or. cdb_cNameUC(icoll)(1:4).eq.'TCLA'.or. &
   &       cdb_cNameUC(icoll)(1:3).eq.'TCT' .or. cdb_cNameUC(icoll)(1:4).eq.'TCLI'.or. cdb_cNameUC(icoll)(1:4).eq.'TCL.'.or.  &
  