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

Fixed issue with initial velocity

Browse source 
The problem was that the number of values for velocity was numNodes and
not 3*numNodes.
This commit is contained in:
Jorge Gonzalez 2023-02-07 09:43:28 +01:00
commit ba2b5bc6f3
1 changed files with 1 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

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

@ -530,7 +530,7 @@ MODULE moduleMeshInputVTU
!Read the species velocity !Read the species velocity
line = findLine(fileID, 'Name="Velocity') line = findLine(fileID, 'Name="Velocity')
ALLOCATE(velocityBlock(1:3*numNodes)) ALLOCATE(velocityBlock(1:3*numNodes))
CALL readDataBlock(fileID, numNodes, velocityBlock) CALL readDataBlock(fileID, 3*numNodes, velocityBlock)
ALLOCATE(velocity(1:numNodes, 1:3)) ALLOCATE(velocity(1:numNodes, 1:3))
DO n = 1, numNodes DO n = 1, numNodes
velocity(n, 1) = velocityBlock(3*(n-1)+1) velocity(n, 1) = velocityBlock(3*(n-1)+1)