@@ -1381,7 +1381,6 @@ contains
13811381 end do
13821382 end if
13831383
1384-
13851384 E_L = gamma_L*pres_L + pi_inf_L + 5.e-1_wp*rho_L*vel_L_rms + qv_L
13861385 E_R = gamma_R*pres_R + pi_inf_R + 5.e-1_wp*rho_R*vel_R_rms + qv_R
13871386
@@ -1706,11 +1705,10 @@ contains
17061705
17071706 ! Initialize all variables
17081707 vel_L_rms = 0._wp; vel_R_rms = 0._wp
1709- rho_L = 0._wp; rho_R = 0._wp;
1710- gamma_L = 0._wp; gamma_R = 0._wp;
1711- pi_inf_L = 0._wp; pi_inf_R = 0._wp;
1712- qv_L = 0._wp; qv_R = 0._wp;
1713-
1708+ rho_L = 0._wp; rho_R = 0._wp;
1709+ gamma_L = 0._wp; gamma_R = 0._wp;
1710+ pi_inf_L = 0._wp; pi_inf_R = 0._wp;
1711+ qv_L = 0._wp; qv_R = 0._wp;
17141712 !$acc loop seq
17151713 do i = 1, contxe
17161714 alpha_rho_L(i) = qL_prim_rs${XYZ}$_vf(j, k, l, i)
@@ -1738,7 +1736,7 @@ contains
17381736 qv_L = qv_L + alpha_rho_L(i)*qvs(i)
17391737 qv_R = qv_R + alpha_rho_R(i)*qvs(i)
17401738 end do
1741-
1739+
17421740 pres_L = qL_prim_rs${XYZ}$_vf(j, k, l, E_idx)
17431741 pres_R = qR_prim_rs${XYZ}$_vf(j + 1, k, l, E_idx)
17441742
@@ -1837,8 +1835,8 @@ contains
18371835 if (bubbles_euler) then
18381836 ! Put p_tilde in
18391837 flux_rs${XYZ}$_vf(j, k, l, contxe + dir_idx(i)) = &
1840- flux_rs${XYZ}$_vf(j, k, l, contxe + dir_idx(i)) - &
1841- dir_flg(dir_idx(i))*(xi_M*ptilde_L + xi_P*ptilde_R)
1838+ flux_rs${XYZ}$_vf(j, k, l, contxe + dir_idx(i)) - &
1839+ dir_flg(dir_idx(i))*(xi_M*ptilde_L + xi_P*ptilde_R)
18421840 end if
18431841 end do
18441842
@@ -1937,11 +1935,10 @@ contains
19371935
19381936 ! Initialize all variables
19391937 vel_L_rms = 0._wp; vel_R_rms = 0._wp
1940- rho_L = 0._wp; rho_R = 0._wp;
1941- gamma_L = 0._wp; gamma_R = 0._wp;
1942- pi_inf_L = 0._wp; pi_inf_R = 0._wp;
1943- qv_L = 0._wp; qv_R = 0._wp;
1944-
1938+ rho_L = 0._wp; rho_R = 0._wp;
1939+ gamma_L = 0._wp; gamma_R = 0._wp;
1940+ pi_inf_L = 0._wp; pi_inf_R = 0._wp;
1941+ qv_L = 0._wp; qv_R = 0._wp;
19451942 !$acc loop seq
19461943 do i = 1, num_fluids
19471944 alpha_L(i) = qL_prim_rs${XYZ}$_vf(j, k, l, E_idx + i)
0 commit comments