[hxaxtma]

Hi guys,
in my code I have a code snippet which causes problems, maybe you can help me out with it. The code compiples, when I run a testcase the following error occurs I am using OF 4.1)

	U=rAU*UEqn.H();
		phi = (fvc::interpolate(U) & mesh.Sf())
                + fvc::ddtCorr(rAU, U, phi); -> This line causes the problem

This is the error message:

Time = 0.05

Courant Number mean: 0.595945 max: 14.185
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 9.7858e-06, No Iterations 140
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 9.63976e-06, No Iterations 148


--> FOAM FATAL ERROR:
dimensions of phi are not correct

From function Foam::tmp<typename Foam::fv::ddtScheme<Type>::fluxFieldType> Foam::fv::EulerDdtScheme<Type>::fvcDdtPhiCorr(cons t volScalarField&, const Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>&, const fluxFieldType&) [with Type = Foam::Vector<double>; typename Foam::fv::ddtScheme<Type>::fluxFieldType = Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>; Foam::volScalarField = Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>; Foam::fv::EulerDdtScheme<Type>::fluxFieldType = Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>]
in file finiteVolume/ddtSchemes/EulerDdtScheme/EulerDdtScheme.C at line 651.

Do you have an idea what I do wrong here?