diff --git a/PRIVATE b/PRIVATE
index 459326e98..dd7f4efd6 160000
--- a/PRIVATE
+++ b/PRIVATE
@@ -1 +1 @@
-Subproject commit 459326e9840c843ade72b04cf28e50889c9779f1
+Subproject commit dd7f4efd64ebadac74cab84f713a738c71d8b737
diff --git a/examples/config/phase/damage/anisobrittle_cubic.yaml b/examples/config/phase/damage/anisobrittle_cubic.yaml
index 372fdc6d0..4f8db50ee 100644
--- a/examples/config/phase/damage/anisobrittle_cubic.yaml
+++ b/examples/config/phase/damage/anisobrittle_cubic.yaml
@@ -9,5 +9,5 @@ s_crit: [0.006666]
 dot_o: 1.e-3
 q: 20
 
-D_11: 1.0
+l_c: 1.0
 mu: 0.001
diff --git a/examples/config/phase/damage/isobrittle_generic.yaml b/examples/config/phase/damage/isobrittle_generic.yaml
index 851302a49..9de78448d 100644
--- a/examples/config/phase/damage/isobrittle_generic.yaml
+++ b/examples/config/phase/damage/isobrittle_generic.yaml
@@ -2,6 +2,6 @@ type: isobrittle
 
 output: [f_phi]
 
-W_crit: 1400000.0
-D_11: 1.0
+G: 1400000.0
+l_c: 1.0
 mu: 0.001
diff --git a/src/phase_damage.f90 b/src/phase_damage.f90
index 9022131e2..5bf351dd0 100644
--- a/src/phase_damage.f90
+++ b/src/phase_damage.f90
@@ -105,9 +105,7 @@ module subroutine damage_init
       damage_active = .true.
       source => sources%get(1)
       param(ph)%mu     = source%get_asFloat('mu')
-      param(ph)%D(1,1) = source%get_asFloat('D_11')
-      if (any(phase_lattice(ph) == ['hP','tI'])) param(ph)%D(3,3) = source%get_asFloat('D_33')
-      param(ph)%D = lattice_symmetrize_33(param(ph)%D,phase_lattice(ph))
+      param(ph)%D  = math_I3 * source%get_asFloat('l_c')**2
     end if
 
   end do
diff --git a/src/phase_damage_isobrittle.f90 b/src/phase_damage_isobrittle.f90
index 2d501b58f..353ce9f00 100644
--- a/src/phase_damage_isobrittle.f90
+++ b/src/phase_damage_isobrittle.f90
@@ -63,7 +63,7 @@ module function isobrittle_init() result(mySources)
       associate(prm => param(ph), dlt => deltaState(ph), stt => state(ph))
         src => sources%get(1)
 
-        prm%W_crit = src%get_asFloat('W_crit')
+        prm%W_crit = src%get_asFloat('G')/src%get_asFloat('l_c')
 
 #if defined (__GFORTRAN__)
         prm%output = output_as1dString(src)