Bug in InterMixingFoam

There is a bug in the code of the InterMixingFoam solver in alphaEqns. The problem is in the ordering. Ordering in current code: Step 1: phase 1-the convection flux is created, bounded, corrected and limited. Step 2: phase 2-the convection flux is created, bounded, corrected and limited. Step 3: phase 1-limited flux constructed. Step 4: phase 2-limited flux constructed. The problem is the coefficients calculated in mules limiter for phase 1 are being overwritten by coefficients calculated for phase 2 and the the SAME coefficients are being used to construct both phases.Order should be: step 1,3,2,4. Line 175 should be put to 110.

Current solver: LEFT picture. Non sharp water-air interface between, due to wrong mules coefficients being used.

Corrected solver: RIGHT picture

(diffusion is set between other and water phase)