Skip to content

Commit

Permalink
Update m_ibm.fpp
Browse files Browse the repository at this point in the history
  • Loading branch information
@anshgupta1234
anshgupta1234 committed Nov 5, 2024
1 parent 9430741 commit 10628d3
Showing 1 changed file with 0 additions and 38 deletions.
38 changes: 0 additions & 38 deletions src/simulation/m_ibm.fpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -578,20 +578,6 @@ contains
inner_points(count_i)%ib_patch_id = &
patch_id
inner_points(count_i)%slip = patch_ib(patch_id)%slip
if ((x_cc(i) - dx(i)) < x_domain%beg .or. &
(x_cc(i) + dx(i)) > x_domain%end) then
ghost_points(count)%DB(1) = 1
else
ghost_points(count)%DB(1) = 0
end if

if ((y_cc(j) - dy(j)) < y_domain%beg .or. &
(y_cc(j) + dy(j)) > y_domain%end) then
ghost_points(count)%DB(2) = 1
else
ghost_points(count)%DB(2) = 0
end if

count_i = count_i + 1

end if
Expand Down Expand Up @@ -642,30 +628,6 @@ contains
ib_markers%sf(i, j, k)
inner_points(count_i)%slip = patch_ib(patch_id)%slip

if ((x_cc(i) - dx(i)) < x_domain%beg) then
ghost_points(count)%DB(1) = -1
else if ((x_cc(i) + dx(i)) > x_domain%end) then
ghost_points(count)%DB(1) = 1
else
ghost_points(count)%DB(1) = 0
end if

if ((y_cc(j) - dy(j)) < y_domain%beg) then
ghost_points(count)%DB(2) = -1
else if ((y_cc(j) + dy(j)) > y_domain%end) then
ghost_points(count)%DB(2) = 1
else
ghost_points(count)%DB(2) = 0
end if

if ((z_cc(k) - dz(k)) < z_domain%beg) then
ghost_points(count)%DB(3) = -1
else if ((z_cc(k) + dz(k)) > z_domain%end) then
ghost_points(count)%DB(3) = 1
else
ghost_points(count)%DB(3) = 0
end if

count_i = count_i + 1
end if
end if
Expand Down

0 comments on commit 10628d3

Please sign in to comment.