diff --git a/code/crystallite.f90 b/code/crystallite.f90
index 230bfcc20..5f1e9e365 100644
--- a/code/crystallite.f90
+++ b/code/crystallite.f90
@@ -261,7 +261,7 @@ do i = 1,material_Ncrystallite
         mySize = 3
       case('grainrotation') ! Deviation from initial grain orientation in axis-angle form (angle in degrees)
         mySize = 4
-      case('defgrad','f','fe','fp','ee','p','firstpiola','1stpiola','s','tstar','secondpiola','2ndpiola')
+      case('defgrad','f','fe','fp','lp','ee','p','firstpiola','1stpiola','s','tstar','secondpiola','2ndpiola')
         mySize = 9
       case default
         mySize = 0      
@@ -2894,6 +2894,10 @@ function crystallite_postResults(&
        mySize = 9_pInt
        crystallite_postResults(c+1:c+1+mySize) = reshape(math_transpose3x3(crystallite_Fp(:,:,g,i,e)),(/mySize/))
        c = c + mySize
+     case ('lp')
+       mySize = 9_pInt
+       crystallite_postResults(c+1:c+1+mySize) = reshape(math_transpose3x3(crystallite_Lp(:,:,g,i,e)),(/mySize/))
+       c = c + mySize
      case ('p','firstpiola','1stpiola')
        mySize = 9_pInt
        crystallite_postResults(c+1:c+1+mySize) = reshape(math_transpose3x3(crystallite_P(:,:,g,i,e)),(/mySize/))
diff --git a/code/material.config b/code/material.config
index a20be0e2e..92e316f48 100644
--- a/code/material.config
+++ b/code/material.config
@@ -73,6 +73,7 @@ crystallite 1
 (output) ee             # elastic strain as Green-Lagrange tensor 
 (output) p              # first Piola-Kichhoff stress tensor; synonyms: "firstpiola", "1stpiola"
 (output) s              # second Piola-Kichhoff stress tensor; synonyms: "tstar", "secondpiola", "2ndpiola"
+(output) lp             # plastic velocity gradient tensor
 
 #-------------------#
 <phase>