add mpp_lim to mixture rule acc

Author: Hyeoksu Lee

src/common/m_variables_conversion.fpp

@@ -339,6 +339,13 @@ contains
             pi_inf_K = pi_infs(1)
             qv_K = qvs(1)
         else
+            if (mpp_lim) then
+                do i = 1, num_fluids
+                    alpha_rho_K(i) = max(0._wp, alpha_rho_K(i))
+                    alpha_K(i) = min(max(0._wp, alpha_K(i)), 1._wp)
+                end do
+                alpha_K = alpha_K/max(sum(alpha_K), sgm_eps)
+            end if
             rho_K = 0._wp; gamma_K = 0._wp; pi_inf_K = 0._wp; qv_K = 0._wp
             do i = 1, num_fluids
                 rho_K = rho_K + alpha_rho_K(i)
@@ -609,8 +616,6 @@ contains
 
                     call s_compute_species_fraction(qK_cons_vf, j, k, l, alpha_rho_K, alpha_K)
 
-                    ! print *, j, k, l, alpha_rho_K, qK_cons_vf(1)%sf(j, k, l), alpha_K, qK_cons_vf(advxb)%sf(j, k, l)
-
                     if (model_eqns /= 4) then
 #ifdef MFC_SIMULATION
                         ! If in simulation, use acc mixture subroutines
@@ -1340,7 +1345,6 @@ contains
                 alpha_rho_K(i) = max(0._wp, alpha_rho_K(i))
                 alpha_K(i) = min(max(0._wp, alpha_K(i)), 1._wp)
             end do
-
             alpha_K = alpha_K/max(sum(alpha_K), 1.e-16_wp)
         end if
 

src/simulation/m_bubbles_EL.fpp

@@ -642,11 +642,7 @@ contains
             call s_get_pinf(k, q_prim_vf, 1, myPinf, cell, aux1, aux2)
 
             ! Obtain liquid density and computing speed of sound from pinf
-            $:GPU_LOOP(parallelism='[seq]')
-            do i = 1, num_fluids
-                myalpha_rho(i) = q_prim_vf(i)%sf(cell(1), cell(2), cell(3))
-                myalpha(i) = q_prim_vf(E_idx + i)%sf(cell(1), cell(2), cell(3))
-            end do
+            call s_compute_species_fraction(q_prim_vf, cell(1), cell(2), cell(3), myalpha_rho, myalpha)
             call s_convert_species_to_mixture_variables_acc(myRho, gamma, pi_inf, qv, myalpha, &
                                                             myalpha_rho, Re)
             call s_compute_cson_from_pinf(q_prim_vf, myPinf, cell, myRho, gamma, pi_inf, myCson)

src/simulation/m_hyperelastic.fpp

@@ -110,11 +110,9 @@ contains
         do l = 0, p
             do k = 0, n
                 do j = 0, m
-                    $:GPU_LOOP(parallelism='[seq]')
-                    do i = 1, num_fluids
-                        alpha_rho_k(i) = q_cons_vf(i)%sf(j, k, l)
-                        alpha_k(i) = q_cons_vf(advxb + i - 1)%sf(j, k, l)
-                    end do
+
+                    call s_compute_species_fraction(q_cons_vf, j, k, l, alpha_rho_k, alpha_k)
+
                     ! If in simulation, use acc mixture subroutines
                     call s_convert_species_to_mixture_variables_acc(rho, gamma, pi_inf, qv, alpha_k, &
                                                                     alpha_rho_k, Re, G_local, Gs_hyper)