Fixed problem with edges with more than 1 node

This commit is contained in:
Jorge Gonzalez 2026-05-08 20:53:30 +02:00
commit 1d099993c9

View file

@ -262,7 +262,7 @@ submodule(moduleMesh) boundaryEM
allocate(mom_nodes(1:edge%nNodes))
do s = 1, nSpecies
mom_center = 0.0d0
do n = 1, self%nNodes
do n = 1, edge%nNodes
node => mesh%nodes(nodes(n))%obj
! Minus sign to get the flux exiting the surface (opposite to the normal)
@ -405,7 +405,7 @@ submodule(moduleMesh) boundaryEM
allocate(mom_nodes(1:edge%nNodes))
do s = 1, nSpecies
do n = 1, self%nNodes
do n = 1, edge%nNodes
node => mesh%nodes(nodes(n))%obj
! Minus sign to get the flux exiting the surface (opposite to the normal)