[wouter]

Hello All,

how can I do something like:

    volScalarField clipValue(1E-18);
    volVectorField U1_clip
    (
        IOobject
        (
            "U1_clip",
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::AUTO_WRITE
        ),
        alpha1 < clipValue?scalar(0):U1 // clip Ux if alphax <clipValue
    );

This does not work so do I need to loop through all cells?

Best regards
Wouter

[akidess]

pos(alpha1-clipValue)*U1

[wouter]

hello All,
Thank you akidess for your answer. I tried various types (scalar, volScalarField(did not compile), label) for the clipvalue but I probable did not find the right one.
The result for U1_clip are (0 0 0) or (0 0.1 0) the first one is what I expected except that it is on unexpected places (alpha1 > 1E-18) and the other one is no good at all because the value of U at those places is (0.00970724 0.0137527 -0.0250987). Something goes wrong with the subtraction or with the multiplication. Could you please give me a direction where to look.

best regards
Wouter

Edit: it did work only I had to place the statements in write.H not in createFields.H. I used : const dimensionedScalar clipValue("clipValue",dimless,1.0e-8)