anchor.diff:      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")
  -    call f_open(unit=coll_allAbsorbUnit, file=coll_allAbsorbFile, formatted=.true.,mode="w")
  -    call f_open(unit=coll_scatterUnit,   file=coll_scatterFile,   formatted=.true.,mode="w")
  -    call f_open(unit=coll_fstImpactUnit, file=coll_fstImpactFile, formatted=.true.,mode="w")
  -    call f_open(unit=coll_impactUnit,    file=coll_impactFile,    formatted=.true.,mode="w")
  -    call f_open(unit=coll_flukImpUnit,   file=coll_flukImpFile,   formatted=.true.,mode="w")
  -    call f_open(unit=coll_flukImpAllUnit,file=coll_flukImpAllFile,formatted=.true.,mode="w")
  -    call f_open(unit=coll_jawProfileUnit,file=coll_jawProfileFile,formatted=.true.,mode="w")
  +    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"
  
positive.code: 
    ! Write orbitchecking.dat
    call f_requestUnit(coll_orbitCheckFile,coll_orbitCheckUnit)
    call f_open(unit=coll_orbitCheckUnit,file=coll_orbitCheckFile,formatted=.true.,mode="w")
    write(coll_orbitCheckUnit,"(a1,1x,a6,3(1x,a15))") "#","s","torbitx","torbity"
    do i=1,iu
      write(coll_orbitCheckUnit,"(i8,3(1x,1pe15.7))") i, dcum(i), torbx(i), torby(i)
  
negative.code:   else
  #endif
      call f_requestUnit(coll_efficDPFile,coll_efficDPUnit)
      call f_open(unit=coll_efficDPUnit,file=coll_efficDPFile,formatted=.true.,mode="w",status="replace")
      if(n_tot_absorbed /= 0) then
        write(coll_efficDPUnit,"(a1,1x,a13,2(1x,a15),2(1x,a8))") "#","dp/p","n_dpop/tot_nabs","n_dpop","tot_nabs","npart"
        do k=1,numeffdpop
  