diff --git a/src/homogenization_RGC.f90 b/src/homogenization_RGC.f90
index 92ea5301d..1d7bc6f86 100644
--- a/src/homogenization_RGC.f90
+++ b/src/homogenization_RGC.f90
@@ -942,10 +942,10 @@ pure function homogenization_RGC_postResults(ip,el,avgP,avgF)
  do o = 1_pInt,homogenization_Noutput(mesh_element(3,el))
    select case(homogenization_RGC_outputID(o,homID))
      case (avgdefgrad_ID)
-       homogenization_RGC_postResults(c+1_pInt:c+9_pInt) = reshape(avgF,[9])
+       homogenization_RGC_postResults(c+1_pInt:c+9_pInt) = reshape(transpose(avgF),[9])
        c = c + 9_pInt
      case (avgfirstpiola_ID)
-       homogenization_RGC_postResults(c+1_pInt:c+9_pInt) = reshape(avgP,[9])
+       homogenization_RGC_postResults(c+1_pInt:c+9_pInt) = reshape(transpose(avgP),[9])
        c = c + 9_pInt
      case (ipcoords_ID)
        homogenization_RGC_postResults(c+1_pInt:c+3_pInt) = mesh_ipCoordinates(1:3,ip,el)             ! current ip coordinates
diff --git a/src/homogenization_isostrain.f90 b/src/homogenization_isostrain.f90
index 8ee0df73d..24aedf75f 100644
--- a/src/homogenization_isostrain.f90
+++ b/src/homogenization_isostrain.f90
@@ -300,10 +300,10 @@ pure function homogenization_isostrain_postResults(ip,el,avgP,avgF)
        homogenization_isostrain_postResults(c+1_pInt) = real(homogenization_isostrain_Ngrains(homID),pReal)
        c = c + 1_pInt
      case (avgdefgrad_ID)
-       homogenization_isostrain_postResults(c+1_pInt:c+9_pInt) = reshape(avgF,[9])
+       homogenization_isostrain_postResults(c+1_pInt:c+9_pInt) = reshape(transpose(avgF),[9])
        c = c + 9_pInt
      case (avgfirstpiola_ID)
-       homogenization_isostrain_postResults(c+1_pInt:c+9_pInt) = reshape(avgP,[9])
+       homogenization_isostrain_postResults(c+1_pInt:c+9_pInt) = reshape(transpose(avgP),[9])
        c = c + 9_pInt
      case (ipcoords_ID)
        homogenization_isostrain_postResults(c+1_pInt:c+3_pInt) = mesh_ipCoordinates(1:3,ip,el)                       ! current ip coordinates