[mkhm]

Dear foamers,


I am trying to create a new post-processing tool. To do that, I followed this tutorial: http://www.tfd.chalmers.se/~hani/kur...CFD_Report.pdf . I can compile successfully but nothing is written in my time folders. I tried to copy one of the case in OpenFOAM-4.x/src/functionObjects/field/, give it another name, change the name in .C and .H file and mention this name in "file" located in OpenFOAM-4.x/src/functionObjects/field/Make/. No success. There is no compilation error and when I run rhoCentralFoam -postProcess -func newFunction as a command line, the necessary fields are read but nothing is written. To be more precise, when we execute for instance MachNo, we have :


Reading field U

Creating turbulence model

Selecting turbulence model type laminar
functionObjects::MachNo MachNo writing field: Ma


In my case withe newFunction, the line "functionObjects::MachNo MachNo writing field: Ma" (with Ma or MachNo = new names) is missing.



I tried also to include the following line in my controlDic:


functions
{
k
{
// Type of functionObject
type k;

// Where to load it from (if not already in solver)
libs ("libfieldFunctionObjects.so");

// Function object enabled flag
enabled true;

// When to output the average fields
writeControl writeTime;
}
}


But neither, nothing is written in my time folders.



I know that by modifying the solver, we can make it write the same variables that we write by executing postProcess commands, but I want to learn where in my procedure of creating new functionObjects, I am wrong.



I appreciate your help dear foamers,



regards

[mkhm]

No one could help me ?

[Alishaha2]

you need to add following lines to your .C - main function

Foam::timeSelector::addOptions();
#include "addRegionOption.H"
#include "addDictOption.H"

#include "setRootCase.H"
#include "createTime.H"
Foam::instantList timeDirs = Foam::timeSelector::select0(runTime, args);
#include "createNamedMesh.H"

forAll(timeDirs, timeI)
{
runTime.setTime(timeDirs[timeI], timeI);

Foam::Info<< "Time = " << runTime.timeName() << Foam::endl;

mesh.readUpdate();

Foam::Yourfunction(args, runTime, mesh);

Foam::Info<< Foam::endl;