CFD Online Discussion Forums - application help

Dear users, I'm using openfoam 2.1 and I need to create a post-processing application to get this filed:

acoSor=rho*fvc::div(fvc::div(U*U))

This field can be used as a source term for a FW-H acoustic solver. I already tried to modify other utilities (for example the one for vorticity) and get it done, but I always get errors.

this is my code actually

Code:

#include "calc.H"

#include "fvc.H"

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

     

void Foam::calc(const argList& args, const Time& runTime, const fvMesh& mesh)

{

    bool writeResults = !args.optionFound("noWrite");



    IOobject Uheader

    (

        "U",

        runTime.timeName(),

        mesh,

        IOobject::MUST_READ

    );



    if (Uheader.headerOk())

    {

        Info<< "    Reading U" << endl;

        volVectorField U(Uheader, mesh);



        Info<< "    Calculating vorticity" << endl;

        volScalarField acoSor

        (

            IOobject

            (

                "acoSor",

                runTime.timeName(),

                mesh,

                IOobject::NO_READ

            ),

            rho*fvc::div(fvc::div(U*U))

        );

        Info<< "    Calculating acoSor" << endl;



        if (writeResults)

        {

            acoSor.write();

        }

    }

    else

    {

        Info<< "    No U" << endl;

    }



    Info<< "\nEnd\n" << endl;

}

Can you help me please? Thank you very much!