diff --git a/code/constitutive_dislotwin.f90 b/code/constitutive_dislotwin.f90
index 0ff685ba0..2eb1a11e1 100644
--- a/code/constitutive_dislotwin.f90
+++ b/code/constitutive_dislotwin.f90
@@ -113,8 +113,8 @@ real(pReal),       dimension(:,:,:),       allocatable :: constitutive_dislotwin
                                                           constitutive_dislotwin_interactionMatrixTwinSlip, &   ! interaction matrix of twin systems with slip systems for each instance
                                                           constitutive_dislotwin_interactionMatrixTwinTwin, &   ! interaction matrix of the different twin systems for each instance
                                                           constitutive_dislotwin_forestProjectionEdge           ! matrix of forest projections of edge dislocations for each instance
-real(pReal), dimension(:,:,:,:,:),        allocatable ::  constitutive_dislotwin_sbSv                                                          
-                                                          
+real(pReal), dimension(:,:,:,:,:),        allocatable ::  constitutive_dislotwin_sbSv
+
 CONTAINS
 !****************************************
 !* - constitutive_dislotwin_init
diff --git a/code/crystallite.f90 b/code/crystallite.f90
index 4c31c7353..19caba172 100644
--- a/code/crystallite.f90
+++ b/code/crystallite.f90
@@ -1039,7 +1039,8 @@ do n = 1,4
           constitutive_RK4dotState(g,i,e)%p = constitutive_RK4dotState(g,i,e)%p + weight(n)*constitutive_dotState(g,i,e)%p
           RK4dotTemperature(g,i,e) = RK4dotTemperature(g,i,e) + weight(n)*crystallite_dotTemperature(g,i,e)
         elseif (n == 4) then
-          constitutive_dotState(g,i,e)%p = (constitutive_RK4dotState(g,i,e)%p + weight(n)*constitutive_dotState(g,i,e)%p) /6.0_pReal  ! use weighted RKdotState for final integration
+          constitutive_dotState(g,i,e)%p = (constitutive_RK4dotState(g,i,e)%p + &
+                                            weight(n)*constitutive_dotState(g,i,e)%p) / 6.0_pReal         ! use weighted RKdotState for final integration
           crystallite_dotTemperature(g,i,e) = (RK4dotTemperature(g,i,e) + weight(n)*crystallite_dotTemperature(g,i,e)) / 6.0_pReal
         endif
       endif
diff --git a/code/lattice.f90 b/code/lattice.f90
index 61ce71c59..f3009f070 100644
--- a/code/lattice.f90
+++ b/code/lattice.f90
@@ -776,10 +776,10 @@ subroutine lattice_init()
  allocate(lattice_NslipSystem(lattice_maxNslipFamily,lattice_Nstructure)); lattice_NslipSystem = 0.0_pReal
  allocate(lattice_NtwinSystem(lattice_maxNtwinFamily,lattice_Nstructure)); lattice_NtwinSystem = 0.0_pReal
 
- allocate(lattice_interactionSlipSlip(lattice_maxNslip,lattice_maxNslip,lattice_Nstructure)); lattice_interactionSlipSlip = 0_pInt  ! other:me
- allocate(lattice_interactionSlipTwin(lattice_maxNtwin,lattice_maxNslip,lattice_Nstructure)); lattice_interactionSlipTwin = 0_pInt  ! other:me
- allocate(lattice_interactionTwinSlip(lattice_maxNslip,lattice_maxNtwin,lattice_Nstructure)); lattice_interactionTwinSlip = 0_pInt  ! other:me
- allocate(lattice_interactionTwinTwin(lattice_maxNtwin,lattice_maxNtwin,lattice_Nstructure)); lattice_interactionTwinTwin = 0_pInt  ! other:me
+ allocate(lattice_interactionSlipSlip(lattice_maxNslip,lattice_maxNslip,lattice_Nstructure)); lattice_interactionSlipSlip = 0_pInt ! other:me
+ allocate(lattice_interactionSlipTwin(lattice_maxNtwin,lattice_maxNslip,lattice_Nstructure)); lattice_interactionSlipTwin = 0_pInt ! other:me
+ allocate(lattice_interactionTwinSlip(lattice_maxNslip,lattice_maxNtwin,lattice_Nstructure)); lattice_interactionTwinSlip = 0_pInt ! other:me
+ allocate(lattice_interactionTwinTwin(lattice_maxNtwin,lattice_maxNtwin,lattice_Nstructure)); lattice_interactionTwinTwin = 0_pInt ! other:me
 
 end subroutine