diff --git a/src/source_damage_anisoBrittle.f90 b/src/source_damage_anisoBrittle.f90
index c380e9790..4a9ae1f68 100644
--- a/src/source_damage_anisoBrittle.f90
+++ b/src/source_damage_anisoBrittle.f90
@@ -91,6 +91,8 @@ subroutine source_damage_anisoBrittle_init(fileUnit)
    IO_error, &
    IO_timeStamp, &
    IO_EOF
+    use math, only: &
+   math_expand
  use material, only: &
    material_allocateSourceState, &
    phase_source, &
@@ -175,7 +177,14 @@ subroutine source_damage_anisoBrittle_init(fileUnit)
    if (prm%sdot_0   <= 0.0_pReal) extmsg = trim(extmsg)//' anisobrittle_sdot0'
    
    prm%Ncleavage = config%getInts('ncleavage',defaultVal=emptyIntArray)
-   
+
+   prm%critDisp = config%getFloats('anisobrittle_criticaldisplacement',requiredSize=size(prm%Ncleavage))
+   prm%critLoad = config%getFloats('anisobrittle_criticalload',        requiredSize=size(prm%Ncleavage))
+
+     ! expand: family => system
+     prm%critDisp  = math_expand(prm%critDisp, prm%Ncleavage)
+     prm%critLoad  = math_expand(prm%critLoad, prm%Ncleavage)
+     
 !--------------------------------------------------------------------------------------------------
 !  exit if any parameter is out of range
    if (extmsg /= '') &
diff --git a/src/source_damage_anisoDuctile.f90 b/src/source_damage_anisoDuctile.f90
index 46898ecf5..94e2b3a4a 100644
--- a/src/source_damage_anisoDuctile.f90
+++ b/src/source_damage_anisoDuctile.f90
@@ -93,6 +93,8 @@ subroutine source_damage_anisoDuctile_init(fileUnit)
    IO_error, &
    IO_timeStamp, &
    IO_EOF
+    use math, only: &
+   math_expand
  use material, only: &
    material_allocateSourceState, &
    phase_source, &
@@ -176,6 +178,13 @@ subroutine source_damage_anisoDuctile_init(fileUnit)
    
    prm%Nslip  = config%getInts('nslip',defaultVal=emptyIntArray)
    
+   prm%critPlasticStrain = config%getFloats('anisoductile_criticalplasticstrain',requiredSize=size(prm%Nslip))
+   prm%critLoad          = config%getFloats('anisoductile_criticalload',         requiredSize=size(prm%Nslip))
+
+     ! expand: family => system
+     prm%critPlasticStrain   = math_expand(prm%critPlasticStrain,  prm%Nslip)
+     prm%critLoad            = math_expand(prm%critLoad,           prm%Nslip)
+   
 !--------------------------------------------------------------------------------------------------
 !  exit if any parameter is out of range
    if (extmsg /= '') &