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

outputID is part of param(instance) array, reading in of output made consistent with kinematic hardeing law/outputtag not needed

This commit is contained in:
Sharan Roongta 2018-05-29 18:29:38 +02:00
parent 7321a36770
commit 1a5ed20708
1 changed files with 42 additions and 34 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

76
src/plastic_dislotwin.f90
View File

@ -118,10 +118,11 @@ module plastic_dislotwin
strain_trans_fraction_ID, & strain_trans_fraction_ID, &
trans_fraction_ID trans_fraction_ID
end enum end enum
integer(kind(undefined_ID)), dimension(:,:), allocatable, private :: &
outputID !< ID of each post result output
type,private :: tParameters type,private :: tParameters
integer(kind(undefined_ID)), dimension(:), allocatable, private :: &
outputID !< ID of each post result output
real(pReal) :: & real(pReal) :: &
CAtomicVolume, & !< atomic volume in Bugers vector unit CAtomicVolume, & !< atomic volume in Bugers vector unit
D0, & !< prefactor for self-diffusion coefficient D0, & !< prefactor for self-diffusion coefficient
@ -267,7 +268,7 @@ subroutine plastic_dislotwin_init(fileUnit)
Nchunks_SlipTrans = 0_pInt, Nchunks_TransSlip = 0_pInt, Nchunks_TransTrans = 0_pInt, & Nchunks_SlipTrans = 0_pInt, Nchunks_TransSlip = 0_pInt, Nchunks_TransTrans = 0_pInt, &
Nchunks_SlipFamilies = 0_pInt, Nchunks_TwinFamilies = 0_pInt, Nchunks_TransFamilies = 0_pInt, & Nchunks_SlipFamilies = 0_pInt, Nchunks_TwinFamilies = 0_pInt, Nchunks_TransFamilies = 0_pInt, &
offset_slip, index_myFamily, index_otherFamily, & offset_slip, index_myFamily, index_otherFamily, &
startIndex, endIndex startIndex, endIndex, output_ID
integer(pInt) :: sizeState, sizeDotState, sizeDeltaState integer(pInt) :: sizeState, sizeDotState, sizeDeltaState
integer(pInt) :: NofMyPhase integer(pInt) :: NofMyPhase
@ -283,6 +284,9 @@ subroutine plastic_dislotwin_init(fileUnit)
character(len=65536) :: & character(len=65536) :: &
tag = '', & tag = '', &
line = '' line = ''
real(pReal), dimension(:), allocatable :: tempPerSlip, tempPerTwin, tempPerTrans real(pReal), dimension(:), allocatable :: tempPerSlip, tempPerTwin, tempPerTrans
write(6,'(/,a)') ' <<<+- constitutive_'//PLASTICITY_DISLOTWIN_label//' init -+>>>' write(6,'(/,a)') ' <<<+- constitutive_'//PLASTICITY_DISLOTWIN_label//' init -+>>>'
@ -305,7 +309,7 @@ subroutine plastic_dislotwin_init(fileUnit)
allocate(plastic_dislotwin_sizePostResult(maxval(phase_Noutput),maxNinstance),source=0_pInt) allocate(plastic_dislotwin_sizePostResult(maxval(phase_Noutput),maxNinstance),source=0_pInt)
allocate(plastic_dislotwin_output(maxval(phase_Noutput),maxNinstance)) allocate(plastic_dislotwin_output(maxval(phase_Noutput),maxNinstance))
plastic_dislotwin_output = '' plastic_dislotwin_output = ''
allocate(outputID(maxval(phase_Noutput),maxNinstance), source=undefined_ID)
allocate(plastic_dislotwin_Noutput(maxNinstance), source=0_pInt) allocate(plastic_dislotwin_Noutput(maxNinstance), source=0_pInt)
allocate(param(maxNinstance)) allocate(param(maxNinstance))
@ -380,6 +384,7 @@ subroutine plastic_dislotwin_init(fileUnit)
if (IO_getTag(line,'[',']') /= '') then ! next phase section if (IO_getTag(line,'[',']') /= '') then ! next phase section
phase = phase + 1_pInt ! advance phase section counter phase = phase + 1_pInt ! advance phase section counter
if (phase_plasticity(phase) == PLASTICITY_DISLOTWIN_ID) then if (phase_plasticity(phase) == PLASTICITY_DISLOTWIN_ID) then
instance = phase_plasticityInstance(phase)
Nchunks_SlipFamilies = count(lattice_NslipSystem(:,phase) > 0_pInt) Nchunks_SlipFamilies = count(lattice_NslipSystem(:,phase) > 0_pInt)
Nchunks_TwinFamilies = count(lattice_NtwinSystem(:,phase) > 0_pInt) Nchunks_TwinFamilies = count(lattice_NtwinSystem(:,phase) > 0_pInt)
Nchunks_TransFamilies = count(lattice_NtransSystem(:,phase)> 0_pInt) Nchunks_TransFamilies = count(lattice_NtransSystem(:,phase)> 0_pInt)
@ -396,6 +401,7 @@ subroutine plastic_dislotwin_init(fileUnit)
allocate(tempPerSlip(Nchunks_SlipFamilies)) allocate(tempPerSlip(Nchunks_SlipFamilies))
allocate(tempPerTwin(Nchunks_TwinFamilies)) allocate(tempPerTwin(Nchunks_TwinFamilies))
allocate(tempPerTrans(Nchunks_TransFamilies)) allocate(tempPerTrans(Nchunks_TransFamilies))
allocate(param(instance)%outputID(phase_Noutput(phase)), source=undefined_ID)
endif endif
cycle ! skip to next line cycle ! skip to next line
endif endif
@ -406,80 +412,82 @@ subroutine plastic_dislotwin_init(fileUnit)
tag = IO_lc(IO_stringValue(line,chunkPos,1_pInt)) ! extract key tag = IO_lc(IO_stringValue(line,chunkPos,1_pInt)) ! extract key
select case(tag) select case(tag)
case ('(output)') case ('(output)')
plastic_dislotwin_Noutput(instance) = plastic_dislotwin_Noutput(instance) + 1_pInt output_ID = undefined_ID
plastic_dislotwin_output(plastic_dislotwin_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,chunkPos,2_pInt))
select case(IO_lc(IO_stringValue(line,chunkPos,2_pInt))) select case(IO_lc(IO_stringValue(line,chunkPos,2_pInt)))
case ('edge_density') case ('edge_density')
outputID(plastic_dislotwin_Noutput(instance),instance) = edge_density_ID output_ID = edge_density_ID
case ('dipole_density') case ('dipole_density')
outputID(plastic_dislotwin_Noutput(instance),instance) = dipole_density_ID output_ID = dipole_density_ID
case ('shear_rate_slip','shearrate_slip') case ('shear_rate_slip','shearrate_slip')
outputID(plastic_dislotwin_Noutput(instance),instance) = shear_rate_slip_ID output_ID = shear_rate_slip_ID
case ('accumulated_shear_slip') case ('accumulated_shear_slip')
outputID(plastic_dislotwin_Noutput(instance),instance) = accumulated_shear_slip_ID output_ID = accumulated_shear_slip_ID
case ('mfp_slip') case ('mfp_slip')
outputID(plastic_dislotwin_Noutput(instance),instance) = mfp_slip_ID output_ID = mfp_slip_ID
case ('resolved_stress_slip') case ('resolved_stress_slip')
outputID(plastic_dislotwin_Noutput(instance),instance) = resolved_stress_slip_ID output_ID = resolved_stress_slip_ID
case ('threshold_stress_slip') case ('threshold_stress_slip')
outputID(plastic_dislotwin_Noutput(instance),instance) = threshold_stress_slip_ID output_ID= threshold_stress_slip_ID
case ('edge_dipole_distance') case ('edge_dipole_distance')
outputID(plastic_dislotwin_Noutput(instance),instance) = edge_dipole_distance_ID output_ID = edge_dipole_distance_ID
case ('stress_exponent') case ('stress_exponent')
outputID(plastic_dislotwin_Noutput(instance),instance) = stress_exponent_ID output_ID = stress_exponent_ID
case ('twin_fraction') case ('twin_fraction')
outputID(plastic_dislotwin_Noutput(instance),instance) = twin_fraction_ID output_ID = twin_fraction_ID
case ('shear_rate_twin','shearrate_twin') case ('shear_rate_twin','shearrate_twin')
outputID(plastic_dislotwin_Noutput(instance),instance) = shear_rate_twin_ID output_ID= shear_rate_twin_ID
case ('accumulated_shear_twin') case ('accumulated_shear_twin')
outputID(plastic_dislotwin_Noutput(instance),instance) = accumulated_shear_twin_ID output_ID = accumulated_shear_twin_ID
case ('mfp_twin') case ('mfp_twin')
outputID(plastic_dislotwin_Noutput(instance),instance) = mfp_twin_ID output_ID = mfp_twin_ID
case ('resolved_stress_twin') case ('resolved_stress_twin')
outputID(plastic_dislotwin_Noutput(instance),instance) = resolved_stress_twin_ID output_ID = resolved_stress_twin_ID
case ('threshold_stress_twin') case ('threshold_stress_twin')
outputID(plastic_dislotwin_Noutput(instance),instance) = threshold_stress_twin_ID output_ID = threshold_stress_twin_ID
case ('resolved_stress_shearband') case ('resolved_stress_shearband')
outputID(plastic_dislotwin_Noutput(instance),instance) = resolved_stress_shearband_ID output_ID = resolved_stress_shearband_ID
case ('shear_rate_shearband','shearrate_shearband') case ('shear_rate_shearband','shearrate_shearband')
outputID(plastic_dislotwin_Noutput(instance),instance) = shear_rate_shearband_ID output_ID = shear_rate_shearband_ID
case ('sb_eigenvalues') case ('sb_eigenvalues')
outputID(plastic_dislotwin_Noutput(instance),instance) = sb_eigenvalues_ID output_ID = sb_eigenvalues_ID
case ('sb_eigenvectors') case ('sb_eigenvectors')
outputID(plastic_dislotwin_Noutput(instance),instance) = sb_eigenvectors_ID output_ID = sb_eigenvectors_ID
case ('stress_trans_fraction') case ('stress_trans_fraction')
outputID(plastic_dislotwin_Noutput(instance),instance) = stress_trans_fraction_ID output_ID = stress_trans_fraction_ID
case ('strain_trans_fraction') case ('strain_trans_fraction')
outputID(plastic_dislotwin_Noutput(instance),instance) = strain_trans_fraction_ID output_ID = strain_trans_fraction_ID
case ('trans_fraction','total_trans_fraction') case ('trans_fraction','total_trans_fraction')
outputID(plastic_dislotwin_Noutput(instance),instance) = trans_fraction_ID output_ID = trans_fraction_ID
case default
plastic_dislotwin_Noutput(instance) = plastic_dislotwin_Noutput(instance) - 1_pInt
end select end select
if (output_ID /= undefined_ID) then
plastic_dislotwin_Noutput(instance) = plastic_dislotwin_Noutput(instance) + 1_pInt
plastic_dislotwin_output(plastic_dislotwin_Noutput(instance),instance) = IO_lc(IO_stringValue(line,chunkPos,2_pInt))
param(instance)%outputID(plastic_dislotwin_Noutput(instance)) = output_ID
endif
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! parameters depending on number of slip system families ! parameters depending on number of slip system families
case ('nslip') case ('nslip')
@ -826,7 +834,7 @@ subroutine plastic_dislotwin_init(fileUnit)
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! Determine size of postResults array ! Determine size of postResults array
outputsLoop: do o = 1_pInt,plastic_dislotwin_Noutput(instance) outputsLoop: do o = 1_pInt,plastic_dislotwin_Noutput(instance)
select case(outputID(o,instance)) select case(param(instance)%outputID(o))
case(edge_density_ID, & case(edge_density_ID, &
dipole_density_ID, & dipole_density_ID, &
shear_rate_slip_ID, & shear_rate_slip_ID, &
@ -2107,7 +2115,7 @@ function plastic_dislotwin_postResults(Tstar_v,Temperature,ipc,ip,el)
c = 0_pInt c = 0_pInt
plastic_dislotwin_postResults = 0.0_pReal plastic_dislotwin_postResults = 0.0_pReal
do o = 1_pInt,plastic_dislotwin_Noutput(instance) do o = 1_pInt,plastic_dislotwin_Noutput(instance)
select case(outputID(o,instance)) select case(param(instance)%outputID(o))
case (edge_density_ID) case (edge_density_ID)
plastic_dislotwin_postResults(c+1_pInt:c+ns) = state(instance)%rhoEdge(1_pInt:ns,of) plastic_dislotwin_postResults(c+1_pInt:c+ns) = state(instance)%rhoEdge(1_pInt:ns,of)