[crubio.abujas]

I'm sorry but I think we are talking about different versions of the code. What loop are you talking about? In the final version there's no loop at all! As we've been back and from with the code a lot this time maybe there could be some misunderstanding.

I've tried this version of the code and worked for me. I've added a bunch of CUSTOM: flags to have some insight of the code behavior.

    inlet
    {
        type            codedFixedValue;
        value           uniform 0;
        name            CardiacImpuse;
        code
        #{
            const scalar &t = this ->db().time().value();
            const scalar dt = this ->db().time().deltaTValue();

            const surfaceScalarField& phi = db().lookupObject<surfaceScalarField>("phi"); 
            const fvsPatchField<scalar>& phip =
                patch().patchField<surfaceScalarField,scalar> (phi);
            const fvsPatchField<scalar>& phip_last =
                patch().patchField<surfaceScalarField,scalar> (phi.oldTime());

            // Constants
            const scalar C=1.16e-05;
            const scalar R1=1.64e+03;
            const scalar R2=1.79e+05;
            const scalar p_0=7.55;
            const scalar Q_0=1.03e-4;
            const scalar T_REF=0.1;

            // Static variables
            static scalar integrand (0);
            static scalar current_time (0);

            scalar tau = std::fmod(t, T_REF);

             // Check if the timestep in new
            if (current_time < t) {

                // Update the stored timestep
                current_time = t;

                Info << "CUSTOM: current_time=" << current_time << nl
                     << "CUSTOM: dt=" << dt << nl
                     << "CUSTOM: tau=" << tau << endl;

                // Update the integrand with the previous value
                integrand += (exp(tau/(R2*C))*gSum(phip_last)*dt);
                Info << "CUSTOM: Updated integrad to " << integrand << endl;

                // Check for reset condition
                if (tau < dt) {
                    integrand = 0;
                    Info << "CUSTOM: Reseting the integrand" << endl;
                }
            }

            operator ==(R1*gSum(phip)+exp(-tau/(R2*C))*((p_0-R1*Q_0)+(1/C)*integrand));
        #};
    }

[Neb]

Oh No, I'm sorry, my mistake, I meant in the IF, not loop!: Eek:
The simulation has now ended and it worked!!!
I still have to fix something, but finally the pressure values ​​are correct!

I don't know how to thank you, you were very kind.

[Neb]

I am very sorry to to disturb again.

Unfortunately there is one last problem.
I checked the values, the code compiles without errors and the pressures are now correct, but it does not read the IF. So it do not apply the correction to t = 1 with integrand=0.

I have tried other order:

if (tau<dt)
{
operator ==(7.55);
integrand=0;
}
and this too:

const volScalarField & p = db (). lookupObject <volScalarField> ("p");
const fvPatchField<scalar> & pip = patch().patchField<volScalarField,scalar> (p);
fvPatchField& P_new=const_cast<fvPatchScalarField&>(pip);
if ( tau<dt)
{
forAll (P_new , facei)
{ P_new[facei]=7.55; } }

but it continues the calculation ignoring them, the final pressure always remains different from 7.55 and it's always the same value. How is it possible?
I want that at time = 1 the pressure is 7.55.

[Neb]

I solved the code. : D: D: D

Thank you very much for all the suggestions have been very useful!