anchor.diff:  #endif
       call units_openUnit(unit=scatter_logFile,fileName="scatter_log.dat",formatted=.true.,mode="w",err=fErr,status="replace")
       write(scatter_logFile,"(a)") "# scatter_log"
  -    write(scatter_logFile,"(a1,a8,1x,a8,2(1x,a20),1x,a8,1x,a4,1x,a13,7(1x,a16))") &
  +    write(scatter_logFile,"(a1,a8,1x,a8,2(1x,a20),1x,a8,1x,a4,2(1x,a12),1x,a9,5(1x,a16))") &
         "#","ID","turn",chr_rPad("bez",20),chr_rPad("generator",20),chr_rPad("process",8),&
  -      "lost","t[MeV^2]","dE/E","dP/P","theta[mrad]","phi[rad]","density","prob","stat_corr"
  +      "lost","t[MeV^2]","theta[mrad]","phi[rad]","dE/E","dP/P","density","prob","stat_corr"
       flush(scatter_logFile)
   #ifdef CR
       scatter_logFilePos = 2
  
positive.code:     write(93,"(a)") "SCATTER> INIT kept already opened file 'scatter_log.dat'"
    end if
  #endif
  
    ! Open scatter_summary.dat
    if(scatter_sumFile == -1) call funit_requestUnit("scatter_summary.dat",scatter_sumFile)
  
negative.code:       character(len=:), allocatable :: colUnits(:)
        logical spErr
        integer nSplit
        call chr_split("turn napx bez generator process nScat nLost scatRat scaling",colNames,nSplit,spErr)
        call chr_split("1 1 text text text 1 1 1 1",colUnits,nSplit,spErr)
        call h5_writeDataSetAttr(scatter_sumDataSet,"colNames",colNames)
        call h5_writeDataSetAttr(scatter_sumDataSet,"colUnits",colUnits)
      end block
  