[Za-ck]

Hi,

I would like to plot the alpha field on a cutting plane using the runTimePostProcessing tool. Everything looks right without interpolation, but with activated interpolation the values seems to get mixed up. Same error on an isoSurface. No problems using a plane instead of a cuttingPlane.

Any ideas?

Cheers,
Za-ck

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1906                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
	version     1906;
	format      ascii;
	class       dictionary;
	location    "system";
	object      visualization;
}


functions
{
    planes
    {
        type            surfaces;
        libs            ("libsampling.so");

        writeControl    writeTime;
        writeFields     true;

        fields          (alpha.water);

        surfaceFormat   none;

        store           true;

        //cell cellPatchConstrained cellPoint cellPointFace cellPointWallModified pointMVC
        interpolationScheme cellPoint; 

        

        surfaces
        {
            xPlane
            {
                type            cuttingPlane; //plane; //
                planeType       pointAndNormal;
                interpolate     true;
//                 regularise      false;
                pointAndNormalDict
                {
                    point   (0 0 0);
                    normal  (1 0 0);
                }
             }            
                
         };

    }



    postPro1
    {
        #includeEtc "caseDicts/postProcessing/visualization/runTimePostPro.cfg"

        // Time control etc
        writeControl    writeTime;

        output
        {
            name        image;
            width       1920;
            height      1200;
        }
        
        
        
        camera
        {
            nFrameTotal     1;
            parallelProjection  no;
            position        (-20 0 0);
            focalPoint      (-1 0 0.0);

            up              (0.0 0.0 1.0);
        }

        // Default colours
        colours
        {
            ${..colourScheme.paraview};
        }

        surfaces
        {
            //A cutting plane from sampled surfaces:
            stored1
            {
                type            functionObjectSurface;
                functionObject  planes.xPlane; //planes.freeSurface; //planes.test2; //
                colourMap       coolToWarm;
                representation  surface;
                visible         yes;
                featureEdges    no;
                colourBy        field;
                field           alpha.water;
                range           (0.0 1.0);
                opacity         1;

            }
  

        }
        
        text
        {
            
        }

    }

}

[tecmul]

I'm confused about what the "interpolate" switch does. Doesn't sampling on a plane always need some form of interpolation?