diff --git a/code/damage_isoDuctile.f90 b/code/damage_isoDuctile.f90 index e483d64bd..53fb58cd8 100644 --- a/code/damage_isoDuctile.f90 +++ b/code/damage_isoDuctile.f90 @@ -299,7 +299,7 @@ subroutine damage_isoDuctile_microstructure(subdt,ipc, ip, el) damageState(phase)%state(1,constituent) = & damageState(phase)%state(2,constituent) + & (damageState(phase)%subState0(1,constituent) - damageState(phase)%state(2,constituent))* & - exp(-subdt/(damageState(phase)%state(2,constituent)*lattice_DamageMobility(phase))) + exp(-subdt/lattice_DamageMobility(phase)) end subroutine damage_isoDuctile_microstructure