Skip to content

Commit fbddfa9

Browse files
sbryngelsonsbryngelson
authored
Raise the amdflang species bound so mechanisms above ten species build (#1852)
1 parent 03d5bbf commit fbddfa9

12 files changed

Lines changed: 101 additions & 82 deletions

‎src/common/include/shared_parallel_macros.fpp‎

Lines changed: 7 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -8,6 +8,13 @@
88
#:set USING_CCE = (MFC_COMPILER == CCE_COMPILER_ID)
99
#:set USING_AMD = (MFC_COMPILER == AMD_COMPILER_ID)
1010

11+
#! Fallback extents the USING_AMD guards substitute for device-global array bounds when case
12+
#! optimization is off. They are not independent: sys_size counts the species, so AMD_SYS_SIZE_MAX
13+
#! must cover 3*num_fluids + num_vels + 1 + AMD_NUM_SPECIES_MAX. Keep them here rather than as
14+
#! literals at each declaration, so raising one cannot silently outgrow the other.
15+
#:set AMD_NUM_SPECIES_MAX = 60
16+
#:set AMD_SYS_SIZE_MAX = 70
17+
1118
#:def ASSERT_LIST(data, datatype)
1219
#:assert data is not None
1320
#:assert isinstance(data, list)

‎src/common/m_checker_common.fpp‎

Lines changed: 10 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -54,7 +54,16 @@ contains
5454
#:if not MFC_CASE_OPTIMIZATION
5555
@:PROHIBIT(num_fluids > 3, "num_fluids <= 3 for AMDFLang when Case optimization is off")
5656
@:PROHIBIT((bubbles_euler .or. bubbles_lagrange) .and. nb > 3, "nb <= 3 for AMDFLang when Case optimization is off")
57-
@:PROHIBIT(chemistry .and. num_species > 10, "num_species > 10 for AMDFLang when Case optimization is off")
57+
@:PROHIBIT(chemistry .and. num_species > ${AMD_NUM_SPECIES_MAX}$, &
58+
& "num_species <= ${AMD_NUM_SPECIES_MAX}$ for AMDFLang when Case optimization is off")
59+
! The sys_size bound is not independent of the one above it. Chemistry pins num_fluids to 1, so with
60+
! num_vels <= 3 the species terminate sys_size at 2*1 + 3 + 1 + 60 = 66 for five equations and
61+
! 3*1 + 3 + 1 + 60 = 67 for six; 70 covers both with room, and hypoelastic stresses would add up to six
62+
! more. It had no check of its own while the species cap was ten, because sys_size could not then reach
63+
! the dimension(20) the guard gives every sys_size array; HLLC's star states have no other bound, so
64+
! raising one cap without the other overruns them with nothing to say so.
65+
@:PROHIBIT(sys_size > ${AMD_SYS_SIZE_MAX}$, &
66+
& "sys_size <= ${AMD_SYS_SIZE_MAX}$ for AMDFLang when Case optimization is off")
5867
#:endif
5968
6069
end subroutine s_check_amd

‎src/common/m_chemistry.fpp‎

Lines changed: 7 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -116,8 +116,8 @@ contains
116116
real(wp) :: rho, omega_m
117117

118118
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
119-
real(wp), dimension(10) :: Ys
120-
real(wp), dimension(10) :: omega
119+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: Ys
120+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: omega
121121
#:else
122122
real(wp), dimension(num_species) :: Ys
123123
real(wp), dimension(num_species) :: omega
@@ -175,7 +175,7 @@ contains
175175
real(wp), parameter :: stiff_target = 0.5_wp
176176
177177
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
178-
real(wp), dimension(10) :: Ys, cdot, ddot, y0, prod0, Lloss, alp
178+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: Ys, cdot, ddot, y0, prod0, Lloss, alp
179179
#:else
180180
real(wp), dimension(num_species) :: Ys, cdot, ddot, y0, prod0, Lloss, alp
181181
#:endif
@@ -327,10 +327,10 @@ contains
327327
type(scalar_field), intent(in) :: q_T_sf
328328
329329
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
330-
real(wp), dimension(10) :: Xs_L, Xs_R, Xs_cell, Ys_L, Ys_R, Ys_cell
331-
real(wp), dimension(10) :: mass_diffusivities_mixavg1, mass_diffusivities_mixavg2
332-
real(wp), dimension(10) :: mass_diffusivities_mixavg_Cell, dXk_dxi, h_l, h_r, h_k
333-
real(wp), dimension(10) :: Mass_Diffu_Flux, dYk_dxi
330+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: Xs_L, Xs_R, Xs_cell, Ys_L, Ys_R, Ys_cell
331+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: mass_diffusivities_mixavg1, mass_diffusivities_mixavg2
332+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: mass_diffusivities_mixavg_Cell, dXk_dxi, h_l, h_r, h_k
333+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: Mass_Diffu_Flux, dYk_dxi
334334
#:else
335335
real(wp), dimension(num_species) :: Xs_L, Xs_R, Xs_cell, Ys_L, Ys_R, Ys_cell
336336
real(wp), dimension(num_species) :: mass_diffusivities_mixavg1, mass_diffusivities_mixavg2

‎src/common/m_variables_conversion.fpp‎

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -468,7 +468,7 @@ contains
468468
#:if USING_AMD and not MFC_CASE_OPTIMIZATION
469469
real(wp), dimension(3) :: alpha_K, alpha_rho_K
470470
real(wp), dimension(3) :: nRtmp
471-
real(wp) :: rhoYks(1:10)
471+
real(wp) :: rhoYks(1:${AMD_NUM_SPECIES_MAX}$)
472472
#:else
473473
real(wp), dimension(num_fluids) :: alpha_K, alpha_rho_K
474474
real(wp), dimension(nb) :: nRtmp
@@ -1021,10 +1021,10 @@ contains
10211021
! functions, the shear and volume Reynolds numbers and the Weber numbers
10221022
10231023
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
1024-
real(wp), dimension(3) :: alpha_rho_K
1025-
real(wp), dimension(3) :: alpha_K
1026-
real(wp), dimension(3) :: vel_K
1027-
real(wp), dimension(10) :: Y_K
1024+
real(wp), dimension(3) :: alpha_rho_K
1025+
real(wp), dimension(3) :: alpha_K
1026+
real(wp), dimension(3) :: vel_K
1027+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: Y_K
10281028
#:else
10291029
real(wp), dimension(num_fluids) :: alpha_rho_K
10301030
real(wp), dimension(num_fluids) :: alpha_K

‎src/simulation/m_cbc.fpp‎

Lines changed: 10 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -475,18 +475,18 @@ contains
475475
real(wp) :: dpres_ds
476476
477477
#:if USING_AMD
478-
real(wp), dimension(20) :: L
478+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$) :: L
479479
#:else
480480
real(wp), dimension(sys_size) :: L
481481
#:endif
482482
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
483-
real(wp), dimension(3) :: alpha_rho, dalpha_rho_ds, mf
484-
real(wp), dimension(3) :: vel, dvel_ds
485-
real(wp), dimension(3) :: adv_local, dadv_ds
486-
real(wp), dimension(3) :: dadv_dt
487-
real(wp), dimension(3) :: dvel_dt
488-
real(wp), dimension(3) :: dalpha_rho_dt
489-
real(wp), dimension(10) :: Ys, h_k, dYs_dt, dYs_ds, Xs, Gamma_i, Cp_i
483+
real(wp), dimension(3) :: alpha_rho, dalpha_rho_ds, mf
484+
real(wp), dimension(3) :: vel, dvel_ds
485+
real(wp), dimension(3) :: adv_local, dadv_ds
486+
real(wp), dimension(3) :: dadv_dt
487+
real(wp), dimension(3) :: dvel_dt
488+
real(wp), dimension(3) :: dalpha_rho_dt
489+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: Ys, h_k, dYs_dt, dYs_ds, Xs, Gamma_i, Cp_i
490490
#:else
491491
real(wp), dimension(num_fluids) :: alpha_rho, dalpha_rho_ds, mf
492492
real(wp), dimension(num_vels) :: vel, dvel_ds
@@ -646,8 +646,8 @@ contains
646646
!> gamma_method = 1: Ref. Section 2.3.1 Formulation of doi:10.7907/ZKW8-ES97.
647647
call get_mole_fractions(Mw, Ys, Xs)
648648
call get_species_specific_heats_r(T, Cp_i)
649-
Gamma_i = Cp_i/(Cp_i - 1.0_wp)
650-
gamma = sum(Xs(:)/(Gamma_i(:) - 1.0_wp))
649+
Gamma_i(1:num_species) = Cp_i(1:num_species)/(Cp_i(1:num_species) - 1.0_wp)
650+
gamma = sum(Xs(1:num_species)/(Gamma_i(1:num_species) - 1.0_wp))
651651
else if (chem_params%gamma_method == 2) then
652652
!> gamma_method = 2: c_p / c_v where c_p, c_v are specific heats.
653653
call get_mixture_specific_heat_cv_mass(T, Ys, Cv)

‎src/simulation/m_compute_cbc.fpp‎

Lines changed: 25 additions & 25 deletions
Original file line numberDiff line numberDiff line change
@@ -38,7 +38,7 @@ contains
3838

3939
$:GPU_ROUTINE(parallelism='[seq]')
4040
#:if USING_AMD
41-
real(wp), dimension(20), intent(inout) :: L
41+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$), intent(inout) :: L
4242
#:else
4343
real(wp), dimension(sys_size), intent(inout) :: L
4444
#:endif
@@ -62,7 +62,7 @@ contains
6262

6363
$:GPU_ROUTINE(parallelism='[seq]')
6464
#:if USING_AMD
65-
real(wp), dimension(20), intent(inout) :: L
65+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$), intent(inout) :: L
6666
#:else
6767
real(wp), dimension(sys_size), intent(inout) :: L
6868
#:endif
@@ -85,7 +85,7 @@ contains
8585

8686
$:GPU_ROUTINE(parallelism='[seq]')
8787
#:if USING_AMD
88-
real(wp), dimension(20), intent(inout) :: L
88+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$), intent(inout) :: L
8989
#:else
9090
real(wp), dimension(sys_size), intent(inout) :: L
9191
#:endif
@@ -108,12 +108,12 @@ contains
108108

109109
$:GPU_ROUTINE(parallelism='[seq]')
110110
#:if USING_AMD
111-
real(wp), dimension(20), intent(inout) :: L
111+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$), intent(inout) :: L
112112
#:else
113113
real(wp), dimension(sys_size), intent(inout) :: L
114114
#:endif
115115
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
116-
real(wp), dimension(10), intent(in) :: dYs_ds
116+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$), intent(in) :: dYs_ds
117117
#:else
118118
real(wp), dimension(num_species), intent(in) :: dYs_ds
119119
#:endif
@@ -135,7 +135,7 @@ contains
135135

136136
real(wp), dimension(3), intent(in) :: lambda
137137
#:if USING_AMD
138-
real(wp), dimension(20), intent(inout) :: L
138+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$), intent(inout) :: L
139139
#:else
140140
real(wp), dimension(sys_size), intent(inout) :: L
141141
#:endif
@@ -159,15 +159,15 @@ contains
159159

160160
real(wp), dimension(3), intent(in) :: lambda
161161
#:if USING_AMD
162-
real(wp), dimension(20), intent(inout) :: L
162+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$), intent(inout) :: L
163163
#:else
164164
real(wp), dimension(sys_size), intent(inout) :: L
165165
#:endif
166166
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
167-
real(wp), dimension(3), intent(in) :: mf, dalpha_rho_ds
168-
real(wp), dimension(3), intent(in) :: dvel_ds
169-
real(wp), dimension(3), intent(in) :: dadv_ds
170-
real(wp), dimension(10), intent(in) :: dYs_ds
167+
real(wp), dimension(3), intent(in) :: mf, dalpha_rho_ds
168+
real(wp), dimension(3), intent(in) :: dvel_ds
169+
real(wp), dimension(3), intent(in) :: dadv_ds
170+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$), intent(in) :: dYs_ds
171171
#:else
172172
real(wp), dimension(num_fluids), intent(in) :: mf, dalpha_rho_ds
173173
real(wp), dimension(num_dims), intent(in) :: dvel_ds
@@ -199,7 +199,7 @@ contains
199199

200200
real(wp), dimension(3), intent(in) :: lambda
201201
#:if USING_AMD
202-
real(wp), dimension(20), intent(inout) :: L
202+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$), intent(inout) :: L
203203
#:else
204204
real(wp), dimension(sys_size), intent(inout) :: L
205205
#:endif
@@ -223,15 +223,15 @@ contains
223223

224224
real(wp), dimension(3), intent(in) :: lambda
225225
#:if USING_AMD
226-
real(wp), dimension(20), intent(inout) :: L
226+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$), intent(inout) :: L
227227
#:else
228228
real(wp), dimension(sys_size), intent(inout) :: L
229229
#:endif
230230
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
231-
real(wp), dimension(3), intent(in) :: mf, dalpha_rho_ds
232-
real(wp), dimension(3), intent(in) :: dvel_ds
233-
real(wp), dimension(3), intent(in) :: dadv_ds
234-
real(wp), dimension(10), intent(in) :: dYs_ds
231+
real(wp), dimension(3), intent(in) :: mf, dalpha_rho_ds
232+
real(wp), dimension(3), intent(in) :: dvel_ds
233+
real(wp), dimension(3), intent(in) :: dadv_ds
234+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$), intent(in) :: dYs_ds
235235
#:else
236236
real(wp), dimension(num_fluids), intent(in) :: mf, dalpha_rho_ds
237237
real(wp), dimension(num_dims), intent(in) :: dvel_ds
@@ -257,7 +257,7 @@ contains
257257

258258
real(wp), dimension(3), intent(in) :: lambda
259259
#:if USING_AMD
260-
real(wp), dimension(20), intent(inout) :: L
260+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$), intent(inout) :: L
261261
#:else
262262
real(wp), dimension(sys_size), intent(inout) :: L
263263
#:endif
@@ -288,7 +288,7 @@ contains
288288

289289
real(wp), dimension(3), intent(in) :: lambda
290290
#:if USING_AMD
291-
real(wp), dimension(20), intent(inout) :: L
291+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$), intent(inout) :: L
292292
#:else
293293
real(wp), dimension(sys_size), intent(inout) :: L
294294
#:endif
@@ -317,7 +317,7 @@ contains
317317

318318
$:GPU_ROUTINE(function_name='s_compute_supersonic_inflow_L', parallelism='[seq]', cray_inline=True)
319319
#:if USING_AMD
320-
real(wp), dimension(20), intent(inout) :: L
320+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$), intent(inout) :: L
321321
#:else
322322
real(wp), dimension(sys_size), intent(inout) :: L
323323
#:endif
@@ -333,15 +333,15 @@ contains
333333

334334
real(wp), dimension(3), intent(in) :: lambda
335335
#:if USING_AMD
336-
real(wp), dimension(20), intent(inout) :: L
336+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$), intent(inout) :: L
337337
#:else
338338
real(wp), dimension(sys_size), intent(inout) :: L
339339
#:endif
340340
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
341-
real(wp), dimension(3), intent(in) :: mf, dalpha_rho_ds
342-
real(wp), dimension(3), intent(in) :: dvel_ds
343-
real(wp), dimension(3), intent(in) :: dadv_ds
344-
real(wp), dimension(10), intent(in) :: dYs_ds
341+
real(wp), dimension(3), intent(in) :: mf, dalpha_rho_ds
342+
real(wp), dimension(3), intent(in) :: dvel_ds
343+
real(wp), dimension(3), intent(in) :: dadv_ds
344+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$), intent(in) :: dYs_ds
345345
#:else
346346
real(wp), dimension(num_fluids), intent(in) :: mf, dalpha_rho_ds
347347
real(wp), dimension(num_dims), intent(in) :: dvel_ds

‎src/simulation/m_ibm.fpp‎

Lines changed: 6 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -157,12 +157,12 @@ contains
157157
real(wp) :: c_IP
158158

159159
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
160-
real(wp), dimension(3) :: Gs
161-
real(wp), dimension(3) :: alpha_rho_IP, alpha_IP
162-
real(wp), dimension(3) :: r_IP, v_IP, pb_IP, mv_IP
163-
real(wp), dimension(18) :: nmom_IP
164-
real(wp), dimension(12) :: presb_IP, massv_IP
165-
real(wp), dimension(10) :: Ys_IP
160+
real(wp), dimension(3) :: Gs
161+
real(wp), dimension(3) :: alpha_rho_IP, alpha_IP
162+
real(wp), dimension(3) :: r_IP, v_IP, pb_IP, mv_IP
163+
real(wp), dimension(18) :: nmom_IP
164+
real(wp), dimension(12) :: presb_IP, massv_IP
165+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: Ys_IP
166166
#:else
167167
real(wp), dimension(num_fluids) :: Gs
168168
real(wp), dimension(num_fluids) :: alpha_rho_IP, alpha_IP

‎src/simulation/m_riemann_solver_hll.fpp‎

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -41,11 +41,11 @@ contains
4141
type(int_bounds_info), intent(in) :: ix, iy, iz
4242

4343
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
44-
real(wp), dimension(3) :: alpha_rho_L, alpha_rho_R
45-
real(wp), dimension(3) :: vel_L, vel_R
46-
real(wp), dimension(3) :: alpha_L, alpha_R
47-
real(wp), dimension(10) :: Ys_L, Ys_R, R_species, h_iL, h_iR
48-
real(wp), dimension(10) :: Cp_iL, Cp_iR, Xs_L, Xs_R, Gamma_iL, Gamma_iR
44+
real(wp), dimension(3) :: alpha_rho_L, alpha_rho_R
45+
real(wp), dimension(3) :: vel_L, vel_R
46+
real(wp), dimension(3) :: alpha_L, alpha_R
47+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: Ys_L, Ys_R, R_species, h_iL, h_iR
48+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: Cp_iL, Cp_iR, Xs_L, Xs_R, Gamma_iL, Gamma_iR
4949
#:else
5050
real(wp), dimension(num_fluids) :: alpha_rho_L, alpha_rho_R
5151
real(wp), dimension(num_vels) :: vel_L, vel_R

‎src/simulation/m_riemann_solver_hllc.fpp‎

Lines changed: 5 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -60,7 +60,8 @@ contains
6060
real(wp) :: E_L, E_R
6161
real(wp) :: H_L, H_R
6262
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
63-
real(wp), dimension(10) :: Ys_L, Ys_R, Xs_L, Xs_R, Gamma_iL, Gamma_iR, Cp_iL, Cp_iR, R_species, h_iL, h_iR
63+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: Ys_L, Ys_R, Xs_L, Xs_R, Gamma_iL, Gamma_iR, Cp_iL, Cp_iR, R_species, &
64+
& h_iL, h_iR
6465
#:else
6566
real(wp), dimension(num_species) :: Ys_L, Ys_R, Xs_L, Xs_R, Gamma_iL, Gamma_iR, Cp_iL, Cp_iR, R_species, h_iL, h_iR
6667
#:endif
@@ -118,8 +119,8 @@ contains
118119

119120
! HLLC star-state helpers
120121
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
121-
real(wp), dimension(20) :: U_L, U_R
122-
real(wp), dimension(20) :: F_L, F_R, F_star_L, F_star_R, F_HLLC
122+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$) :: U_L, U_R
123+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$) :: F_L, F_R, F_star_L, F_star_R, F_HLLC
123124
#:else
124125
real(wp), dimension(sys_size) :: U_L, U_R
125126
real(wp), dimension(sys_size) :: F_L, F_R, F_star_L, F_star_R, F_HLLC
@@ -140,7 +141,7 @@ contains
140141

141142
! ADC (HLL -> HLLC)
142143
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
143-
real(wp), dimension(20) :: F_HLL
144+
real(wp), dimension(${AMD_SYS_SIZE_MAX}$) :: F_HLL
144145
#:else
145146
real(wp), dimension(sys_size) :: F_HLL
146147
#:endif

‎src/simulation/m_riemann_solver_lf.fpp‎

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -37,11 +37,11 @@ contains
3737
type(int_bounds_info), intent(in) :: ix, iy, iz
3838

3939
#:if not MFC_CASE_OPTIMIZATION and USING_AMD
40-
real(wp), dimension(3) :: alpha_rho_L, alpha_rho_R
41-
real(wp), dimension(3) :: vel_L, vel_R
42-
real(wp), dimension(3) :: alpha_L, alpha_R
43-
real(wp), dimension(10) :: Ys_L, Ys_R
44-
real(wp), dimension(10) :: Cp_iL, Cp_iR, Xs_L, Xs_R, Gamma_iL, Gamma_iR
40+
real(wp), dimension(3) :: alpha_rho_L, alpha_rho_R
41+
real(wp), dimension(3) :: vel_L, vel_R
42+
real(wp), dimension(3) :: alpha_L, alpha_R
43+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: Ys_L, Ys_R
44+
real(wp), dimension(${AMD_NUM_SPECIES_MAX}$) :: Cp_iL, Cp_iR, Xs_L, Xs_R, Gamma_iL, Gamma_iR
4545
real(wp), dimension(3, 3) :: vel_grad_L, vel_grad_R !< Averaged velocity gradient tensor `d(vel_i)/d(coord_j)`.
4646
#:else
4747
real(wp), dimension(num_fluids) :: alpha_rho_L, alpha_rho_R

0 commit comments

Comments
 (0)