[Teresa Sun]

The issue happens while implementing an eddy-viscosity model that is based on Smagorinsky model.

We want to compare the dimensioned volScalarField(m^2/s^2) to a constant in .C file, for instance, Smagorinsky.C.

template<class BasicTurbulenceModel>		
tmp<volScalarField> QR<BasicTurbulenceModel>::k 
(			
    const tmp<volTensorField>& gradU		
) const
{                                             
    volSymmTensorField  D(symm(gradU));       
    volScalarField 	r(max(-det(D),dimensionedScalar("r",dimensionSet(0, 0, -3, 0, 0, 0, 0), 0.0)));
    new volScalarField     
    (
      IOobject
          (
             "q1",
             this->runTime_.timeName(),
             this->mesh_,
             IOobject::NO_READ,
             IOobject::NO_WRITE
          ),
             this->mesh_,
             dimensionedScalar("q1", dimensionSet(0,0,-2,0,0,0,0), 1e-6)
     );
    if ( q < q1 )  
    {         
    volScalarField	r1(0.0*(-det(D)));
    return tmp<volScalarField>
	(
	new	volScalarField		
	(
		IOobject
		(
			IOobject::groupName("k",this->alphaRhoPhi_.group()),
			this->runTime_.timeName(),
			this->mesh_
		),
		r1		                     
	);
    }	

    else
    {
     return tmp<volScalarField>
      (
	new volScalarField	
	(
		IOobject
		(
			IOobject::groupName("k", this->alphaRhoPhi_.group()),
			this->runTime_.timeName(),		
			this->mesh_				
		),
		mag(r)/q				
	)
      );
     }	
}

The error message is

'q1' was not declared in this scope

Any point would be great appreciated!
Teresa