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

prevent use of unallocated array 

interaction matrices are used even if no slip/twin systems defined
This commit is contained in:
Martin Diehl 2018-10-04 07:17:00 +02:00
parent f494803c76
commit cd7f285686
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
src/plastic_phenopowerlaw.f90
View File

@ -234,6 +234,7 @@ subroutine plastic_phenopowerlaw_init
prm%xi_slip_sat = math_expand(prm%xi_slip_sat,prm%Nslip) prm%xi_slip_sat = math_expand(prm%xi_slip_sat,prm%Nslip)
prm%H_int = math_expand(prm%H_int,prm%Nslip) prm%H_int = math_expand(prm%H_int,prm%Nslip)
else slipActive else slipActive
allocate(prm%interaction_SlipSlip(0,0))
allocate(prm%xi_slip_0(0)) allocate(prm%xi_slip_0(0))
endif slipActive endif slipActive
@ -266,6 +267,7 @@ subroutine plastic_phenopowerlaw_init
! expand slip related parameters from system => family ! expand slip related parameters from system => family
prm%xi_twin_0 = math_expand(prm%xi_twin_0,prm%Ntwin) prm%xi_twin_0 = math_expand(prm%xi_twin_0,prm%Ntwin)
else twinActive else twinActive
allocate(prm%interaction_TwinTwin(0,0))
allocate(prm%xi_twin_0(0)) allocate(prm%xi_twin_0(0))
endif twinActive endif twinActive
@ -276,8 +278,9 @@ subroutine plastic_phenopowerlaw_init
prm%interaction_TwinSlip = lattice_interactionTwinSlip2(prm%Ntwin,prm%Nslip,& prm%interaction_TwinSlip = lattice_interactionTwinSlip2(prm%Ntwin,prm%Nslip,&
config_phase(p)%getFloats('interaction_twinslip'), & config_phase(p)%getFloats('interaction_twinslip'), &
structure(1:3)) structure(1:3))
prm%h0_TwinSlip = config_phase(p)%getFloat('h0_twinslip')
else slipAndTwinActive else slipAndTwinActive
allocate(prm%interaction_SlipTwin(prm%totalNslip,prm%TotalNtwin)) ! at least one dimension 0
allocate(prm%interaction_TwinSlip(prm%totalNtwin,prm%TotalNslip)) ! at least one dimension 0
prm%h0_TwinSlip = 0.0_pReal prm%h0_TwinSlip = 0.0_pReal
endif slipAndTwinActive endif slipAndTwinActive