[AnjaMiehe]

Hello Foamers,

I want to code an if condition for a volScalarField and I do get errors in compiling that the if statement does not work.

volScalarField T

volScalarField lambda
(
    IOobject
    (
        "lambda",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    if (T>Tf)? 0 : pow((Tf-T)/(Tf-Tl),-0.35)      
)

As the pow coding of a negative base with a negative exponent does not work, I need a workaround.
Any other coding idea is welcome as well.

Thanks in advance,
Regards,
Anja

[akidess]

Use pos() and neg(): pos(T - Tf)* pow((Tf-T)/(Tf-Tl),-0.35)

[AnjaMiehe]

Hello akidess,

and thank you very much. That was the basic idea I needed. To make it work finally, I changed it a little bit using the magnitude:

 pos(Tf-T)*pow(mag(Tf-T)/(Tf-Tl),-0.35)

in order not to have a negative base.

Thanks again, Anja