JGonzalez/fpakc
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

All EM boundaries should contain the nodes they apply to, regardless of the type

Browse source
This commit is contained in:
Jorge Gonzalez 2026-04-14 10:50:40 +02:00
commit fa5f037d41
1 changed files with 7 additions and 33 deletions
Download patch file Download diff file Expand all files Collapse all files

40
src/modules/init/moduleInput.f90
View file

@ -1346,40 +1346,14 @@ MODULE moduleInput
! EM Boundaries
do b = 1, nBoundariesEM
associate(bound => boundariesEM(b)%obj)
select type(bound)
type is(boundaryEMDirichlet)
! Loop over all physical surfaces
do ps = 1, nPhysicalSurfaces
! If the boundary for the species is linked to the one analysing, add the edges
if (associated(physicalSurfaces(ps)%EM, bound)) then
bound%nodes = [bound%nodes, physicalSurfaces(ps)%nodes]
end if
! Loop over all physical surfaces
do ps = 1, nPhysicalSurfaces
! If the boundary for the species is linked to the one analysing, add the edges
if (associated(physicalSurfaces(ps)%EM, bound)) then
bound%nodes = [bound%nodes, physicalSurfaces(ps)%nodes]
end if
end do
type is(boundaryEMDirichletTime)
! Loop over all physical surfaces
do ps = 1, nPhysicalSurfaces
! If the boundary for the species is linked to the one analysing, add the edges
if (associated(physicalSurfaces(ps)%EM, bound)) then
bound%nodes = [bound%nodes, physicalSurfaces(ps)%nodes]
end if
end do
type is(boundaryEMFloating)
! Loop over all physical surfaces
do ps = 1, nPhysicalSurfaces
! If the boundary for the species is linked to the one analysing, add the edges
if (associated(physicalSurfaces(ps)%EM, bound)) then
bound%nodes = [bound%nodes, physicalSurfaces(ps)%nodes]
bound%edges = [bound%edges, physicalSurfaces(ps)%edges]
end if
end do
end select
end do
bound%nNodes = size(bound%nodes)