From e9ac9bdd80668e75214e43d4ac6b21884573e086 Mon Sep 17 00:00:00 2001 From: Jorge Gonzalez Date: Tue, 7 Feb 2023 15:19:13 +0100 Subject: [PATCH] Finishing implementation of vtu mesh format I have tested all geometries and cases and all seem to work perfectly with .vtu meshes. Input and output works great, starting from a previous case also works. Everything I was able to test is okay. Still, what I want to do is to change a few things in the output (e.g., change OUTPUT prefix to Step) and try to improve reading gmsh meshes to make them more compact as vtu is now. --- src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 b/src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 index f58564e..763517f 100644 --- a/src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 +++ b/src/modules/mesh/inout/vtu/moduleMeshInputVTU.f90 @@ -334,7 +334,7 @@ MODULE moduleMeshInputVTU CASE(2) IF (types(n) == 3) THEN - e = e+1 + e = e + 1 ALLOCATE(p(1:2)) p(1) = connectivity(offsets(n) - 1) p(2) = connectivity(offsets(n)) @@ -359,7 +359,7 @@ MODULE moduleMeshInputVTU END IF CASE(1) - IF (types(n) == 3) THEN + IF (types(n) == 1) THEN e = e + 1 ALLOCATE(p(1:1)) p(1) = connectivity(offsets(n))