From 152ae6b097422d5b5e924835627b5d52473c5c5b Mon Sep 17 00:00:00 2001 From: JGonzalez Date: Thu, 11 Jul 2024 18:19:26 +0200 Subject: [PATCH] No longer needed The correction in the node volume is no longer needed as now things are being calculated right with the last change. Still, at some point I should review the calculation of the node volume in 2DCyl. --- .../mesh/inout/gmsh2/moduleMeshInputGmsh2.f90 | 14 -------------- src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 | 14 -------------- 2 files changed, 28 deletions(-) diff --git a/src/modules/mesh/inout/gmsh2/moduleMeshInputGmsh2.f90 b/src/modules/mesh/inout/gmsh2/moduleMeshInputGmsh2.f90 index 1c8de50..1ad6a36 100644 --- a/src/modules/mesh/inout/gmsh2/moduleMeshInputGmsh2.f90 +++ b/src/modules/mesh/inout/gmsh2/moduleMeshInputGmsh2.f90 @@ -297,20 +297,6 @@ MODULE moduleMeshInputGmsh2 CLOSE(10) - ! Adjust node volume at axis - SELECT CASE(self%geometry) - CASE("Cyl") - DO n = 1, self%numNodes - r = self%nodes(n)%obj%getCoordinates() - IF (r(2) == 0.D0) THEN - self%nodes(n)%obj%v = self%nodes(n)%obj%v * 3.0D0/2.0D0 - - END IF - - END DO - - END SELECT - !Call mesh connectivity CALL self%connectMesh diff --git a/src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 b/src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 index d408387..184d645 100644 --- a/src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 +++ b/src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 @@ -497,20 +497,6 @@ MODULE moduleMeshInputVTU END DO - ! Adjust node volume at axis - SELECT CASE(self%geometry) - CASE("Cyl") - DO n = 1, self%numNodes - r = self%nodes(n)%obj%getCoordinates() - IF (r(2) == 0.D0) THEN - self%nodes(n)%obj%v = self%nodes(n)%obj%v * 3.0D0/2.0D0 - - END IF - - END DO - - END SELECT - !Call mesh connectivity CALL self%connectMesh