[houss]

Hi Foamers!
i am simulating the flow in a axial diffuser with rhosimplefoam. i definde a coded functionobject to get values like (ptot,rho,u,p,pressure loss between inlet and outlet) on certain patches (INLET AND OUTLET). Now i am able to display these values through the solver running but i want actually to save these in a separate file for each outputinterval..briefly it should look like this:
ui pi ri ptoti uo po ro ptoto pl
x x x x x x x x
...
..

controldict

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

application rhoSimpleFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 6;

deltaT 1;

writeControl timeStep;

writeInterval 3;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression off;

timeFormat general;

timePrecision 6;

graphFormat raw;

runTimeModifiable true;

functions
{
#include "pressure_loss"
}






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

pressure_loss

pressure_loss
{
type coded;
functionObjectLibs ("libutilityFunctionObjects.so");
redirectType error;
enabled true;
outputControl timeStep;
outputInterval 1;
fields
(
phi
rho
U
p
);
code
#{
//LABEL PATCHES
label IN = mesh().boundaryMesh().findPatchID("INLET");
label OUT = mesh().boundaryMesh().findPatchID("OUTLET");

//const fvPatchVectorField& in_u = U.boundaryField()[IN];
if (IN == -1)
{
Info <<"patch named INLET DOES NOT EXIST!"<< endl;
}
if (OUT == -1)
{
Info <<"patch named OUTLET DOES NOT EXIST!"<< endl;
}
const volVectorField& U = mesh().lookupObject<volVectorField>("U");
scalar u_i = 0.0;
u_i = mag(average(U.boundaryField()[IN]));
Info<<"AVERAGED VELOCITY INLET =\n"<<u_i<< endl;

const volScalarField& p = mesh().lookupObject<volScalarField>("p");
scalar p_i = 0.0;
p_i = average(p.boundaryField()[IN]);
Info<<"AVERAGED STATIC PRESSURE INLET =\n"<<p_i<< endl;

const volScalarField& rho = mesh().lookupObject<volScalarField>("rho");
scalar r_i = 0.0;
r_i = average(rho.boundaryField()[IN]);
Info<<"AVERAGED DENSITY INLET =\n"<<r_i<< endl;

const surfaceScalarField& phi = mesh().lookupObject<surfaceScalarField>("phi");
scalar m_i = 0.0;
m_i = sum(phi.boundaryField()[IN]);
Info<<"MASSFLOW INLET =\n" <<m_i<<"\n\n" <<endl;

scalar u_o = 0.0;
u_o = mag(average(U.boundaryField()[OUT]));
Info<<"AVERAGED VELOCITY OUTLET =\n"<<u_o<< endl;
scalar p_o = 0.0;
p_o = average(p.boundaryField()[OUT]);
Info<<"AVERAGED STATIC PRESSURE OUTLET =\n"<<p_o<< endl;
scalar r_o = 0.0;
r_o = average(rho.boundaryField()[OUT]);
Info<<"AVERAGED DENSITY OUTLET =\n"<<r_o<< endl;

scalar m_o = 0.0;
m_o = sum(phi.boundaryField()[OUT]);
Info<<"MASSFLOW OUTLET =\n" <<m_o<<"\n\n" <<endl;

//CALCULATE TOTAL PRESSURE DROP
scalar p_l = 0.0;
p_l = p_i + r_i*pow(u_i,2)/2 -(p_o + r_o*pow(u_o,2)/2);
Info<<"PRESSURE DROP =\n"<<p_l<<endl;

#};

}

i would really appreciate your help.