diff --git a/PRIVATE b/PRIVATE
index 1e8c66897..0ce46aeb2 160000
--- a/PRIVATE
+++ b/PRIVATE
@@ -1 +1 @@
-Subproject commit 1e8c66897820468ab46958d995005e2b69204d0e
+Subproject commit 0ce46aeb2dc74566642a4076b4218f809977e503
diff --git a/examples/SpectralMethod/Polycrystal/material.yaml b/examples/SpectralMethod/Polycrystal/material.yaml
index 0717d4eab..04775be4c 100644
--- a/examples/SpectralMethod/Polycrystal/material.yaml
+++ b/examples/SpectralMethod/Polycrystal/material.yaml
@@ -110,7 +110,7 @@ phase:
   Aluminum:
     elasticity: {C_11: 106.75e9, C_12: 60.41e9, C_44: 28.34e9, type: hooke}
     generic:
-      output: [F, P, Fe, Fp, Lp, O]
+      output: [F, P, F_e, F_p, L_p, O]
     lattice: fcc
     plasticity:
       N_sl: [12]
diff --git a/src/crystallite.f90 b/src/crystallite.f90
index c00b076a3..f2b305aa2 100644
--- a/src/crystallite.f90
+++ b/src/crystallite.f90
@@ -759,23 +759,23 @@ subroutine crystallite_results
           selected_tensors = select_tensors(crystallite_partitionedF,p)
           call results_writeDataset(group,selected_tensors,output_constituent(p)%label(o),&
                                    'deformation gradient','1')
-        case('Fe')
+        case('F_e')
           selected_tensors = select_tensors(crystallite_Fe,p)
           call results_writeDataset(group,selected_tensors,output_constituent(p)%label(o),&
                                    'elastic deformation gradient','1')
-        case('Fp')
+        case('F_p')
           selected_tensors = select_tensors(crystallite_Fp,p)
           call results_writeDataset(group,selected_tensors,output_constituent(p)%label(o),&
                                    'plastic deformation gradient','1')
-        case('Fi')
+        case('F_i')
           selected_tensors = select_tensors(crystallite_Fi,p)
           call results_writeDataset(group,selected_tensors,output_constituent(p)%label(o),&
                                    'inelastic deformation gradient','1')
-        case('Lp')
+        case('L_p')
           selected_tensors = select_tensors(crystallite_Lp,p)
           call results_writeDataset(group,selected_tensors,output_constituent(p)%label(o),&
                                    'plastic velocity gradient','1/s')
-        case('Li')
+        case('L_i')
           selected_tensors = select_tensors(crystallite_Li,p)
           call results_writeDataset(group,selected_tensors,output_constituent(p)%label(o),&
                                    'inelastic velocity gradient','1/s')