From 9c30051663f40d7b4313c76b44f9a583d90b5612 Mon Sep 17 00:00:00 2001 From: JGonzalez Date: Thu, 26 Feb 2026 15:23:10 +0100 Subject: [PATCH] Forgot to link edge to boundary procedures. 1D case working right now --- src/modules/mesh/1DCart/moduleMesh1DCart.f90 | 7 +++++++ src/modules/mesh/1DRad/moduleMesh1DRad.f90 | 7 +++++++ src/modules/mesh/2DCart/moduleMesh2DCart.f90 | 7 +++++++ src/modules/mesh/2DCyl/moduleMesh2DCyl.f90 | 7 +++++++ src/modules/mesh/3DCart/moduleMesh3DCart.f90 | 7 +++++++ 5 files changed, 35 insertions(+) diff --git a/src/modules/mesh/1DCart/moduleMesh1DCart.f90 b/src/modules/mesh/1DCart/moduleMesh1DCart.f90 index 41519e3..1240c03 100644 --- a/src/modules/mesh/1DCart/moduleMesh1DCart.f90 +++ b/src/modules/mesh/1DCart/moduleMesh1DCart.f90 @@ -130,6 +130,13 @@ MODULE moduleMesh1DCart call meshEdgePointer_add(physicalSurfaces(ps_index)%edges, self%n) call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n1%n) + ! Link edge to particle boundaries + allocate(self%boundariesParticle(1:nSpecies)) + do s = 1, nSpecies + self%boundariesParticle(s)%obj => physicalSurfaces(ps_index)%particles(s)%obj + + end do + end subroutine initEdge1DCart !Get nodes from edge diff --git a/src/modules/mesh/1DRad/moduleMesh1DRad.f90 b/src/modules/mesh/1DRad/moduleMesh1DRad.f90 index d32874e..e250bb2 100644 --- a/src/modules/mesh/1DRad/moduleMesh1DRad.f90 +++ b/src/modules/mesh/1DRad/moduleMesh1DRad.f90 @@ -130,6 +130,13 @@ MODULE moduleMesh1DRad call meshEdgePointer_add(physicalSurfaces(ps_index)%edges, self%n) call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n1%n) + ! Link edge to particle boundaries + allocate(self%boundariesParticle(1:nSpecies)) + do s = 1, nSpecies + self%boundariesParticle(s)%obj => physicalSurfaces(ps_index)%particles(s)%obj + + end do + end subroutine initEdge1DRad !Get nodes from edge diff --git a/src/modules/mesh/2DCart/moduleMesh2DCart.f90 b/src/modules/mesh/2DCart/moduleMesh2DCart.f90 index 1a998bc..7a84d15 100644 --- a/src/modules/mesh/2DCart/moduleMesh2DCart.f90 +++ b/src/modules/mesh/2DCart/moduleMesh2DCart.f90 @@ -170,6 +170,13 @@ MODULE moduleMesh2DCart call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n1%n) call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n2%n) + ! Link edge to particle boundaries + allocate(self%boundariesParticle(1:nSpecies)) + do s = 1, nSpecies + self%boundariesParticle(s)%obj => physicalSurfaces(ps_index)%particles(s)%obj + + end do + END SUBROUTINE initEdge2DCart !Get nodes from edge diff --git a/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90 b/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90 index a85bb4f..f14b2c8 100644 --- a/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90 +++ b/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90 @@ -179,6 +179,13 @@ MODULE moduleMesh2DCyl call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n1%n) call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n2%n) + ! Link edge to particle boundaries + allocate(self%boundariesParticle(1:nSpecies)) + do s = 1, nSpecies + self%boundariesParticle(s)%obj => physicalSurfaces(ps_index)%particles(s)%obj + + end do + END SUBROUTINE initEdge2DCyl !Get nodes from edge diff --git a/src/modules/mesh/3DCart/moduleMesh3DCart.f90 b/src/modules/mesh/3DCart/moduleMesh3DCart.f90 index a9310e7..a89f51c 100644 --- a/src/modules/mesh/3DCart/moduleMesh3DCart.f90 +++ b/src/modules/mesh/3DCart/moduleMesh3DCart.f90 @@ -153,6 +153,13 @@ MODULE moduleMesh3DCart call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n2%n) call meshNodePointer_add(physicalSurfaces(ps_index)%nodes, self%n3%n) + ! Link edge to particle boundaries + allocate(self%boundariesParticle(1:nSpecies)) + do s = 1, nSpecies + self%boundariesParticle(s)%obj => physicalSurfaces(ps_index)%particles(s)%obj + + end do + END SUBROUTINE initEdge3DCartTria !Get nodes from surface