OpenFOAM-4.x-lab/applications/solvers/combustion/NOxFoam_thermalNOx/createFields.H

volScalarField T
(
    IOobject
    (
        "T",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

volScalarField rho
(
    IOobject
    (
        "rho",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

Info<< "\nReading field U\n" << endl;
volVectorField U
(
    IOobject
    (
        "U",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

Info<< "Reading/calculating face flux field phi\n" << endl;
surfaceScalarField phi
(
    IOobject
    (
        "phi",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    linearInterpolate(rho*U) & mesh.Sf()
);

volScalarField nut
(
    IOobject
    (
        "nut",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

volScalarField mut
(
    IOobject
    (
        "mut",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    dimensionedScalar("mut",dimensionSet(1,-1,-1,0,0,0,0),0.0)
);

volScalarField Xno
(
    IOobject
    (
        "Xno",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT,
        IOobject::NO_WRITE
    ),
    mesh, 
    dimensionedScalar("Xno",dimensionSet(0,-3,0,0,1,0,0),0.0)
);

volScalarField Xn2
(
    IOobject
    (
        "Xn2",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::NO_WRITE
    ),
    mesh,
    dimensionedScalar("Xn2",dimensionSet(0,-3,0,0,1,0,0),0.0)
);

volScalarField Xo2
(
    IOobject
    (
        "Xo2",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::NO_WRITE
    ),
    mesh,
    dimensionedScalar("Xo2",dimensionSet(0,-3,0,0,1,0,0),0.0)
);

volScalarField Xh2o
(
    IOobject
    (
        "Xh2o",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::NO_WRITE
    ),
    mesh,
    dimensionedScalar("Xh2o",dimensionSet(0,-3,0,0,1,0,0),0.0)
);

volScalarField Xo
(
    IOobject
    (
        "Xo",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::NO_WRITE
    ),
    mesh,
    dimensionedScalar("Xo",dimensionSet(0,-3,0,0,1,0,0),0.0)
);

volScalarField Xoh
(
    IOobject
    (
        "Xoh",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::NO_WRITE
    ),
    mesh,
    dimensionedScalar("Xoh",dimensionSet(0,-3,0,0,1,0,0),0.0)
);

volScalarField NO
(
    IOobject
    (
        "NO",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

volScalarField N2
(
    IOobject
    (
        "N2",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
   mesh
);

volScalarField O2
(
    IOobject
    (
        "O2",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

volScalarField H2O
(
    IOobject
    (
        "H2O",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

volScalarField O
(
    IOobject
    (
        "O",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

volScalarField OH
(
    IOobject
    (
        "OH",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

volScalarField Sno
(
    IOobject
    (
        "Sno",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    dimensionedScalar("Sno",dimensionSet(1,-3,-1,0,0,0,0),0.0)
);

volScalarField kf1
(
    IOobject
    (
        "kf1",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    dimensionedScalar("kf1",dimensionSet(0,3,-1,0,-1,0,0),0.0)
);

volScalarField kf2
(
    IOobject
    (
        "kf2",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    dimensionedScalar("kf2",dimensionSet(0,3,-1,0,-1,0,0),0.0)
);

volScalarField kf3
(
    IOobject
    (
        "kf3",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    dimensionedScalar("kf3",dimensionSet(0,3,-1,0,-1,0,0),0.0)
);

volScalarField kr1
(
    IOobject
    (
        "kr1",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    dimensionedScalar("kr1",dimensionSet(0,3,-1,0,-1,0,0),0.0)
);

volScalarField kr2
(
    IOobject
    (
        "kr2",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    dimensionedScalar("kr2",dimensionSet(0,3,-1,0,-1,0,0),0.0)
);

IOdictionary modelParameter
(
    IOobject
    (
        "modelParameter",
        runTime.constant(),
        mesh,
        IOobject::MUST_READ_IF_MODIFIED,
        IOobject::NO_WRITE
    )
);

bool instantaneousRadicals(modelParameter.lookupOrDefault("instantaneousRadicals", true));