src/lagrangian/molecularDynamics/molecule/mdTools/resetMDFields.H

   1 if (runTime.outputTime())
   2 {
   3     allSpeciesN_RU = List< scalarField >
   4     (
   5         molecules.potential().nIds(),
   6         scalarField(mesh.nCells(), 0.0)
   7     );
   8
   9     allSpeciesM_RU = List< scalarField >
  10     (
  11         molecules.potential().nIds(),
  12         scalarField(mesh.nCells(), 0.0)
  13     );
  14
  15     allSpeciesVelocitySum_RU = List< vectorField >
  16     (
  17         molecules.potential().nIds(),
  18         vectorField(mesh.nCells(), vector::zero)
  19     );
  20
  21     allSpeciesVelocityMagSquaredSum_RU = List< scalarField >
  22     (
  23         molecules.potential().nIds(),
  24         scalarField(mesh.nCells(), 0.0)
  25     );
  26 }