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

Still working on it

Browse source 
No uniform density yet...
This commit is contained in:
Jorge Gonzalez 2024-07-06 10:12:03 +02:00
commit 6b5ac16e4b
2 changed files with 1 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

1
src/modules/mesh/2DCyl/moduleMesh2DCyl.f90
View file

@ -235,7 +235,6 @@ MODULE moduleMesh2DCyl
p1 = (/self%z(1), self%r(1) /) p1 = (/self%z(1), self%r(1) /)
p2 = (/self%z(2), self%r(2) /) p2 = (/self%z(2), self%r(2) /)
r(1:2) = (1.D0 - rnd)*p1 + rnd*p2 r(1:2) = (1.D0 - rnd)*p1 + rnd*p2
r(2) = (self%r(2) + self%r(1)) * 0.5D0
r(3) = 0.D0 r(3) = 0.D0
END FUNCTION randPosEdge END FUNCTION randPosEdge

2
src/modules/moduleInject.f90
View file

@ -318,7 +318,7 @@ MODULE moduleInject
!Random position in edge !Random position in edge
partInj(n)%r = randomEdge%randPos() partInj(n)%r = randomEdge%randPos()
!Assign weight to particle. !Assign weight to particle.
partInj(n)%weight = self%species%weight * sqrt(partInj(n)%r(2) / (0.1D0*L_ref)) partInj(n)%weight = self%species%weight * sqrt(partInj(n)%r(2) / (0.1D0/L_ref))
!Volume associated to the edge: !Volume associated to the edge:
IF (ASSOCIATED(randomEdge%e1)) THEN IF (ASSOCIATED(randomEdge%e1)) THEN
partInj(n)%cell = randomEdge%e1%n partInj(n)%cell = randomEdge%e1%n