From 650fbd75099abcdd63916bbbad0f549947292f7d Mon Sep 17 00:00:00 2001 From: Abisheik Panneerselvam Date: Mon, 31 May 2021 15:01:42 +0200 Subject: [PATCH] Initial nodal and ip coordinates added to the HDF5 output file --- src/mesh/discretization_mesh.f90 | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/src/mesh/discretization_mesh.f90 b/src/mesh/discretization_mesh.f90 index 349e66ff2..af3021be0 100644 --- a/src/mesh/discretization_mesh.f90 +++ b/src/mesh/discretization_mesh.f90 @@ -144,10 +144,10 @@ subroutine discretization_mesh_init(restart) CHKERRQ(ierr) ! Get initial nodal coordinates - call DMGetCoordinates(geomMesh,coords,ierr) + call DMGetCoordinates(geomMesh,coords_node0,ierr) CHKERRQ(ierr) allocate(mesh_node0_temp(dimPlex*mesh_Nnodes)) - call VecGetArrayF90(coords, mesh_node0_temp,ierr) + call VecGetArrayF90(coords_node0, mesh_node0_temp,ierr) CHKERRQ(ierr) mesh_maxNips = FEM_nQuadrature(dimPlex,integrationOrder) @@ -167,7 +167,6 @@ subroutine discretization_mesh_init(restart) allocate(mesh_node0(3,mesh_Nnodes),source=0.0_pReal) mesh_node0 = reshape(mesh_node0_temp,[dimPlex,mesh_Nnodes]) - deallocate(mesh_node0_temp) call discretization_init(materialAt,& reshape(mesh_ipCoordinates,[3,mesh_maxNips*mesh_NcpElems]), &