volScalargpuField& p = thermo.p();
const volScalargpuField& T = thermo.T();
const volScalargpuField& psi = thermo.psi();
const volScalargpuField& mu = thermo.mu();

bool inviscid(true);
if (max(mu.primitiveField()) > 0.0)
{
    inviscid = false;
}

PtrList<volScalargpuField> DiffsY(Y.size());
forAll(Y, i)
{
    DiffsY.set
    (
        i,
	new volScalargpuField 
	(
	    IOobject
    	    (
                "DiffsY",
                mesh.time().timeName(),
         	mesh,
         	IOobject::NO_READ,
         	IOobject::NO_WRITE
     	    ),
            devicemesh,
            dimensionedScalar("zero", mu.dimensions(), 0.0),
            zeroGradientFvPatchScalargpuField::typeName
	)
    );    
}
PtrList<volScalargpuField> HsY(Y.size()); 
forAll(Y, i)
{
    HsY.set
    (
        i,
	new volScalargpuField 
	(
	    IOobject
    	    (
                "HsY",
                mesh.time().timeName(),
         	mesh,
         	IOobject::NO_READ,
         	IOobject::NO_WRITE
     	    ),
            devicemesh,
            dimensionedScalar("zero", e.dimensions(), 0.0),
            zeroGradientFvPatchScalargpuField::typeName
	)
    );    
}