/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2112                                  |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    arch        "LSB;label=32;scalar=64";
    class       volVectorField;
    location    "17";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);
boundaryField
{
    inlet
    {
        type                turbulentDFSEMInlet;
        delta               1;  //Characteristic length scale
        U                   uniform (0 0 1); //mean velocity
        R                   uniform (0.2 0 0 0.2 0 0.2); // Reynolds stress: <Rxx> <Rxy> <Rxz> <Ryy> <Ryz> <Rzz>
        L                   uniform 0.4; //Integral length scale
        nCellPerEddy        1;  //Minimum eddy length in units of number of cells
        value               uniform (0 0 1);
	}
	//~ inlet
	//~ {
		//~ type            codedFixedValue;
		//~ value           uniform (0 0 0);
		//~ name            myInlet;
		//~ code
		//~ #{
			//~ scalar U_max = 2;
			//~ const fvPatch& boundaryPatch = this->patch();
			//~ const vectorField& Cf = boundaryPatch.Cf();
			//~ vectorField& field = *this;
			//~ forAll(boundaryPatch, i)
			//~ {
				//~ scalar r = sqrt(Cf[i].y()*Cf[i].y() + Cf[i].x()*Cf[i].x())/0.5;
				//~ field[i] = vector(0, 0, U_max*Foam::pow(1.0-r, 1.0/7.0));
			//~ }
		//~ #};
	//~ }
    outlet
    {
        type            advective;
        value           uniform (0 0 0);
    }
    walls
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
}


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