resetMDFields.H 624 Bytes
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
if (runTime.outputTime())
{
    allSpeciesN_RU = List< scalarField >
    (
        molecules.pairPotentials().nIds(),
        scalarField (mesh.nCells(), 0.0)
    );

    allSpeciesM_RU = List< scalarField >
    (
        molecules.pairPotentials().nIds(),
        scalarField (mesh.nCells(), 0.0)
    );

    allSpeciesVelocitySum_RU = List< vectorField >
    (
        molecules.pairPotentials().nIds(),
        vectorField (mesh.nCells(), vector::zero)
    );

    allSpeciesVelocityMagSquaredSum_RU = List< scalarField >
    (
        molecules.pairPotentials().nIds(),
        scalarField (mesh.nCells(), 0.0)
    );
}