Finally! Dividing by the volume of the node gave the right value

This commit is contained in:
Jorge Gonzalez 2026-04-15 16:56:36 +02:00
commit 06a797de49

View file

@ -163,7 +163,7 @@ submodule(moduleMesh) boundaryEM
do n = 1, self%nNodes
associate(node => self%nodes(n)%obj)
vectorF(node%n) = vectorF(node%n) - self%electricField
vectorF(node%n) = vectorF(node%n) - self%electricField/node%v
end associate