[tests][mpnc] add molardensity in fluidstate for equilibrium fluidstate calculation

test/porousmediumflow/mpnc/implicit/combustionproblem1c.hh

@@ -333,6 +333,10 @@ public:
                 dummyCache,
                 wPhaseIdx);
         fluidState.setDensity(wPhaseIdx, density);
+        const Scalar molarDensity = FluidSystem::molarDensity(fluidState,
+                dummyCache,
+                wPhaseIdx);
+        fluidState.setMolarDensity(wPhaseIdx, molarDensity);
 
         for(int phaseIdx=0; phaseIdx<numPhases; phaseIdx++) {
             const Scalar h = FluidSystem::enthalpy(fluidState,

test/porousmediumflow/mpnc/implicit/evaporationatmosphereproblem.hh

@@ -333,9 +333,13 @@ public:
 
         // compute density of injection phase
         const Scalar density = FluidSystem::density(fluidState,
-                                                     dummyCache,
-                                                     gasPhaseIdx);
+                                                    dummyCache,
+                                                    gasPhaseIdx);
         fluidState.setDensity(gasPhaseIdx, density);
+        const Scalar molarDensity = FluidSystem::molarDensity(fluidState,
+                                                              dummyCache,
+                                                              gasPhaseIdx);
+        fluidState.setMolarDensity(gasPhaseIdx, molarDensity);
 
         for(int phaseIdx=0; phaseIdx<numPhases; phaseIdx++)
         {