diff --git a/src/plastic_disloUCLA.f90 b/src/plastic_disloUCLA.f90
index 36c785b24..44cc2bac0 100644
--- a/src/plastic_disloUCLA.f90
+++ b/src/plastic_disloUCLA.f90
@@ -198,7 +198,8 @@ subroutine plastic_disloUCLA_init()
                                                        config%getString('lattice_structure'))
       prm%forestProjectionEdge = lattice_forestProjection_edge(prm%N_sl,config%getString('lattice_structure'),&
                                                                config%getFloat('c/a',defaultVal=0.0_pReal))
- 
+      prm%forestProjectionEdge = transpose(prm%forestProjectionEdge)
+      
       prm%rho_mob_0   = config%getFloats('rhoedge0',       requiredSize=size(prm%N_sl))
       prm%rho_dip_0   = config%getFloats('rhoedgedip0',    requiredSize=size(prm%N_sl))
       prm%v0          = config%getFloats('v0',             requiredSize=size(prm%N_sl))
@@ -458,7 +459,7 @@ subroutine plastic_disloUCLA_dependentState(instance,of)
  
   associate(prm => param(instance), stt => state(instance),dst => dependentState(instance))
  
-  dislocationSpacing = sqrt(matmul(transpose(prm%forestProjectionEdge), &
+  dislocationSpacing = sqrt(matmul(prm%forestProjectionEdge, &
                                    stt%rho_mob(:,of)+stt%rho_dip(:,of)))
   dst%threshold_stress(:,of) = prm%mu*prm%b_sl &
                              * sqrt(matmul(prm%h_sl_sl,stt%rho_mob(:,of)+stt%rho_dip(:,of)))
diff --git a/src/plastic_dislotwin.f90 b/src/plastic_dislotwin.f90
index 09d7737a5..4e085478f 100644
--- a/src/plastic_dislotwin.f90
+++ b/src/plastic_dislotwin.f90
@@ -253,6 +253,7 @@ subroutine plastic_dislotwin_init
                                                   config%getString('lattice_structure'))
      prm%forestProjection     = lattice_forestProjection_edge(prm%N_sl,config%getString('lattice_structure'),&
                                                               config%getFloat('c/a',defaultVal=0.0_pReal))
+     prm%forestProjection     = transpose(prm%forestProjection)
 
      prm%n0_sl                = lattice_slip_normal(prm%N_sl,config%getString('lattice_structure'),&
                                                     config%getFloat('c/a',defaultVal=0.0_pReal))  
@@ -884,7 +885,7 @@ subroutine plastic_dislotwin_dependentState(T,instance,of)
  f_over_t_tr = sumf_trans/prm%t_tr                   ! but this not
                                                      ! ToDo ...Physically correct, but naming could be adjusted
 
- inv_lambda_sl_sl = sqrt(matmul(transpose(prm%forestProjection), &
+ inv_lambda_sl_sl = sqrt(matmul(prm%forestProjection, &
                                 stt%rho_mob(:,of)+stt%rho_dip(:,of)))/prm%CLambdaSlip
 
  if (prm%sum_N_tw > 0 .and. prm%sum_N_sl > 0) &