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

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

internalField   uniform (0 0 0);

boundaryField {

  up
  {
    type            symmetryPlane;
  }
  down
  {
    type            symmetryPlane;
  }
  left
  {
    type            fixedValue;
    value           uniform (0 0 0);
  }
  right
  {
    type            fixedValue;
    value           uniform (0 0 0);
  }
  flow_to_porousMat
  {


    type            codedFixedValue;
    name   	        ablation;

    code
#{
    const vector axis(0, 0, 0);
    tmp<vectorField> nf_tmp =
    - this->patch().Sf()
    / this->patch().magSf();
    const vectorField& nf = nf_tmp();
    const fvPatchField<scalar>& YO2 =
    patch().lookupPatchField<volScalarField, scalar>("O2");
    operator==(-nf*YO2*10/2000);
#};
    value           $internalField;

  }
}



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