diff --git a/code/lattice.f90 b/code/lattice.f90
index 786f67ebe..a97e124c2 100644
--- a/code/lattice.f90
+++ b/code/lattice.f90
@@ -16,13 +16,13 @@ module lattice
  implicit none
  private
  integer(pInt), parameter, public :: &
-   LATTICE_maxNslipFamily     =  13_pInt, &                                                          !< max # of slip system families over lattice structures
+   LATTICE_maxNslipFamily     =  13_pInt, &                                                         !< max # of slip system families over lattice structures
    LATTICE_maxNtwinFamily     =  4_pInt, &                                                          !< max # of twin system families over lattice structures
    LATTICE_maxNtransFamily    =  2_pInt, &                                                          !< max # of transformation system families over lattice structures
    LATTICE_maxNcleavageFamily =  3_pInt, &                                                          !< max # of transformation system families over lattice structures
    LATTICE_maxNslip           = 52_pInt, &                                                          !< max # of slip systems over lattice structures
    LATTICE_maxNtwin           = 24_pInt, &                                                          !< max # of twin systems over lattice structures
-   LATTICE_maxNinteraction    = 182_pInt, &                                                          !< max # of interaction types (in hardening matrix part)
+   LATTICE_maxNinteraction    = 182_pInt, &                                                         !< max # of interaction types (in hardening matrix part)
    LATTICE_maxNnonSchmid      = 6_pInt, &                                                           !< max # of non schmid contributions over lattice structures
    LATTICE_maxNtrans          = 12_pInt, &                                                          !< max # of transformations over lattice structures
    LATTICE_maxNcleavage       = 9_pInt                                                              !< max # of cleavage over lattice structures
@@ -1628,36 +1628,6 @@ subroutine lattice_initializeStructure(myPhase,CoverA,CoverA_trans,a_fcc,a_bcc)
      call IO_error(135_pInt,el=i,ip=myPhase,ext_msg='matrix diagonal "el"ement of phase "ip"')
  enddo
 
- ! Elasticity matrices for transformed phase
- select case(lattice_structure(myPhase))
-   case (LATTICE_fcc_ID)
-     select case(trans_lattice_structure(myPhase))
-       case (LATTICE_bcc_ID)
-         lattice_trans_C66(1:6,1:6,myPhase) = lattice_C66(1:6,1:6,myPhase)
-         lattice_trans_mu(myPhase) = lattice_mu(myPhase)
-         lattice_trans_nu(myPhase) = lattice_nu(myPhase)
-         lattice_trans_C3333(1:3,1:3,1:3,1:3,myPhase) = lattice_C3333(1:3,1:3,1:3,1:3,myPhase)
-       case (LATTICE_hex_ID)
-         lattice_trans_C66(1:6,1:6,myPhase) = lattice_symmetrizeC66(trans_lattice_structure(myPhase),&
-                                                                   lattice_trans_C66(1:6,1:6,myPhase))
-         lattice_trans_mu(myPhase) = 0.2_pReal *(  lattice_trans_C66(1,1,myPhase) &
-                                                 - lattice_trans_C66(1,2,myPhase) &
-                                                 + 3.0_pReal*lattice_trans_C66(4,4,myPhase))
-         lattice_trans_nu(myPhase) = (  lattice_trans_C66(1,1,myPhase) &
-                                      + 4.0_pReal*lattice_trans_C66(1,2,myPhase) &
-                                      - 2.0_pReal*lattice_trans_C66(4,4,myPhase)) &
-                                                     /(  4.0_pReal*lattice_trans_C66(1,1,myPhase) &
-                                                       + 6.0_pReal*lattice_trans_C66(1,2,myPhase) &
-                                                       + 2.0_pReal*lattice_trans_C66(4,4,myPhase))
-         lattice_trans_C3333(1:3,1:3,1:3,1:3,myPhase) = math_Voigt66to3333(lattice_trans_C66(1:6,1:6,myPhase))
-         lattice_trans_C66(1:6,1:6,myPhase) = math_Mandel3333to66(lattice_trans_C3333(1:3,1:3,1:3,1:3,myPhase))
-         do i = 1_pInt, 6_pInt
-           if (abs(lattice_trans_C66(i,i,myPhase))<tol_math_check) &
-           call IO_error(136_pInt,el=i,ip=myPhase,ext_msg='matrix diagonal "el"ement of phase "ip"')
-         enddo
-     end select
- end select
-
  lattice_thermalConductivity33(1:3,1:3,myPhase) = lattice_symmetrize33(lattice_structure(myPhase),&
                                                                        lattice_thermalConductivity33(1:3,1:3,myPhase))
  lattice_thermalExpansion33(1:3,1:3,myPhase) = lattice_symmetrize33(lattice_structure(myPhase),&
diff --git a/code/plastic_dislotwin.f90 b/code/plastic_dislotwin.f90
index 23d592364..6497b3d8e 100644
--- a/code/plastic_dislotwin.f90
+++ b/code/plastic_dislotwin.f90
@@ -1033,7 +1033,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_C3333(p,q,r,s,instance) * &
                lattice_Qtrans(l,p,index_otherFamily+j,phase) * &
                lattice_Qtrans(m,q,index_otherFamily+j,phase) * &
                lattice_Qtrans(n,r,index_otherFamily+j,phase) * &