diff --git a/code/crystallite.f90 b/code/crystallite.f90 index a49a9b890..1ca40bb52 100644 --- a/code/crystallite.f90 +++ b/code/crystallite.f90 @@ -3271,7 +3271,7 @@ subroutine crystallite_integrateStateFPI() - damageState(p)%subState0(1:mySizeDamageDotState,c) & - (damageState(p)%dotState(1:mySizeDamageDotState,c) * damageStateDamper & + damageState(p)%previousDotState(1:mySizeDamageDotState,c) & - * (1.0_pReal - damageStatedamper)) * crystallite_subdt(g,i,e) + * (1.0_pReal - damageStateDamper)) * crystallite_subdt(g,i,e) thermalStateResiduum(1:mySizeThermalDotState) = thermalState(p)%state(1:mySizeThermalDotState,c) & - thermalState(p)%subState0(1:mySizeThermalDotState,c) & diff --git a/code/damage_phaseField.f90 b/code/damage_phaseField.f90 index d7a565a56..a7cdb4753 100644 --- a/code/damage_phaseField.f90 +++ b/code/damage_phaseField.f90 @@ -214,7 +214,7 @@ subroutine damage_phaseField_init(fileUnit) allocate(damageState(phase)%state_backup (sizeState,NofMyPhase), source=0.0_pReal) allocate(damageState(phase)%dotState (sizeDotState,NofMyPhase), source=0.0_pReal) - allocate(damageState(phase)%deltaState (sizeDotState,NofMyPhase), source=0.0_pReal) + allocate(damageState(phase)%deltaState (sizeDotState,NofMyPhase), source=0.0_pReal) allocate(damageState(phase)%dotState_backup (sizeDotState,NofMyPhase), source=0.0_pReal) if (any(numerics_integrator == 1_pInt)) then allocate(damageState(phase)%previousDotState (sizeDotState,NofMyPhase), source=0.0_pReal)