Structure for Neumann BC

src/modules/mesh/moduleMesh.f90

@@ -1056,6 +1056,33 @@ MODULE moduleMesh
 
   end interface
 
+  ! Neumann boundary condition (constant electric field normal to edge)
+  type, extends(boundaryEMGeneric):: boundaryEMNeumann
+    real(8):: electricField ! Electric field normal to the edge
+    contains
+      ! boundaryEMGeneric deferred procedures
+      procedure, pass:: apply => applyNeumann
+
+  end type boundaryEMNeumann
+
+  interface
+    module subroutine initNeumann(self, config, object)
+      use json_module
+
+      class(boundaryEMGeneric), allocatable, intent(inout):: self
+      type(json_file), intent(inout):: config
+      character(:), allocatable, intent(in):: object
+
+    end subroutine initNeumann
+
+    module subroutine applyNeumann(self, vectorF)
+      class(boundaryEMNeumann), intent(in):: self
+      real(8), intent(inout):: vectorF(:)
+
+    end subroutine applyNeumann
+
+  end interface
+
   ! Floating: Floating surface, ie, zero net current
   type, extends(boundaryEMGeneric):: boundaryEMFloating
     real(8):: potential

src/modules/mesh/moduleMesh@boundaryEM.f90

@@ -123,6 +123,27 @@ submodule(moduleMesh) boundaryEM
 
     end subroutine updateDirichletTime
 
+    ! Neumann constant value
+    ! Init
+    module subroutine initNeumann(self, config, object)
+      use json_module
+      implicit none
+
+      class(boundaryEMGeneric), allocatable, intent(inout):: self
+      type(json_file), intent(inout):: config
+      character(:), allocatable, intent(in):: object
+
+    end subroutine initNeumann
+
+    ! Apply
+    module subroutine applyNeumann(self, vectorF)
+      implicit none
+
+      class(boundaryEMNeumann), intent(in):: self
+      real(8), intent(inout):: vectorF(:)
+
+    end subroutine applyNeumann
+
     ! Floating surface
     ! Init
     module subroutine initFloating(self, config, object)