Skip to content

Commit 1bc7a1d

Browse files
mrodrig6mrodrig6
committed
formatted wp in riemann solver
1 parent 719a164 commit 1bc7a1d

File tree

1 file changed

+8
-9
lines changed

1 file changed

+8
-9
lines changed

‎src/simulation/m_riemann_solvers.fpp‎

Lines changed: 8 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -522,13 +522,13 @@ contains
522522
call get_mixture_energy_mass(T_L, Ys_L, E_L)
523523
call get_mixture_energy_mass(T_R, Ys_R, E_R)
524524

525-
E_L = rho_L*E_L + 0.5*rho_L*vel_L_rms
526-
E_R = rho_R*E_R + 0.5*rho_R*vel_R_rms
525+
E_L = rho_L*E_L + 0.5_wp*rho_L*vel_L_rms
526+
E_R = rho_R*E_R + 0.5_wp*rho_R*vel_R_rms
527527
H_L = (E_L + pres_L)/rho_L
528528
H_R = (E_R + pres_R)/rho_R
529529
else
530-
E_L = gamma_L*pres_L + pi_inf_L + 0.5*rho_L*vel_L_rms + qv_L
531-
E_R = gamma_R*pres_R + pi_inf_R + 0.5*rho_R*vel_R_rms + qv_R
530+
E_L = gamma_L*pres_L + pi_inf_L + 0.5_wp*rho_L*vel_L_rms + qv_L
531+
E_R = gamma_R*pres_R + pi_inf_R + 0.5_wp*rho_R*vel_R_rms + qv_R
532532
H_L = (E_L + pres_L)/rho_L
533533
H_R = (E_R + pres_R)/rho_R
534534
end if
@@ -2359,14 +2359,13 @@ contains
23592359
call get_mixture_energy_mass(T_L, Ys_L, E_L)
23602360
call get_mixture_energy_mass(T_R, Ys_R, E_R)
23612361
2362-
E_L = rho_L*E_L + 0.5*rho_L*vel_L_rms
2363-
E_R = rho_R*E_R + 0.5*rho_R*vel_R_rms
2362+
E_L = rho_L*E_L + 0.5_wp*rho_L*vel_L_rms
2363+
E_R = rho_R*E_R + 0.5_wp*rho_R*vel_R_rms
23642364
H_L = (E_L + pres_L)/rho_L
23652365
H_R = (E_R + pres_R)/rho_R
23662366
else
2367-
E_L = gamma_L*pres_L + pi_inf_L + 0.5*rho_L*vel_L_rms + qv_L
2368-
2369-
E_R = gamma_R*pres_R + pi_inf_R + 0.5*rho_R*vel_R_rms + qv_R
2367+
E_L = gamma_L*pres_L + pi_inf_L + 0.5_wp*rho_L*vel_L_rms + qv_L
2368+
E_R = gamma_R*pres_R + pi_inf_R + 0.5_wp*rho_R*vel_R_rms + qv_R
23702369
23712370
H_L = (E_L + pres_L)/rho_L
23722371
H_R = (E_R + pres_R)/rho_R

0 commit comments

Comments
 (0)