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

Some parallel issues

Browse source 
Forgot to check Gmsh2 format with new changes.

Cartesian coordinates were not calculating things properly.
This commit is contained in:
Jorge Gonzalez 2024-07-07 11:53:56 +02:00
commit 5bc064d018
2 changed files with 3 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

3
src/modules/mesh/2DCart/moduleMesh2DCart.f90
View file

@ -556,6 +556,7 @@ MODULE moduleMesh2DCart
!Compute element volume !Compute element volume
PURE SUBROUTINE volumeQuad(self) PURE SUBROUTINE volumeQuad(self)
USE moduleRefParam, ONLY: L_ref
IMPLICIT NONE IMPLICIT NONE
CLASS(meshCell2DCartQuad), INTENT(inout):: self CLASS(meshCell2DCartQuad), INTENT(inout):: self
@ -573,7 +574,7 @@ MODULE moduleMesh2DCart
fPsi = self%fPsi(Xi, 4) fPsi = self%fPsi(Xi, 4)
!Compute total volume of the cell !Compute total volume of the cell
self%volume = detJ*4.D0 self%volume = detJ*4.D0/L_ref
!Compute volume per node !Compute volume per node
self%n1%v = self%n1%v + fPsi(1)*self%volume self%n1%v = self%n1%v + fPsi(1)*self%volume
self%n2%v = self%n2%v + fPsi(2)*self%volume self%n2%v = self%n2%v + fPsi(2)*self%volume

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

@ -299,7 +299,7 @@ MODULE moduleMeshInputGmsh2
! Adjust node volume at axis ! Adjust node volume at axis
SELECT CASE(self%geometry) SELECT CASE(self%geometry)
CASE("Cyl") CASE("Cyl")
DO n = 1, numNodes DO n = 1, self%numNodes
r = self%nodes(n)%obj%getCoordinates() r = self%nodes(n)%obj%getCoordinates()
IF (r(2) == 0.D0) THEN IF (r(2) == 0.D0) THEN
self%nodes(n)%obj%v = self%nodes(n)%obj%v * 2.0D0!2.0D0/3.0D0 self%nodes(n)%obj%v = self%nodes(n)%obj%v * 2.0D0!2.0D0/3.0D0