/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  9
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       volScalarField;
    object      scalar1;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            codedFixedValue;
        value           uniform 1;
        
        name halfBore;
        code
        #{
        	const vectorField& Cf = patch().Cf();

		scalarField& field = *this;

		forAll(Cf, faceI)
		{
			const scalar r = Cf[faceI][1];
			if (r>0)
			{
			field[faceI] = 1;
			}
			
			else
			{
			field[faceI] = 0;
			}
		}
        #};
        
    }

    outlet
    {
        type            zeroGradient;
    }

    wall
    {
        type            zeroGradient;
    }
}

// ************************************************************************* //
