diff --git a/src/plastic_dislotwin.f90 b/src/plastic_dislotwin.f90
index 50b14bdf9..d6c73e8f3 100644
--- a/src/plastic_dislotwin.f90
+++ b/src/plastic_dislotwin.f90
@@ -1029,7 +1029,7 @@ subroutine plastic_dislotwin_init(fileUnit)
            do p = 1_pInt,3_pInt; do q = 1_pInt,3_pInt; do r = 1_pInt,3_pInt; do s = 1_pInt,3_pInt
              plastic_dislotwin_Ctwin3333(l,m,n,o,index_myFamily+j,instance) = &
              plastic_dislotwin_Ctwin3333(l,m,n,o,index_myFamily+j,instance) + &
-               lattice_C3333(p,q,r,s,instance) * &
+               lattice_C3333(p,q,r,s,phase) * &
                lattice_Qtwin(l,p,index_otherFamily+j,phase) * &
                lattice_Qtwin(m,q,index_otherFamily+j,phase) * &
                lattice_Qtwin(n,r,index_otherFamily+j,phase) * &
@@ -1087,7 +1087,7 @@ subroutine plastic_dislotwin_init(fileUnit)
            do p = 1_pInt,3_pInt; do q = 1_pInt,3_pInt; do r = 1_pInt,3_pInt; do s = 1_pInt,3_pInt
              plastic_dislotwin_Ctrans3333(l,m,n,o,index_myFamily+j,instance) = &
              plastic_dislotwin_Ctrans3333(l,m,n,o,index_myFamily+j,instance) + &
-               lattice_trans_C3333(p,q,r,s,instance) * &
+               lattice_trans_C3333(p,q,r,s,phase) * &
                lattice_Qtrans(l,p,index_otherFamily+j,phase) * &
                lattice_Qtrans(m,q,index_otherFamily+j,phase) * &
                lattice_Qtrans(n,r,index_otherFamily+j,phase) * &