[DeanM]

Dear Foamers,

I would like to read a time-dependent source term, which can be defined by a function in a file in the case directory and read from the solver itself.
So far I managed to set the source term using funkySetFields, but only as initial condition. My approach is the following: In createFields.H the field of the mass source is defined as follows:

       
        volScalarField massSource
        (
        IOobject
        (
            "massSource",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ_IF_MODIFIED,
            IOobject::NO_WRITE
        ),
        mesh
        );

The equation to solve looks like

    fvScalarMatrix rhoEqn
    (
      fvm::ddt(rho) 
    + some terms
    == massSource
    );
    rhoEqn.solve();

Finally the user should be able to define "massSource" with a mathematical expression as function of time and space in a file in the case directory. I tried to do it with funkySetFields, it works well for time 0, but I need it for all times. Example:

    massSource
    {
        field massS;
        expression    "sin(pi*time())*(pos().x+pos().y"; //arbitrary function 
        //condition    "...";
    }

In conclusion: How can I define and read time-dependent source terms from the case directory?
Thank you for any ideas.
D.