//Info<< "Reading thermophysical properties\n" << endl;

//autoPtr<psiThermo> pThermo
//(
//    psiThermo::New(mesh)
//);
//psiThermo& thermo = pThermo();

//volScalarField& e = thermo.he();

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

//volScalarField rho
//(
//    IOobject
//    (
//        "rho",
//        runTime.timeName(),
//        mesh,
//        IOobject::NO_READ,
//        IOobject::AUTO_WRITE
//    ),
//    thermo.rho()
//);

//volVectorField rhoU
//(
//    IOobject
//    (
//        "rhoU",
//        runTime.timeName(),
//        mesh,
//        IOobject::NO_READ,
//        IOobject::NO_WRITE
//    ),
//    rho*U
//);

//volScalarField rhoE
//(
//    IOobject
//    (
//        "rhoE",
//        runTime.timeName(),
//        mesh,
//        IOobject::NO_READ,
//        IOobject::NO_WRITE
//    ),
//    rho*(e + 0.5*magSqr(U))
//);

//surfaceScalarField pos
//(
//    IOobject
//    (
//        "pos",
//        runTime.timeName(),
//        mesh
//    ),
//    mesh,
//    dimensionedScalar("pos", dimless, 1.0)
//);

//surfaceScalarField neg
//(
//    IOobject
//    (
//        "neg",
//        runTime.timeName(),
//        mesh
//    ),
//    mesh,
//    dimensionedScalar("neg", dimless, -1.0)
//);

//surfaceScalarField phi("phi", fvc::flux(rhoU));

//Info<< "Creating turbulence model\n" << endl;
//autoPtr<compressible::turbulenceModel> turbulence
//(
//    compressible::turbulenceModel::New
//    (
//        rho,
//        U,
//        phi,
//        thermo
//    )
//);
