[be_inspired]

Hi all,

Using OF-v1612+:
For an axi-symmetric simulation, I want to calculate the distance from every cell to the axis of rotation. After that, I want to generate an iso-surface of radius equal to 50 centimeters, and over that isosurface, then I could have the velocity field that is what I want lastly (using sampling).

The piece of code to calculate the cell coordinates would be something like this but nothing is written:


functions
{
radius
{
type coded;
functionObjectLibs ("libutilityFunctionObjects.so");
name radius;
writeControl writeTime;
redirectType radius;
codeExecute
#{
// Lookup U

const volVectorField& C = mesh().C();
//const volScalarField zcoord = mesh.C().component(vector::Z);
//const volScalarField ycoord = mesh.C().component(vector::Y);
//volScalarField ycoord = ;
//zcoord = mesh.C().component(vector::Z);
//const volScalarField& radius = ((y**2+z**2)**0.5);
// Write
C.write();
Info<< endl << "hello test"<<endl;
#};

}

Any suggestion would be great.
Best Regards