diff --git a/PRIVATE b/PRIVATE
index 15bb4a182..70310b426 160000
--- a/PRIVATE
+++ b/PRIVATE
@@ -1 +1 @@
-Subproject commit 15bb4a18292bf4f5d3b371dc1edbf9a33606eb60
+Subproject commit 70310b426d03ef07dc053e080729d35e5eb3f0cb
diff --git a/examples/config/phase/damage/anisobrittle_cubic.yaml b/examples/config/phase/damage/anisobrittle_cubic.yaml
index 714be6195..5320a4aa9 100644
--- a/examples/config/phase/damage/anisobrittle_cubic.yaml
+++ b/examples/config/phase/damage/anisobrittle_cubic.yaml
@@ -7,5 +7,5 @@ q: 20
 
 output: [f_phi]
 
-D_11: 1.0
-M: 0.001
+K_11: 1.0
+mu: 0.001
diff --git a/examples/config/phase/damage/isobrittle_generic.yaml b/examples/config/phase/damage/isobrittle_generic.yaml
index b4416f995..98c9b9763 100644
--- a/examples/config/phase/damage/isobrittle_generic.yaml
+++ b/examples/config/phase/damage/isobrittle_generic.yaml
@@ -5,5 +5,5 @@ isoBrittle_atol: 0.01
 
 output: [f_phi]
 
-D_11: 1.0
-M: 0.001
+K_11: 1.0
+mu: 0.001
diff --git a/src/grid/grid_damage_spectral.f90 b/src/grid/grid_damage_spectral.f90
index 26e0a909e..7851f4f50 100644
--- a/src/grid/grid_damage_spectral.f90
+++ b/src/grid/grid_damage_spectral.f90
@@ -188,9 +188,9 @@ function grid_damage_spectral_solution(timeinc) result(solution)
   stagNorm = maxval(abs(phi_current - phi_stagInc))
   solnNorm = maxval(abs(phi_current))
   call MPI_Allreduce(MPI_IN_PLACE,stagNorm,1,MPI_DOUBLE,MPI_MAX,PETSC_COMM_WORLD,ierr)
-  call MPI_Allreduce(MPI_IN_PLACE,solnNorm,1,MPI_DOUBLE,MPI_MAX,PETSC_COMM_WORLD,ierr)
-  phi_stagInc = phi_current
   solution%stagConverged = stagNorm < max(num%eps_damage_atol, num%eps_damage_rtol*solnNorm)
+  call MPI_Allreduce(MPI_IN_PLACE,solution%stagConverged,1,MPI_LOGICAL,MPI_LAND,PETSC_COMM_WORLD,ierr)
+  phi_stagInc = phi_current
 
 !--------------------------------------------------------------------------------------------------
 ! updating damage state
diff --git a/src/grid/grid_thermal_spectral.f90 b/src/grid/grid_thermal_spectral.f90
index 5a1745668..6af80825b 100644
--- a/src/grid/grid_thermal_spectral.f90
+++ b/src/grid/grid_thermal_spectral.f90
@@ -183,9 +183,9 @@ function grid_thermal_spectral_solution(timeinc) result(solution)
   stagNorm = maxval(abs(T_current - T_stagInc))
   solnNorm = maxval(abs(T_current))
   call MPI_Allreduce(MPI_IN_PLACE,stagNorm,1,MPI_DOUBLE,MPI_MAX,PETSC_COMM_WORLD,ierr)
-  call MPI_Allreduce(MPI_IN_PLACE,solnNorm,1,MPI_DOUBLE,MPI_MAX,PETSC_COMM_WORLD,ierr)
-  T_stagInc = T_current
   solution%stagConverged = stagNorm < max(num%eps_thermal_atol, num%eps_thermal_rtol*solnNorm)
+  call MPI_Allreduce(MPI_IN_PLACE,solution%stagConverged,1,MPI_LOGICAL,MPI_LAND,PETSC_COMM_WORLD,ierr)
+  T_stagInc = T_current
 
 !--------------------------------------------------------------------------------------------------
 ! updating thermal state
diff --git a/src/phase_damage.f90 b/src/phase_damage.f90
index e9b62b702..d8ac5046a 100644
--- a/src/phase_damage.f90
+++ b/src/phase_damage.f90
@@ -74,7 +74,7 @@ contains
 module subroutine damage_init
 
   integer :: &
-    ph, &                                                                                           !< counter in phase loop
+    ph, &
     Nmembers
   class(tNode), pointer :: &
    phases, &
@@ -105,10 +105,10 @@ module subroutine damage_init
     if (sources%length == 1) then
       damage_active = .true.
       source => sources%get(1)
-      param(ph)%mu     = source%get_asFloat('M',defaultVal=0.0_pReal)
-      param(ph)%K(1,1) = source%get_asFloat('D_11',defaultVal=0.0_pReal)
-      param(ph)%K(3,3) = source%get_asFloat('D_33',defaultVal=0.0_pReal)
-      param(ph)%K = lattice_applyLatticeSymmetry33(param(ph)%K,phase%get_asString('lattice'))
+      param(ph)%mu     = source%get_asFloat('mu',defaultVal=0.0_pReal)                              ! ToDo: make mandatory?
+      param(ph)%K(1,1) = source%get_asFloat('K_11',defaultVal=0.0_pReal)                            ! ToDo: make mandatory?
+      param(ph)%K(3,3) = source%get_asFloat('K_33',defaultVal=0.0_pReal)                            ! ToDo: depends on symmetry
+      param(ph)%K = lattice_applyLatticeSymmetry33(param(ph)%K,phase_lattice(ph))
     endif
 
   enddo
diff --git a/src/phase_thermal.f90 b/src/phase_thermal.f90
index 808cd7cf2..8d2915ef7 100644
--- a/src/phase_thermal.f90
+++ b/src/phase_thermal.f90
@@ -100,10 +100,10 @@ module subroutine thermal_init(phases)
     allocate(current(ph)%dot_T(Nmembers),source=0.0_pReal)
     phase => phases%get(ph)
     thermal => phase%get('thermal',defaultVal=emptyDict)
-    param(ph)%C_p = thermal%get_asFloat('C_p',defaultVal=0.0_pReal)
+    param(ph)%C_p = thermal%get_asFloat('C_p',defaultVal=0.0_pReal)                                 ! ToDo: make mandatory?
     param(ph)%K(1,1) = thermal%get_asFloat('K_11',defaultVal=0.0_pReal)                             ! ToDo: make mandatory?
     param(ph)%K(3,3) = thermal%get_asFloat('K_33',defaultVal=0.0_pReal)                             ! ToDo: depends on symmtery
-    param(ph)%K = lattice_applyLatticeSymmetry33(param(ph)%K,phase%get_asString('lattice'))
+    param(ph)%K = lattice_applyLatticeSymmetry33(param(ph)%K,phase_lattice(ph))
 
     sources => thermal%get('source',defaultVal=emptyList)
     thermal_Nsources(ph) = sources%length