From a2552d6241543b4587d8e1bb913b474afaf190cd Mon Sep 17 00:00:00 2001 From: Haiming Zhang Date: Sun, 11 Jan 2015 19:10:16 +0000 Subject: [PATCH] Polishing. --- code/crystallite.f90 | 2 +- code/damage_phaseField.f90 | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) 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)