Skip to content

Commit 31fd75c

Browse files
mohdsaid497566Malmahrouqi3
mohdsaid497566
authored and
Malmahrouqi3
committed
removed dir_idx indexing from the definition of wave speed subroutine
1 parent 53caf2e commit 31fd75c

File tree

1 file changed

+13
-13
lines changed

1 file changed

+13
-13
lines changed

‎src/simulation/m_riemann_solvers.fpp‎

Lines changed: 13 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -668,9 +668,9 @@ contains
668668
end do
669669
end if
670670

671-
call s_compute_wave_speed(wave_speeds, vel_L, vel_R, pres_L, pres_R, rho_L, rho_R, rho_avg, &
671+
call s_compute_wave_speed(wave_speeds, vel_L(dir_idx(1)), vel_R(dir_idx(1)), pres_L, pres_R, rho_L, rho_R, rho_avg, &
672672
c_L, c_R, c_avg, c_fast%L, c_fast%R, G_L, G_R, &
673-
tau_e_L, tau_e_R, &
673+
tau_e_L(dir_idx_tau(1)), tau_e_R(dir_idx_tau(1)), &
674674
s_L, s_R, s_S)
675675

676676
xi_M = (5e-1_wp + sign(5e-1_wp, s_L)) &
@@ -1402,9 +1402,9 @@ contains
14021402
end if
14031403

14041404
! COMPUTING THE DIRECT WAVE SPEEDS
1405-
call s_compute_wave_speed(wave_speeds, vel_L, vel_R, pres_L, pres_R, rho_L, rho_R, rho_avg, &
1405+
call s_compute_wave_speed(wave_speeds, vel_L(dir_idx(1)), vel_R(dir_idx(1)), pres_L, pres_R, rho_L, rho_R, rho_avg, &
14061406
c_L, c_R, c_avg, c_fast%L, c_fast%R, G_L, G_R, &
1407-
tau_e_L, tau_e_R, &
1407+
tau_e_L(dir_idx_tau(1)), tau_e_R(dir_idx_tau(1)), &
14081408
s_L, s_R, s_S)
14091409

14101410
! goes with q_star_L/R = xi_L/R * (variable)
@@ -1637,10 +1637,10 @@ contains
16371637
call s_compute_speed_of_sound(pres_R, rho_avg, gamma_avg, pi_inf_R, H_avg, alpha_R, &
16381638
vel_avg_rms, 0._wp, c_avg)
16391639

1640-
call s_compute_wave_speed(wave_speeds, vel_L, vel_R, pres_L, pres_R, rho_L, rho_R, rho_avg, &
1640+
call s_compute_wave_speed(wave_speeds, vel_L(dir_idx(1)), vel_R(dir_idx(1)), pres_L, pres_R, rho_L, rho_R, rho_avg, &
16411641
c_L, c_R, c_avg, c_fast%L, c_fast%R, G_L, G_R, &
1642-
tau_e_L, tau_e_R, &
1643-
s_L, s_R, s_S, s_M, s_P)
1642+
tau_e_L(dir_idx_tau(1)), tau_e_R(dir_idx_tau(1)), &
1643+
s_L, s_R, s_S)
16441644

16451645
! goes with q_star_L/R = xi_L/R * (variable)
16461646
! xi_L/R = ( ( s_L/R - u_L/R )/(s_L/R - s_star) )
@@ -1963,10 +1963,10 @@ contains
19631963
@:compute_low_Mach_correction()
19641964
end if
19651965

1966-
call s_compute_wave_speed(wave_speeds, vel_L, vel_R, pres_L, pres_R, rho_L, rho_R, rho_avg, &
1966+
call s_compute_wave_speed(wave_speeds, vel_L(dir_idx(1)), vel_R(dir_idx(1)), pres_L, pres_R, rho_L, rho_R, rho_avg, &
19671967
c_L, c_R, c_avg, c_fast%L, c_fast%R, G_L, G_R, &
1968-
tau_e_L, tau_e_R, &
1969-
s_L, s_R, s_S, s_M, s_P)
1968+
tau_e_L(dir_idx_tau(1)), tau_e_R(dir_idx_tau(1)), &
1969+
s_L, s_R, s_S)
19701970

19711971
! goes with q_star_L/R = xi_L/R * (variable)
19721972
! xi_L/R = ( ( s_L/R - u_L/R )/(s_L/R - s_star) )
@@ -2329,10 +2329,10 @@ contains
23292329
@:compute_low_Mach_correction()
23302330
end if
23312331

2332-
call s_compute_wave_speed(wave_speeds, vel_L, vel_R, pres_L, pres_R, rho_L, rho_R, rho_avg, &
2332+
call s_compute_wave_speed(wave_speeds, vel_L(dir_idx(1)), vel_R(dir_idx(1)), pres_L, pres_R, rho_L, rho_R, rho_avg, &
23332333
c_L, c_R, c_avg, c_fast%L, c_fast%R, G_L, G_R, &
2334-
tau_e_L, tau_e_R, &
2335-
s_L, s_R, s_S, s_M, s_P)
2334+
tau_e_L(dir_idx_tau(1)), tau_e_R(dir_idx_tau(1)), &
2335+
s_L, s_R, s_S)
23362336

23372337
! goes with q_star_L/R = xi_L/R * (variable)
23382338
! xi_L/R = ( ( s_L/R - u_L/R )/(s_L/R - s_star) )

0 commit comments

Comments
 (0)