Skip to content

Commit af68ebb

Browse files
anandrdbzAnand
andauthored
Fixing an IBM bug on Frontier (#892)
Co-authored-by: Anand <anand@ipsec-10-2-195-71.vpn.gatech.edu>
1 parent e89abc6 commit af68ebb

File tree

1 file changed

+1
-27
lines changed

1 file changed

+1
-27
lines changed

src/simulation/m_ibm.fpp

Lines changed: 1 addition & 27 deletions
Original file line numberDiff line numberDiff line change
@@ -313,40 +313,14 @@ contains
313313
!$acc parallel loop gang vector private(physical_loc, dyn_pres, alpha_rho_IP, alpha_IP, vel_g, rho, gamma, pi_inf, Re_K, innerp, j, k, l, q)
314314
do i = 1, num_inner_gps
315315

316-
vel_g = 0._wp
317316
innerp = inner_points(i)
318317
j = innerp%loc(1)
319318
k = innerp%loc(2)
320319
l = innerp%loc(3)
321-
patch_id = inner_points(i)%ib_patch_id
322-
323-
! Calculate physical location of GP
324-
if (p > 0) then
325-
physical_loc = [x_cc(j), y_cc(k), z_cc(l)]
326-
else
327-
physical_loc = [x_cc(j), y_cc(k), 0._wp]
328-
end if
329-
330-
!$acc loop seq
331-
do q = 1, num_fluids
332-
q_prim_vf(q)%sf(j, k, l) = alpha_rho_IP(q)
333-
q_prim_vf(advxb + q - 1)%sf(j, k, l) = alpha_IP(q)
334-
end do
335-
336-
if (surface_tension) then
337-
q_prim_vf(c_idx)%sf(j, k, l) = c_IP
338-
end if
339-
340-
call s_convert_species_to_mixture_variables_acc(rho, gamma, pi_inf, qv_K, alpha_IP, &
341-
alpha_rho_IP, Re_K)
342-
343-
dyn_pres = 0._wp
344320

345321
!$acc loop seq
346322
do q = momxb, momxe
347-
q_cons_vf(q)%sf(j, k, l) = rho*vel_g(q - momxb + 1)
348-
dyn_pres = dyn_pres + q_cons_vf(q)%sf(j, k, l)* &
349-
vel_g(q - momxb + 1)/2._wp
323+
q_cons_vf(q)%sf(j, k, l) = 0._wp
350324
end do
351325
end do
352326

0 commit comments

Comments
 (0)