[marupio]

I don't get it. When I use the "basePoint" / "normal" to specify a plane, it fails silently whenever the Y-coordinate of "normal" is zero. That is, it runs, reports all the times, and completes, but there are no files in the surfaces/[time]/ directories. Here's my sampleDict:

setFormat raw;
surfaceFormat raw;

interpolationScheme cellPointFace;

fields
(
U
);

sets
(
);

surfaces
(
newPlane
{
type plane;
basePoint (5 2.5 5);
normalVector (0 0 1);
interpolate true;
triangulate true;
}
);

This produces no output. If I change the normalVector to:

normalVector (0 0.00001 1);

it succeeds.

Does anyone have any idea what is going on here?

(OpenFOAM 1.5-dev)

[olesen]

Originally Posted by marupio

I don't get it. When I use the "basePoint" / "normal" to specify a plane, it fails silently whenever the Y-coordinate of "normal" is zero.

...
This produces no output. If I change the normalVector to:
normalVector (0 0.00001 1);
it succeeds.
Does anyone have any idea what is going on here?
(OpenFOAM 1.5-dev)

The 1.6.x version has some code to handle cases where the cuttingPlane aligned fairly closely to a cell face (where the normal cell cutting algorithm would otherwise fail horribly) -- this is what the zeroish and positive local variables are about.
I don't know if the sourceforge svn version has this, do you experience the same problem in 1.6.x?

[marupio]

That's the problem! Thanks Mark. I tried it in 1.6 and it works fine. So when I'm using 1.5-dev I'll just misalign the plane ever so slightly.

-Dave