[Vitor Monteiro]

The code was exactly this one I was using later. But instead it was implemented at the multiphaseEulerFoam.C file and recompiled again.

forAll(fluid.multiComponentPhases(), multiComponentPhasei) // fases
      {
        phaseModel& phase = fluid.multiComponentPhases()[multiComponentPhasei];
        UPtrList<volScalarField>& Y = phase.YActiveRef();
        const volScalarField& alpha = phase;
        const volScalarField& rho = phase.rho();
    const scalarField& v = mesh().V();

        forAll(Y, i) // espécies
        {

     const double& m = gSum(alpha*rho*v);
     outputFilePtr() << mesh().time().timeName() << "    "
            << Y[i].name() << "    "
            << m  << endl;

         Info<< "coded totalMass output:" << endl
            << "    name: " << Y[i].name() << endl
        << "    mass: " << m << endl
            << endl;

     }; // fim for
       }; //fim for