diff --git a/code/crystallite.f90 b/code/crystallite.f90
index 40fa0fec2..d87d0407f 100644
--- a/code/crystallite.f90
+++ b/code/crystallite.f90
@@ -323,6 +323,7 @@ close(file)
         crystallite_localConstitution(g,i,e) = phase_localConstitution(material_phase(g,i,e))
         !$OMP FLUSH(crystallite_Fp0)
         crystallite_Fe(1:3,1:3,g,i,e)  = math_transpose33(crystallite_Fp0(1:3,1:3,g,i,e))
+        crystallite_Fp(1:3,1:3,g,i,e)  = crystallite_Fp0(1:3,1:3,g,i,e)
       enddo
     enddo
   enddo