From af075aaa9d64372e55c1b60ff450c1e259112291 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Fri, 20 Dec 2013 08:33:46 +0000 Subject: [PATCH] found bug (array index out of bounds), Philip/Franz please check if correct (Gfortran was complaining) --- code/crystallite.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/code/crystallite.f90 b/code/crystallite.f90 index 8d974a205..f1dbca371 100644 --- a/code/crystallite.f90 +++ b/code/crystallite.f90 @@ -2680,7 +2680,7 @@ subroutine crystallite_integrateStateFPI() if (crystallite_todo(g,i,e)) then mySizeDotState = constitutive_sizeDotState(g,i,e) constitutive_state(g,i,e)%p(1:mySizeDotState) = constitutive_subState0(g,i,e)%p(1:mySizeDotState) & - + constitutive_dotState(g,i,e)%p & + + constitutive_dotState(g,i,e)%p(1:mySizeDotState) & * crystallite_subdt(g,i,e) endif enddo; enddo; enddo @@ -2789,8 +2789,8 @@ subroutine crystallite_integrateStateFPI() mySizeDotState = constitutive_sizeDotState(g,i,e) stateResiduum(1:mySizeDotState) = constitutive_state(g,i,e)%p(1:mySizeDotState) & - constitutive_subState0(g,i,e)%p(1:mySizeDotState) & - - (constitutive_dotState(g,i,e)%p * statedamper & - + constitutive_previousDotState(g,i,e)%p & + - (constitutive_dotState(g,i,e)%p(1:mySizeDotState) * statedamper & + + constitutive_previousDotState(g,i,e)%p(1:mySizeDotState) & * (1.0_pReal - statedamper)) * crystallite_subdt(g,i,e) ! --- correct state with residuum ---