anchor.diff:      setFields(5)  = h5_dataField(name="PROCESS", type=h5_typeChar, size=8)
       setFields(6)  = h5_dataField(name="LOST",    type=h5_typeInt)
       setFields(7)  = h5_dataField(name="T",       type=h5_typeReal)
  -    setFields(8)  = h5_dataField(name="dE/E",    type=h5_typeReal)
  -    setFields(9)  = h5_dataField(name="dP/P",    type=h5_typeReal)
  -    setFields(10) = h5_dataField(name="THETA",   type=h5_typeReal)
  -    setFields(11) = h5_dataField(name="PHI",     type=h5_typeReal)
  +    setFields(8)  = h5_dataField(name="THETA",   type=h5_typeReal)
  +    setFields(9)  = h5_dataField(name="PHI",     type=h5_typeReal)
  +    setFields(10) = h5_dataField(name="dE/E",    type=h5_typeReal)
  +    setFields(11) = h5_dataField(name="dP/P",    type=h5_typeReal)
       setFields(12) = h5_dataField(name="DENSITY", type=h5_typeReal)
       setFields(13) = h5_dataField(name="PROB",    type=h5_typeReal)
       setFields(14) = h5_dataField(name="STATCORR",type=h5_typeReal)
  
positive.code: #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))") &
        "#","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"
      flush(scatter_logFile)
  #ifdef CR
      scatter_logFilePos = 2
  
negative.code:     call scatter_writeReport
      write(scatter_sumFile,"(a)") "#  Summary Log"
      write(scatter_sumFile,"(a)") "# ============="
      write(scatter_sumFile,"(a1,a8,1x,a8,2(1x,a20),3(1x,a8),1x,a9,2(1x,a13))") &
        "#","turn","napx",chr_rPad("bez",20),chr_rPad("generator",20),chr_rPad("process",8), &
        "nScat","nLost","scatRat","crossSec[mb]","scaling"
      flush(scatter_sumFile)
  #ifdef CR
      scatter_sumFilePos = scatter_sumFilePos + 3
  