From 975c113ae8c2db3dbed9fb8e7e2c3770d25cf084 Mon Sep 17 00:00:00 2001
From: Luc	Hantcherli <l.hantcherli@mpie.de>
Date: Fri, 7 Dec 2007 10:34:21 +0000
Subject: [PATCH] Homogenised C per grain can now depend on microstructure
 (important for deformation twinning)

---
 trunk/CPFEM.f90 | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/trunk/CPFEM.f90 b/trunk/CPFEM.f90
index 3a6217dfa..9652a096a 100644
--- a/trunk/CPFEM.f90
+++ b/trunk/CPFEM.f90
@@ -403,7 +403,7 @@ if(mod(CPFEM_cn,2)==0) then
     return
  endif
 
- C_66 = constitutive_HomogenizedC(grain, CPFEM_in, cp_en)
+
  A = matmul(Fg_new,invFp_old)  ! actually Fe
  A = matmul(transpose(A), A)
 
@@ -422,6 +422,7 @@ state: do                ! outer iteration: state
            return
          endif
 		 call constitutive_Microstructure(state_new,CPFEM_Temperature(CPFEM_in,cp_en),grain,CPFEM_in,cp_en)
+		 C_66 = constitutive_HomogenizedC(state_new, grain, CPFEM_in, cp_en)
          iStress = 0_pInt
 stress:  do              ! inner iteration: stress
            iStress = iStress+1