max size of crystallite result is now based on actual use within the mesh, not just picking the largest chunk from the material.config file...

This commit is contained in:
Philip Eisenlohr 2011-06-08 09:16:48 +00:00
parent 7aed307bfc
commit a26d9c844b
3 changed files with 9 additions and 3 deletions

View File

@ -47,7 +47,7 @@ type(p_vec), dimension(:,:,:,:), allocatable :: constitutive_RKCK45dotState
integer(pInt), dimension(:,:,:), allocatable :: constitutive_sizeDotState, & ! size of dotState array
constitutive_sizeState, & ! size of state array per grain
constitutive_sizePostResults ! size of postResults array per grain
integer(pInt) constitutive_maxSizeDotState, &
integer(pInt) constitutive_maxSizeDotState, &
constitutive_maxSizeState, &
constitutive_maxSizePostResults

View File

@ -284,7 +284,13 @@ do i = 1,material_Ncrystallite
endif
enddo
enddo
crystallite_maxSizePostResults = maxval(crystallite_sizePostResults)
crystallite_maxSizePostResults = 0_pInt
do j = 1,material_Nmicrostructure
if (microstructure_active(j)) &
crystallite_maxSizePostResults = max(crystallite_maxSizePostResults,&
crystallite_sizePostResults(microstructure_crystallite(j)))
enddo
! write description file for crystallite output

View File

@ -3352,7 +3352,7 @@ if (mesh_maxValStateVar(2) < 1_pInt) call IO_error(120) ! no microstructure spec
allocate (mesh_HomogMicro(mesh_maxValStateVar(1),mesh_maxValStateVar(2))); mesh_HomogMicro = 0_pInt
do e = 1,mesh_NcpElems
if (mesh_element(3,e) < 1_pInt) call IO_error(110,e) ! no homogenization specified
if (mesh_element(4,e) < 1_pInt) call IO_error(120,e) ! no homogenization specified
if (mesh_element(4,e) < 1_pInt) call IO_error(120,e) ! no microstructure specified
mesh_HomogMicro(mesh_element(3,e),mesh_element(4,e)) = &
mesh_HomogMicro(mesh_element(3,e),mesh_element(4,e)) + 1 ! count combinations of homogenization and microstructure
enddo