Info<< "Reading field T\n" << endl;

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

    // Add overset specific interpolations
    {
        dictionary oversetDict;
        oversetDict.add("T", true);

        const_cast<dictionary&>
        (
            mesh.schemesDict()
        ).add
        (
            "oversetInterpolationRequired",
            oversetDict,
            true
        );
    }


    Info<< "Reading transportProperties\n" << endl;

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


    Info<< "Reading diffusivity DT\n" << endl;

    dimensionedScalar DT("DT", dimViscosity, transportProperties);