Deleted the constraing of flux > 0

This commit is contained in:
Jorge Gonzalez 2026-04-13 14:54:40 +02:00
commit e5519cee2d

View file

@ -224,9 +224,6 @@ submodule(moduleMesh) boundaryEM
mom_center = edge%gatherF(edge%centerXi(), edge%nNodes, mom_nodes)
! Only account for charge exiting the surface
mom_center = max(mom_center, 0.0d0)
edgeDensityCurrent = edgeDensityCurrent + qSpecies(s) * mom_center
end do