diff --git a/code/lattice.f90 b/code/lattice.f90
index 5060b4b26..3e353e7ef 100644
--- a/code/lattice.f90
+++ b/code/lattice.f90
@@ -433,33 +433,33 @@ integer(pInt), allocatable, dimension(:,:,:) :: lattice_interactionSlipSlip, &
  real(pReal), dimension(4+4,lattice_hex_Nslip), parameter :: lattice_hex_systemSlip = &
  reshape((/&
 ! Basal systems <1120>{0001} (independent of c/a-ratio, Bravais notation (4 coordinate base))
-  2, -1, -1,  0,     0,  0,  0,  1, &
- -1,  2, -1,  0,     0,  0,  0,  1, &
- -1, -1,  2,  0,     0,  0,  0,  1, &
+  2, -1, -1,  0,     0,  0,  0,  1, & !A1
+ -1,  2, -1,  0,     0,  0,  0,  1, & !A2
+ -1, -1,  2,  0,     0,  0,  0,  1, & !A3
 ! 1st type prismatic systems <1120>{1010}  (independent of c/a-ratio)
-  2, -1, -1,  0,     0,  1, -1,  0, &
- -1,  2, -1,  0,    -1,  0,  1,  0, &
- -1, -1,  2,  0,     1, -1,  0,  0, &
+  2, -1, -1,  0,     0,  1, -1,  0, & !B1
+ -1,  2, -1,  0,    -1,  0,  1,  0, & !C2
+ -1, -1,  2,  0,     1, -1,  0,  0, & !D3
 ! 1st type 1st order pyramidal systems <1120>{1011} -- plane normals depend on the c/a-ratio
-  2, -1, -1,  0,     0,  1, -1,  1, &
-  1,  1, -2,  0,    -1,  1,  0,  1, &
- -1,  2, -1,  0,    -1,  0,  1,  1, &
- -2,  1,  1,  0,     0, -1,  1,  1, &
- -1, -1,  2,  0,     1, -1,  0,  1, &
-  1, -2,  1,  0,     1,  0, -1,  1, &
+  2, -1, -1,  0,     0,  1, -1,  1, & !E1
+  1,  1, -2,  0,    -1,  1,  0,  1, & !F(-3)
+ -1,  2, -1,  0,    -1,  0,  1,  1, & !G2
+ -2,  1,  1,  0,     0, -1,  1,  1, & !H(-1)
+ -1, -1,  2,  0,     1, -1,  0,  1, & !I3
+  1, -2,  1,  0,     1,  0, -1,  1, & !J(-2)
 ! pyramidal system: c+a slip <2113>{1011} -- plane normals depend on the c/a-ratio
- -1,  2, -1,  3,     0,  1, -1,  1, &
-  1,  1, -2,  3,     0,  1, -1,  1, &
- -2,  1,  1,  3,    -1,  1,  0,  1, &
- -1,  2, -1,  3,    -1,  1,  0,  1, &
- -1, -1,  2,  3,    -1,  0,  1,  1, &
- -2,  1,  1,  3,    -1,  0,  1,  1, &
-  1, -2,  1,  3,     0, -1,  1,  1, &
- -1, -1,  2,  3,     0, -1,  1,  1, &
-  2, -1, -1,  3,     1, -1,  0,  1, &
-  1, -2,  1,  3,     1, -1,  0,  1, &
-  1,  1, -2,  3,     1,  0, -1,  1, &
-  2, -1, -1,  3,     1,  0, -1,  1, &
+ -1,  2, -1,  3,     0,  1, -1,  1, & !E4
+  1,  1, -2,  3,     0,  1, -1,  1, & !E5
+ -2,  1,  1,  3,    -1,  1,  0,  1, & !F6
+ -1,  2, -1,  3,    -1,  1,  0,  1, & !F4
+ -1, -1,  2,  3,    -1,  0,  1,  1, & !G7
+ -2,  1,  1,  3,    -1,  0,  1,  1, & !G6
+  1, -2,  1,  3,     0, -1,  1,  1, & !H8
+ -1, -1,  2,  3,     0, -1,  1,  1, & !H7
+  2, -1, -1,  3,     1, -1,  0,  1, & !I9
+  1, -2,  1,  3,     1, -1,  0,  1, & !I8
+  1,  1, -2,  3,     1,  0, -1,  1, & !J5
+  2, -1, -1,  3,     1,  0, -1,  1 &  !J9
    /),(/4+4,lattice_hex_Nslip/))
 
  real(pReal), dimension(4+4,lattice_hex_Ntwin), parameter :: lattice_hex_systemTwin = &