anchor.diff:    else
   #endif
       call f_requestUnit(coll_summaryFile,coll_summaryUnit)
  -    call f_open(unit=coll_summaryUnit,file=coll_summaryFile,formatted=.true.,mode="w")
  +    call f_open(unit=coll_summaryUnit,file=coll_summaryFile,formatted=.true.,mode="w",status="replace")
       write(coll_summaryUnit,"(a1,1x,a5,1x,a20,2(1x,a8),2(1x,a15),1x,a6)") "#","icoll",chr_rPad("collname",20),&
         "nimp","nabs","imp_av","imp_sig","length"
       do icoll = 1, cdb_nColl
  
positive.code:   else
  #endif
      call f_requestUnit(coll_efficFile,coll_efficUnit)
      call f_open(unit=coll_efficUnit,file=coll_efficFile,formatted=.true.,mode="w")
      if(n_tot_absorbed /= 0) then
        write(coll_efficUnit,"(a1,1x,a13,6(1x,a15),1x,a8)") "#","rad_sigma",&
          "frac_x","frac_y","frac_r","eff_x","eff_y","eff_r","n_abs"
  
negative.code:     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",status="replace")
      call f_open(unit=coll_allAbsorbUnit, file=coll_allAbsorbFile, formatted=.true.,mode="w",status="replace")
      call f_open(unit=coll_scatterUnit,   file=coll_scatterFile,   formatted=.true.,mode="w",status="replace")
      call f_open(unit=coll_fstImpactUnit, file=coll_fstImpactFile, formatted=.true.,mode="w",status="replace")
      call f_open(unit=coll_impactUnit,    file=coll_impactFile,    formatted=.true.,mode="w",status="replace")
      call f_open(unit=coll_flukImpUnit,   file=coll_flukImpFile,   formatted=.true.,mode="w",status="replace")
      call f_open(unit=coll_flukImpAllUnit,file=coll_flukImpAllFile,formatted=.true.,mode="w",status="replace")
      call f_open(unit=coll_jawProfileUnit,file=coll_jawProfileFile,formatted=.true.,mode="w",status="replace")
  
      write(coll_allImpactUnit,"(a)") "# 1=name 2=turn 3=s"
      write(coll_allAbsorbUnit,"(a)") "# 1=name 2=turn 3=s"
  