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

expressions
 (
	sedimentvelocity
	{
	field Ua; //field to initialise
	fields
	(
	z Ua
	);
	lookuptables (
	{
	name inletprofilesedvel;
	outOfBounds clamp;
	file "$FOAM_CASE/1d_profil/Ua0.xy";
	}
	{
	name inletprofilesedvelv;
	outOfBounds clamp;
	file "$FOAM_CASE/1d_profil/Ua2.xy";
	}
	{
	name inletprofilesedvelw;
	outOfBounds clamp;
	file "$FOAM_CASE/1d_profil/Ua1.xy";
	}
	);
        value uniform (1 0 0);
	expression "vector(inletprofilesedvel(pos().z),inletprofilesedvelv(pos().z),inletprofilesedvelw(pos().z))";
	keepPatches 1;
}
	fluidvelocity
	{
	field Ub; //field to initialise
	fields
	(
	z Ub
	);
	lookuptables (
	{
	name inletprofilefvel;
	outOfBounds clamp;
	file "$FOAM_CASE/1d_profil/Ub0.xy";
	}
	{
	name inletprofilefvelv;
	outOfBounds clamp;
	file "$FOAM_CASE/1d_profil/Ub2.xy";
	}
	{
	name inletprofilefvelw;
	outOfBounds clamp;
	file "$FOAM_CASE/1d_profil/Ub1.xy";
	}
	);
        value uniform (1 0 0);
	expression "vector(inletprofilefvel(pos().z),inletprofilefvelv(pos().z),inletprofilefvelw(pos().z))";
	keepPatches 1;
}

       alpha
        {
        field alpha_a; //field to initialise
        fields
        (
        z alpha
        );
        lookuptables (
        {
        name inletprofilealpha;
        outOfBounds clamp;
        file "$FOAM_CASE/1d_profil/alpha_a.xy";
        }
        );
        expression "inletprofilealpha(pos().z)";
        keepPatches 1;
        }
      pb
        {
        field p_rbgh; //field to initialise
        fields
        (
        z p_rbgh
        );
        lookuptables (
        {
        name inletprofilepb;
        outOfBounds clamp;
        file "$FOAM_CASE/1d_profil/p_rbgh.xy";
        }
        );

        expression "inletprofilepb(pos().z)-3.998*(pos().x-0.8)";
        keepPatches 1;
        }

       k
        {
        field k; //field to initialise
        fields
        (
        z k
        );
        lookuptables (
        {
        name inletprofilek;
        outOfBounds clamp;
        file "$FOAM_CASE/1d_profil/k.xy";
        }
        );
        expression "inletprofilek(pos().z)";
        keepPatches 1;
        }
       omega
        {
        field omega; //field to initialise
        fields
        (
        z omega
        );
        lookuptables (
        {
        name inletprofileomega;
        outOfBounds clamp;
        file "$FOAM_CASE/1d_profil/omega.xy";
        }
        );
        expression "inletprofileomega(pos().z)";
        keepPatches 1;
        }
 );

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