diff --git a/code/homogenization.f90 b/code/homogenization.f90
index 71e1e1b32..1e88f063e 100644
--- a/code/homogenization.f90
+++ b/code/homogenization.f90
@@ -662,11 +662,12 @@ function homogenization_postResults(&
    ip, &            ! integration point
    el  &            ! element
   )
- use prec,     only: pReal,pInt
- use mesh,     only: mesh_element
- use material, only: homogenization_type
+ use prec,        only: pReal,pInt
+ use mesh,        only: mesh_element
+ use material,    only: homogenization_type
+ use crystallite, only: crystallite_partionedF
  use homogenization_isostrain
- use homogenization_RGC             ! RGC homogenization added <<<updated 31.07.2009>>>
+ use homogenization_RGC             ! RGC homogenization added <<<updated 22.10.2009>>>
  implicit none
 
 !* Definition of variables
@@ -678,9 +679,10 @@ function homogenization_postResults(&
 !* isostrain
    case (homogenization_isostrain_label)
      homogenization_postResults = homogenization_isostrain_postResults(homogenization_state(ip,el),ip,el)
-!* RGC homogenization added <<<updated 31.07.2009>>>
+!* RGC homogenization added <<<updated 22.10.2009>>>
    case (homogenization_RGC_label)
-     homogenization_postResults = homogenization_RGC_postResults(homogenization_state(ip,el),ip,el)
+     homogenization_postResults = homogenization_RGC_postResults(crystallite_partionedF(:,:,:,ip,el), &
+                                                                 homogenization_state(ip,el),ip,el)
  end select
 
  return
diff --git a/code/homogenization_RGC.f90 b/code/homogenization_RGC.f90
index d4c86e6b7..7a55eb400 100644
--- a/code/homogenization_RGC.f90
+++ b/code/homogenization_RGC.f90
@@ -125,6 +125,11 @@ subroutine homogenization_RGC_init(&
        case('magnitudemismatch')
          homogenization_RGC_sizePostResults(i) = &
          homogenization_RGC_sizePostResults(i) + 1
+       case('grainsdeformation')
+         homogenization_RGC_sizePostResults(i) = &
+         homogenization_RGC_sizePostResults(i) + 9_pInt*(homogenization_RGC_Ngrains(1,i)* &
+                                                         homogenization_RGC_Ngrains(2,i)* &
+                                                         homogenization_RGC_Ngrains(3,i))
      end select
    enddo
 
@@ -582,7 +587,7 @@ subroutine homogenization_RGC_averageStressAndItsTangent(&
 
  use prec, only: pReal,pInt,p_vec
  use mesh, only: mesh_element,mesh_NcpElems,mesh_maxNips
- use material, only: homogenization_maxNgrains, homogenization_Ngrains,homogenization_typeInstance
+ use material, only: homogenization_maxNgrains,homogenization_Ngrains,homogenization_typeInstance
  use math, only: math_Plain3333to99
  implicit none
 
@@ -652,6 +657,7 @@ endfunction
 ! return array of homogenization results for post file inclusion
 !********************************************************************
 pure function homogenization_RGC_postResults(&
+   F, &             ! array of current grain deformation gradients
    state, &         ! my state
    ip, &            ! my integration point
    el  &            ! my element
@@ -659,14 +665,16 @@ pure function homogenization_RGC_postResults(&
 
  use prec, only: pReal,pInt,p_vec
  use mesh, only: mesh_element
- use material, only: homogenization_typeInstance,homogenization_Noutput,homogenization_Ngrains
+ use material, only: homogenization_typeInstance,homogenization_Noutput,homogenization_Ngrains, &
+                     homogenization_maxNgrains
  implicit none
 
 !* Definition of variables
+ real(pReal), dimension (3,3,homogenization_maxNgrains), intent(in) :: F
  type(p_vec), intent(in) :: state
  integer(pInt), intent(in) :: ip,el
 !
- integer(pInt) homID,o,c,nIntFaceTot
+ integer(pInt) homID,o,c,nIntFaceTot,i,j,k
  real(pReal), dimension(homogenization_RGC_sizePostResults(homogenization_typeInstance(mesh_element(3,el)))) :: &
    homogenization_RGC_postResults
 
@@ -688,6 +696,15 @@ pure function homogenization_RGC_postResults(&
      case('magnitudemismatch')
        homogenization_RGC_postResults(c+1) = state%p(3*nIntFaceTot+3)
        c = c + 1
+     case('grainsdeformation')
+       do k = 1,homogenization_RGC_Ngrains(1,homID)*homogenization_RGC_Ngrains(2,homID)*homogenization_RGC_Ngrains(3,homID)
+         do i = 1,3
+         do j = 1,3
+           homogenization_RGC_postResults(c+1) = F(i,j,k)
+           c = c + 1
+         enddo
+         enddo
+       enddo
    end select
  enddo
  
diff --git a/code/material.config b/code/material.config
index 2eb817490..fa2dd1790 100644
--- a/code/material.config
+++ b/code/material.config
@@ -11,12 +11,13 @@ Ngrains	1
 [Grain_Cluster]
 type	RGC
 ngrains			8
-clustersize		2	2	2        # product of these numbers must be equal to ngrains(!)
-grainsizeparameter	3.2e+3	3.2e+3	3.2e+3   # 	3.9e3	3.6e3	3.9e3
-overproportionality	1.6e+0	1.6e+0	1.6e+0	 #	1.6e0	1.6e0	1.6e0
+clustersize		2	2	2                    # product of these numbers must be equal to ngrains(!)
+grainsizeparameter	3.2e+3	3.2e+3	3.2e+3
+overproportionality	1.6e+0	1.6e+0	1.6e+0
 (output)		constitutivework
 (output)		penaltyenergy
 (output)		magnitudemismatch
+(output)		grainsdeformation            # warning: tensorial output, size = ngrains x 3 x 3 components
 
 [Taylor2]
 type	isostrain