anchor.diff:  
         case(101)
           ! Format 101:
  -        ! # particleID turn s[m] x[mm] xp[mrad] y[mm] yp[mrad] z[mm] (E-E0)/E0[1] ktrack
  -        if(dump_hdf5Format(3) == 0) then
  +        if(dump_hdf5Format(101) == 0) then
             allocate(setFields(19))
  -          setFields(1)  = h5_dataField(name="ID",         type=h5_typeInt)
  -          setFields(2)  = h5_dataField(name="TURN",       type=h5_typeInt)
  -          setFields(3)  = h5_dataField(name="S",          type=h5_typeReal)
  -          setFields(4)  = h5_dataField(name="X",          type=h5_typeReal)
  -          setFields(5)  = h5_dataField(name="XP",         type=h5_typeReal)
  -          setFields(6)  = h5_dataField(name="Y",          type=h5_typeReal)
  -          setFields(7)  = h5_dataField(name="YP",         type=h5_typeReal)
  -          setFields(8)  = h5_dataField(name="dE/E",       type=h5_typeReal)
  -          setFields(9)  = h5_dataField(name="SIGMA",      type=h5_typeReal)
  -          setFields(10) = h5_dataField(name="KTRACK",     type=h5_typeInt)
  -          setFields(11) = h5_dataField(name="E",          type=h5_typeReal)
  -          setFields(12) = h5_dataField(name="PC",         type=h5_typeReal)
  -          setFields(13) = h5_dataField(name="P/P0",       type=h5_typeReal)
  -          setFields(14) = h5_dataField(name="P0/P)",      type=h5_typeReal)
  -          setFields(15) = h5_dataField(name="BETA0/BETA", type=h5_typeReal)
  -          setFields(16) = h5_dataField(name="MASS",       type=h5_typeReal)
  -          setFields(17) = h5_dataField(name="M/M0/Q0/Q",  type=h5_typeReal)
  -          setFields(18) = h5_dataField(name="ENERGY0",    type=h5_typeReal)
  -          setFields(19) = h5_dataField(name="PC0",        type=h5_typeReal)
  -          call h5_createFormat("dumpFormat3", setFields, dump_hdf5Format(3))
  +          setFields(1)  = h5_dataField(name="ID",      type=h5_typeInt)
  +          setFields(2)  = h5_dataField(name="TURN",    type=h5_typeInt)
  +          setFields(3)  = h5_dataField(name="S",       type=h5_typeReal)
  +          setFields(4)  = h5_dataField(name="X",       type=h5_typeReal)
  +          setFields(5)  = h5_dataField(name="XP",      type=h5_typeReal)
  +          setFields(6)  = h5_dataField(name="Y",       type=h5_typeReal)
  +          setFields(7)  = h5_dataField(name="YP",      type=h5_typeReal)
  +          setFields(8)  = h5_dataField(name="dE/E",    type=h5_typeReal)
  +          setFields(9)  = h5_dataField(name="SIGMA",   type=h5_typeReal)
  +          setFields(10) = h5_dataField(name="KTRACK",  type=h5_typeInt)
  +          setFields(11) = h5_dataField(name="E",       type=h5_typeReal)
  +          setFields(12) = h5_dataField(name="PC",      type=h5_typeReal)
  +          setFields(13) = h5_dataField(name="P/P0",    type=h5_typeReal)
  +          setFields(14) = h5_dataField(name="P0/P)",   type=h5_typeReal)
  +          setFields(15) = h5_dataField(name="RV",      type=h5_typeReal)
  +          setFields(16) = h5_dataField(name="MASS",    type=h5_typeReal)
  +          setFields(17) = h5_dataField(name="MTC",     type=h5_typeReal)
  +          setFields(18) = h5_dataField(name="ENERGY0", type=h5_typeReal)
  +          setFields(19) = h5_dataField(name="PC0",     type=h5_typeReal)
  +          call h5_createFormat("dumpFormat101", setFields, dump_hdf5Format(101))
           end if
  -        call h5_createDataSet(dump_fname(i), h5_dumpID, dump_hdf5Format(3), dump_hdf5DataSet(i), napx)
  +        call h5_createDataSet(dump_fname(i), h5_dumpID, dump_hdf5Format(101), dump_hdf5DataSet(i), napx)
   
         end select
   
  
positive.code: 
  #ifdef HDF5
    call alloc(dump_hdf5DataSet,      nele_new,0,        "dump_hdf5DataSet",-1)
    call alloc(dump_hdf5Format,       9,       0,        "dump_hdf5Format")
  #endif
  
  end subroutine dump_expand_arrays
  
negative.code: #endif
      end if
  
    ! ------------------------------------------------------------------
    !  Format #101
    !  Same as fmt 3, but with additional variables
    ! ------------------------------------------------------------------
    else if(fmt == 101) then
      if(i == 0 .and. ix == 0) then
        localDcum   = zero
  