223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

Initial nodal and ip coordinates added to the HDF5 output file

This commit is contained in:
Abisheik Panneerselvam 2021-05-31 15:01:42 +02:00
parent 718880f177
commit 650fbd7509
1 changed files with 2 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
src/mesh/discretization_mesh.f90
View File

@ -144,10 +144,10 @@ subroutine discretization_mesh_init(restart)
CHKERRQ(ierr) CHKERRQ(ierr)
! Get initial nodal coordinates ! Get initial nodal coordinates
call DMGetCoordinates(geomMesh,coords,ierr) call DMGetCoordinates(geomMesh,coords_node0,ierr)
CHKERRQ(ierr) CHKERRQ(ierr)
allocate(mesh_node0_temp(dimPlex*mesh_Nnodes)) allocate(mesh_node0_temp(dimPlex*mesh_Nnodes))
call VecGetArrayF90(coords, mesh_node0_temp,ierr) call VecGetArrayF90(coords_node0, mesh_node0_temp,ierr)
CHKERRQ(ierr) CHKERRQ(ierr)
mesh_maxNips = FEM_nQuadrature(dimPlex,integrationOrder) 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) allocate(mesh_node0(3,mesh_Nnodes),source=0.0_pReal)
mesh_node0 = reshape(mesh_node0_temp,[dimPlex,mesh_Nnodes]) mesh_node0 = reshape(mesh_node0_temp,[dimPlex,mesh_Nnodes])
deallocate(mesh_node0_temp)
call discretization_init(materialAt,& call discretization_init(materialAt,&
reshape(mesh_ipCoordinates,[3,mesh_maxNips*mesh_NcpElems]), & reshape(mesh_ipCoordinates,[3,mesh_maxNips*mesh_NcpElems]), &