Still not working

This commit is contained in:
Jorge Gonzalez 2026-02-04 20:15:05 +01:00
commit a1d5e20a09

View file

@ -232,13 +232,8 @@ MODULE moduleMeshBoundary
select type(bound => edge%boundary%bTypes(part%species%n)%obj)
type is(boundaryOutflowAdaptive)
if (dot_product(part%v,-edge%normal) > v_cut) then
! print *,part%v(1), v_cut
! part%v = part%v + v_cut*edge%normal
part%v(1) = part%v(1) - v_cut
! print *,part%v(1)
part%v = part%v + v_cut*edge%normal
call reflection(edge, part)
! print *,part%v(1)
! print *
else
call transparent(edge, part)