diff --git a/PRIVATE b/PRIVATE
index 073819b26..226d7b627 160000
--- a/PRIVATE
+++ b/PRIVATE
@@ -1 +1 @@
-Subproject commit 073819b26248e2915bd5f3b6de7d4e0d27a3d423
+Subproject commit 226d7b627e44247b800bce0d9eb7bef1aac6f537
diff --git a/examples/config/phase/damage/isobrittle_generic.yaml b/examples/config/phase/damage/isobrittle_generic.yaml
index 9de78448d..64a51fa42 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]
 
-G: 1400000.0
+G_crit: 1400000.0
 l_c: 1.0
 mu: 0.001
diff --git a/src/homogenization_damage.f90 b/src/homogenization_damage.f90
index 6ae4ac07a..53d6622ac 100644
--- a/src/homogenization_damage.f90
+++ b/src/homogenization_damage.f90
@@ -166,6 +166,7 @@ module subroutine damage_results(ho,group)
 
   integer :: o
 
+
   associate(prm => param(ho))
       outputsLoop: do o = 1,size(prm%output)
         select case(prm%output(o))
diff --git a/src/phase_damage_anisobrittle.f90 b/src/phase_damage_anisobrittle.f90
index f193119eb..167d0829b 100644
--- a/src/phase_damage_anisobrittle.f90
+++ b/src/phase_damage_anisobrittle.f90
@@ -157,7 +157,7 @@ module subroutine anisobrittle_results(phase,group)
     outputsLoop: do o = 1,size(prm%output)
       select case(trim(prm%output(o)))
         case ('f_phi')
-          call results_writeDataset(stt,group,trim(prm%output(o)),'driving force','J/m³')
+          call results_writeDataset(stt,group,trim(prm%output(o)),'driving force','-')
       end select
     end do outputsLoop
   end associate
diff --git a/src/phase_damage_isobrittle.f90 b/src/phase_damage_isobrittle.f90
index 353ce9f00..36f26abd3 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('G')/src%get_asFloat('l_c')
+        prm%W_crit = src%get_asFloat('G_crit')/src%get_asFloat('l_c')
 
 #if defined (__GFORTRAN__)
         prm%output = output_as1dString(src)
@@ -139,7 +139,7 @@ module subroutine isobrittle_results(phase,group)
     outputsLoop: do o = 1,size(prm%output)
       select case(trim(prm%output(o)))
         case ('f_phi')
-          call results_writeDataset(stt,group,trim(prm%output(o)),'driving force','J/m³') ! Wrong, this is dimensionless
+          call results_writeDataset(stt,group,trim(prm%output(o)),'driving force','-')
       end select
     end do outputsLoop