[Linse]

Hello All,

I would like to implement an additional computation into some solvers. Basically, an additional action should be taken if a certain testvalue has been passed.
Basically it works, but only with the threshold value hardcoded within the code.

In order to avoid recompiling the specific solver every time the threshold value is changed, I would like to make a dictionary file such as the "transportProperties" residing in the constant-folder.

How can I get the values inside such files into my solver for comparison reason?

For example it works to write:

if (p[celli] > scalar(12))
{
action that is desired;
}

What do I have to write instead of "scalar(12)" for the routine to take over the value "testvalue" from the dictionary?

What I have working by now (in simplefied form) are following entries:

Within the solver I have:

forAll(mesh.cells(),celli)
{
if (p[celli] > scalar(12))
{action1;}
else
{action2;}
}

Within my createFields.H I have

IOdictionary dictionaryName
(
    IOobject
    (
        "dictionaryName",
        runTime.constant(),
        mesh,
        IOobject::MUST_READ,
        IOobject::NO_WRITE
    )
);

dimensionedScalar  testvalue(dictionaryName.lookup("testvalue"));

As I said: The first snippet works fine, the second complains if there is no file "dictionaryName" within the constant-folder and is fine if that is there. But all my tries in taking over a value from the dictionary-file failed so far.
Compilation is fine, though...

[andyru]

Hi,

for example put in constant/myProperties
put in there:
myconst myconst [0 0 0 0 0 0 0] 1.0;

and read in with:
Info<< "\nReading myProperites" << endl;

IOdictionary myProperties
(
IOobject
(
"myProperties",
runTime.constant(),
mesh,
IOobject::MUST_READ,
IOobject::NO_WRITE
)
);
const dimensionedScalar myconst(myProperties.lookup("myconst"));

hope, this helps!

[Linse]

Hi Andyru!

Thanks for the tip, but up to now that did not succeed...
What I had in the createFields.H is:

IOdictionary evaluationProperties
(
IOobject
(
"evaluationProperties",
runTime.constant(),
mesh,
IOobject::MUST_READ,
IOobject::NO_WRITE
)
);
const dimensionedScalar evalp(evaluationProperties.lookup("evalp"));

This far it compiles fine.

But as soon as I put in the evaluation-routine where it should check if the pressure is higher as evalp, it throws a compilation error as below:

forAll(mesh.cells(),celli)
{
if (p[celli] > evalp )
{action1;}
else
{action2;}
}

throws:
error: no match for 'operator>' in 'p.Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::<anonymous>.Foam:imensionedField <double, Foam::volMesh>::anonymous>.Foam::Field<double>::<a nonymous>,Foam::List<double>::<anonymous>.Foam::UL ist<T>:perator[] [with T = double, Foam::label = int] (celli) > evalp'

I would have guesse it is a simple thing just to find my mistake, but I guess I am always again missing to take the wheel before pushing down the accelerator in front of the wall...

Thanks for any further suggestions!

[marupio]

if (p[celli] > evalp.value() )

[Linse]

Thank you, David!
That did the trick! :-)