JGonzalez/fpakc
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

No longer needed

Browse source 
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.
This commit is contained in:
Jorge Gonzalez 2024-07-11 18:19:26 +02:00
commit 152ae6b097
2 changed files with 0 additions and 28 deletions
Download patch file Download diff file Expand all files Collapse all files

14
src/modules/mesh/inout/gmsh2/moduleMeshInputGmsh2.f90
View file

@ -297,20 +297,6 @@ MODULE moduleMeshInputGmsh2
CLOSE(10) 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 mesh connectivity
CALL self%connectMesh CALL self%connectMesh

14
src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90
View file

@ -497,20 +497,6 @@ MODULE moduleMeshInputVTU
END DO 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 mesh connectivity
CALL self%connectMesh CALL self%connectMesh