From 9320d483052ad345cd9c0eed755b87a67c5e3d07 Mon Sep 17 00:00:00 2001
From: Christoph	Kords <c.kords@mpie.de>
Date: Wed, 27 Oct 2010 08:48:04 +0000
Subject: [PATCH] in RKCK45 state integration: dotState was not saved in last
 Runge-Kutta step

---
 code/crystallite.f90 | 7 ++++---
 1 file changed, 4 insertions(+), 3 deletions(-)

diff --git a/code/crystallite.f90 b/code/crystallite.f90
index f40230d2f..f162cf888 100644
--- a/code/crystallite.f90
+++ b/code/crystallite.f90
@@ -1142,7 +1142,7 @@ endif
 
 
 ! --- FIRST RUNGE KUTTA STEP ---
-
+if (verboseDebugger) write(6,'(a,x,i1)') '<<<RUNGE KUTTA STEP',1
 !$OMP PARALLEL DO
   do e=eIter(1),eIter(2); do i=iIter(1,e),iIter(2,e); do g=gIter(1,e),gIter(2,e)                          ! iterate over elements, ips and grains
     if (crystallite_todo(g,i,e)) then
@@ -1171,7 +1171,6 @@ endif
 
 do n = 1,5
 
-
   ! --- state update ---
   
   !$OMP PARALLEL DO
@@ -1229,7 +1228,7 @@ do n = 1,5
   
 
   ! --- dot state and RK dot state---
-  
+  if (verboseDebugger) write(6,'(a,x,i1)') '<<<RUNGE KUTTA STEP',n+1
   !$OMP PARALLEL DO
     do e=eIter(1),eIter(2); do i=iIter(1,e),iIter(2,e); do g=gIter(1,e),gIter(2,e)                        ! iterate over elements, ips and grains
       selectiveDebugger = (e == debug_e .and. i == debug_i .and. g == debug_g)
@@ -1268,6 +1267,8 @@ relTemperatureResiduum = 0.0_pReal
     if (crystallite_todo(g,i,e)) then
       selectiveDebugger = (e == debug_e .and. i == debug_i .and. g == debug_g)
       sizeDotState = constitutive_sizeDotState(g,i,e)
+      constitutive_RKCK45dotState(6,g,i,e)%p = constitutive_dotState(g,i,e)%p                             ! store Runge-Kutta dotState
+      RKCK45dotTemperature(6,g,i,e) = crystallite_dotTemperature(g,i,e)                                   ! store Runge-Kutta dotTemperature
       
       
       ! --- absolute residuum in state and temperature ---