From ba2b5bc6f3f9c7ee4b1a3f9f627c34b42eee6fef Mon Sep 17 00:00:00 2001 From: Jorge Gonzalez Date: Tue, 7 Feb 2023 09:43:28 +0100 Subject: [PATCH] Fixed issue with initial velocity The problem was that the number of values for velocity was numNodes and not 3*numNodes. --- src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 b/src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 index c62dec1..f58564e 100644 --- a/src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 +++ b/src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 @@ -530,7 +530,7 @@ MODULE moduleMeshInputVTU !Read the species velocity line = findLine(fileID, 'Name="Velocity') ALLOCATE(velocityBlock(1:3*numNodes)) - CALL readDataBlock(fileID, numNodes, velocityBlock) + CALL readDataBlock(fileID, 3*numNodes, velocityBlock) ALLOCATE(velocity(1:numNodes, 1:3)) DO n = 1, numNodes velocity(n, 1) = velocityBlock(3*(n-1)+1)