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

writing initial coordinates to DADF5

This commit is contained in:
Martin Diehl 2019-10-12 15:03:26 +02:00
parent 513f1c6726
commit 734e6ef15f
1 changed files with 12 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/mesh_marc.f90
View File

@ -124,6 +124,9 @@ subroutine mesh_init(ip,el)
real(pReal),dimension(:,:,:,:), allocatable :: & real(pReal),dimension(:,:,:,:), allocatable :: &
mesh_ipAreaNormal !< area normal of interface to neighboring IP (initially!) mesh_ipAreaNormal !< area normal of interface to neighboring IP (initially!)
real(pReal), dimension(:,:), allocatable :: &
ip_reshaped
write(6,'(/,a)') ' <<<+- mesh init -+>>>' write(6,'(/,a)') ' <<<+- mesh init -+>>>'
mesh_unitlength = numerics_unitlength ! set physical extent of a length unit in mesh mesh_unitlength = numerics_unitlength ! set physical extent of a length unit in mesh
@ -197,10 +200,18 @@ subroutine mesh_init(ip,el)
calcMode = .false. ! pretend to have collected what first call is asking (F = I) calcMode = .false. ! pretend to have collected what first call is asking (F = I)
calcMode(ip,mesh_FEasCP('elem',el)) = .true. ! first ip,el needs to be already pingponged to "calc" calcMode(ip,mesh_FEasCP('elem',el)) = .true. ! first ip,el needs to be already pingponged to "calc"
ip_reshaped = reshape(mesh_ipCoordinates,[3,theMesh%elem%nIPs*theMesh%nElems])
call discretization_init(microstructureAt,homogenizationAt,& call discretization_init(microstructureAt,homogenizationAt,&
reshape(mesh_ipCoordinates,[3,theMesh%elem%nIPs*theMesh%nElems]),& ip_reshaped,&
mesh_node0) mesh_node0)
call results_openJobFile
call HDF5_closeGroup(results_addGroup('geometry'))
call results_writeDataset('geometry',ip_reshaped,'x_c', &
'cell center coordinates','m')
call results_writeDataset('geometry',mesh_node0,'x_n', &
'nodal coordinates','m')
call results_closeJobFile()
call geometry_plastic_nonlocal_setIPvolume(IPvolume()) call geometry_plastic_nonlocal_setIPvolume(IPvolume())
call geometry_plastic_nonlocal_setIPneighborhood(mesh_ipNeighborhood2) call geometry_plastic_nonlocal_setIPneighborhood(mesh_ipNeighborhood2)
call geometry_plastic_nonlocal_setIParea(mesh_ipArea) call geometry_plastic_nonlocal_setIParea(mesh_ipArea)
@ -667,7 +678,6 @@ subroutine mesh_marc_buildElements(microstructureAt,homogenizationAt, &
#if defined(DAMASK_HDF5) #if defined(DAMASK_HDF5)
call results_openJobFile call results_openJobFile
call HDF5_closeGroup(results_addGroup('geometry'))
call results_writeDataset('geometry',mesh_FEnodes,'C',& call results_writeDataset('geometry',mesh_FEnodes,'C',&
'connectivity of the elements','-') 'connectivity of the elements','-')
call results_closeJobFile call results_closeJobFile