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

fixed buggy state width assignment (endindex += %nSlip corrected to %totalNslip)

This commit is contained in:
Philip Eisenlohr 2018-08-14 01:07:06 +02:00
parent fd578ab8eb
commit 1493b33a03
1 changed files with 17 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

34
src/plastic_phenopowerlaw.f90
View File

@ -440,25 +440,25 @@ subroutine plastic_phenopowerlaw_init
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! locally defined state aliases and initialization of state0 and aTolState ! locally defined state aliases and initialization of state0 and aTolState
startIndex = 1_pInt startIndex = 1_pInt
endIndex = plasticState(p)%nSlip endIndex = prm%totalNslip
state (instance)%s_slip=>plasticState(p)%state (startIndex:endIndex,:) state (instance)%s_slip => plasticState(p)%state (startIndex:endIndex,:)
dotState(instance)%s_slip=>plasticState(p)%dotState(startIndex:endIndex,:) dotState(instance)%s_slip => plasticState(p)%dotState(startIndex:endIndex,:)
plasticState(p)%state0(startIndex:endIndex,:) = & plasticState(p)%state0(startIndex:endIndex,:) = &
spread(math_expand(prm%tau0_slip, prm%Nslip), 2, NipcMyPhase) spread(math_expand(prm%tau0_slip, prm%Nslip), 2, NipcMyPhase)
plasticState(p)%aTolState(startIndex:endIndex) = prm%aTolResistance plasticState(p)%aTolState(startIndex:endIndex) = prm%aTolResistance
startIndex = endIndex + 1_pInt startIndex = endIndex + 1_pInt
endIndex = endIndex + plasticState(p)%nTwin endIndex = endIndex + prm%totalNtwin
state (instance)%s_twin=>plasticState(p)%state (startIndex:endIndex,:) state (instance)%s_twin => plasticState(p)%state (startIndex:endIndex,:)
dotState(instance)%s_twin=>plasticState(p)%dotState(startIndex:endIndex,:) dotState(instance)%s_twin => plasticState(p)%dotState(startIndex:endIndex,:)
plasticState(p)%state0(startIndex:endIndex,:) = & plasticState(p)%state0(startIndex:endIndex,:) = &
spread(math_expand(prm%tau0_twin, prm%Ntwin), 2, NipcMyPhase) spread(math_expand(prm%tau0_twin, prm%Ntwin), 2, NipcMyPhase)
plasticState(p)%aTolState(startIndex:endIndex) = prm%aTolResistance plasticState(p)%aTolState(startIndex:endIndex) = prm%aTolResistance
startIndex = endIndex + 1_pInt startIndex = endIndex + 1_pInt
endIndex = endIndex + 1_pInt endIndex = endIndex + 1_pInt
state (instance)%sumGamma=>plasticState(p)%state (startIndex,:) state (instance)%sumGamma => plasticState(p)%state (startIndex,:)
dotState(instance)%sumGamma=>plasticState(p)%dotState(startIndex,:) dotState(instance)%sumGamma => plasticState(p)%dotState(startIndex,:)
plasticState(p)%aTolState(startIndex:endIndex) = prm%aTolShear plasticState(p)%aTolState(startIndex:endIndex) = prm%aTolShear
startIndex = endIndex + 1_pInt startIndex = endIndex + 1_pInt
@ -468,21 +468,21 @@ subroutine plastic_phenopowerlaw_init
plasticState(p)%aTolState(startIndex:endIndex) = prm%aTolTwinFrac plasticState(p)%aTolState(startIndex:endIndex) = prm%aTolTwinFrac
startIndex = endIndex + 1_pInt startIndex = endIndex + 1_pInt
endIndex = endIndex + plasticState(p)%nSlip endIndex = endIndex + prm%totalNslip
state (instance)%accshear_slip=>plasticState(p)%state (startIndex:endIndex,:) state (instance)%accshear_slip => plasticState(p)%state (startIndex:endIndex,:)
dotState(instance)%accshear_slip=>plasticState(p)%dotState(startIndex:endIndex,:) dotState(instance)%accshear_slip => plasticState(p)%dotState(startIndex:endIndex,:)
plasticState(p)%aTolState(startIndex:endIndex) = prm%aTolShear plasticState(p)%aTolState(startIndex:endIndex) = prm%aTolShear
! global alias ! global alias
plasticState(p)%slipRate =>plasticState(p)%dotState(startIndex:endIndex,:) plasticState(p)%slipRate => plasticState(p)%dotState(startIndex:endIndex,:)
plasticState(p)%accumulatedSlip =>plasticState(p)%state(startIndex:endIndex,:) plasticState(p)%accumulatedSlip => plasticState(p)%state(startIndex:endIndex,:)
startIndex = endIndex + 1_pInt startIndex = endIndex + 1_pInt
endIndex = endIndex + plasticState(p)%nTwin endIndex = endIndex + prm%totalNtwin
state (instance)%accshear_twin=>plasticState(p)%state (startIndex:endIndex,:) state (instance)%accshear_twin => plasticState(p)%state (startIndex:endIndex,:)
dotState(instance)%accshear_twin=>plasticState(p)%dotState(startIndex:endIndex,:) dotState(instance)%accshear_twin => plasticState(p)%dotState(startIndex:endIndex,:)
plasticState(p)%aTolState(startIndex:endIndex) = prm%aTolShear plasticState(p)%aTolState(startIndex:endIndex) = prm%aTolShear
dotState(instance)%whole =>plasticState(p)%dotState dotState(instance)%whole => plasticState(p)%dotState
end associate end associate
enddo enddo