diff --git a/src/source_damage_anisoBrittle.f90 b/src/source_damage_anisoBrittle.f90
index a8f9de6f6..d9ec6f34c 100644
--- a/src/source_damage_anisoBrittle.f90
+++ b/src/source_damage_anisoBrittle.f90
@@ -56,6 +56,8 @@ module source_damage_anisoBrittle
      totalNcleavage
    integer(pInt), dimension(:), allocatable :: &
      Ncleavage
+   integer(kind(undefined_ID)), allocatable, dimension(:) :: &
+     outputID                                                                                       !< ID of each post result output
  end type tParameters
 
  type(tParameters), dimension(:), allocatable, private :: param                                     !< containers of constitutive parameters (len Ninstance)
@@ -112,8 +114,6 @@ subroutine source_damage_anisoBrittle_init(fileUnit)
    config_phase, &
    material_Nphase, &
    MATERIAL_partPhase
- use numerics,only: &
-   numerics_integrator
  use lattice, only: &
    lattice_maxNcleavageFamily, &
    lattice_NcleavageSystem
@@ -123,17 +123,22 @@ subroutine source_damage_anisoBrittle_init(fileUnit)
 
  integer(pInt), allocatable, dimension(:) :: chunkPos
  integer(pInt) :: Ninstance,mySize=0_pInt,phase,instance,source,sourceOffset,o
- integer(pInt) :: sizeState, sizeDotState, sizeDeltaState
- integer(pInt) :: NofMyPhase,p   
+ integer(pInt) :: NofMyPhase,p   ,i
  integer(pInt) :: Nchunks_CleavageFamilies = 0_pInt, j
- character(len=pStringLen) :: &
-   extmsg = ''
  character(len=65536) :: &
    tag  = '', &
    line = ''
  integer(pInt),          dimension(0), parameter :: emptyIntArray    = [integer(pInt)::]
+ character(len=65536),   dimension(0), parameter :: emptyStringArray = [character(len=65536)::]
+ integer(kind(undefined_ID)) :: &
+   outputID
 
- write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_damage_anisoBrittle_LABEL//' init  -+>>>'
+ character(len=pStringLen) :: &
+   extmsg = ''
+ character(len=65536), dimension(:), allocatable :: &
+   outputs
+
+ write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_DAMAGE_ANISOBRITTLE_LABEL//' init  -+>>>'
  write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"
 
@@ -185,6 +190,24 @@ subroutine source_damage_anisoBrittle_init(fileUnit)
    
    prm%Ncleavage = config%getInts('ncleavage',defaultVal=emptyIntArray)
    
+!--------------------------------------------------------------------------------------------------
+!  exit if any parameter is out of range
+   if (extmsg /= '') &
+     call IO_error(211_pInt,ext_msg=trim(extmsg)//'('//SOURCE_DAMAGE_ANISOBRITTLE_LABEL//')')
+
+!--------------------------------------------------------------------------------------------------
+!  output pararameters
+   outputs = config%getStrings('(output)',defaultVal=emptyStringArray)
+   allocate(prm%outputID(0))
+   do i=1_pInt, size(outputs)
+     outputID = undefined_ID
+     select case(outputs(i))
+       case ('anisobrittle_drivingforce')
+
+     end select
+
+   enddo
+
    end associate
    
  enddo
diff --git a/src/source_damage_anisoDuctile.f90 b/src/source_damage_anisoDuctile.f90
index 7abc751d6..925588594 100644
--- a/src/source_damage_anisoDuctile.f90
+++ b/src/source_damage_anisoDuctile.f90
@@ -57,6 +57,8 @@ module source_damage_anisoDuctile
      totalNslip
    integer(pInt), dimension(:), allocatable :: &
      Nslip
+   integer(kind(undefined_ID)), allocatable, dimension(:) :: &
+     outputID
  end type tParameters
 
  type(tParameters), dimension(:), allocatable, private :: param                                     !< containers of constitutive parameters (len Ninstance)
@@ -113,8 +115,6 @@ subroutine source_damage_anisoDuctile_init(fileUnit)
    config_phase, &
    material_Nphase, &
    MATERIAL_partPhase
- use numerics,only: &
-   numerics_integrator
  use lattice, only: &
    lattice_maxNslipFamily, &
    lattice_NslipSystem
@@ -124,17 +124,22 @@ subroutine source_damage_anisoDuctile_init(fileUnit)
 
  integer(pInt), allocatable, dimension(:) :: chunkPos
  integer(pInt) :: Ninstance,mySize=0_pInt,phase,instance,source,sourceOffset,o
- integer(pInt) :: sizeState, sizeDotState, sizeDeltaState
- integer(pInt) :: NofMyPhase,p   
- integer(pInt) :: Nchunks_SlipFamilies = 0_pInt, j   
- character(len=pStringLen) :: &
-   extmsg = ''
+ integer(pInt) :: NofMyPhase,p ,i
+ integer(pInt) :: Nchunks_SlipFamilies = 0_pInt, j
  character(len=65536) :: &
    tag  = '', &
    line = ''
  integer(pInt),          dimension(0), parameter :: emptyIntArray    = [integer(pInt)::]
+ character(len=65536),   dimension(0), parameter :: emptyStringArray = [character(len=65536)::]
+ integer(kind(undefined_ID)) :: &
+   outputID
 
- write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_damage_anisoDuctile_LABEL//' init  -+>>>'
+ character(len=pStringLen) :: &
+   extmsg = ''
+ character(len=65536), dimension(:), allocatable :: &
+   outputs
+
+ write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_DAMAGE_ANISODUCTILE_LABEL//' init  -+>>>'
  write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"
 
@@ -184,6 +189,24 @@ subroutine source_damage_anisoDuctile_init(fileUnit)
    
    prm%Nslip  = config%getInts('nslip',defaultVal=emptyIntArray)
    
+!--------------------------------------------------------------------------------------------------
+!  exit if any parameter is out of range
+   if (extmsg /= '') &
+     call IO_error(211_pInt,ext_msg=trim(extmsg)//'('//SOURCE_DAMAGE_ANISODUCTILE_LABEL//')')
+
+!--------------------------------------------------------------------------------------------------
+!  output pararameters
+   outputs = config%getStrings('(output)',defaultVal=emptyStringArray)
+   allocate(prm%outputID(0))
+   do i=1_pInt, size(outputs)
+     outputID = undefined_ID
+     select case(outputs(i))
+       case ('anisoductile_drivingforce')
+
+     end select
+
+   enddo
+
    end associate
    
  enddo
diff --git a/src/source_damage_isoBrittle.f90 b/src/source_damage_isoBrittle.f90
index 3a2481639..e09d79056 100644
--- a/src/source_damage_isoBrittle.f90
+++ b/src/source_damage_isoBrittle.f90
@@ -39,6 +39,8 @@ module source_damage_isoBrittle
      critStrainEnergy, &
      N, &
      aTol
+   integer(kind(undefined_ID)), allocatable, dimension(:) :: &
+     outputID
  end type tParameters
 
  type(tParameters), dimension(:), allocatable, private :: param                                     !< containers of constitutive parameters (len Ninstance)
@@ -95,23 +97,26 @@ subroutine source_damage_isoBrittle_init(fileUnit)
    config_phase, &
    material_Nphase, &
    MATERIAL_partPhase
- use numerics,only: &
-   numerics_integrator
 
  implicit none
  integer(pInt), intent(in) :: fileUnit
 
  integer(pInt), allocatable, dimension(:) :: chunkPos
  integer(pInt) :: Ninstance,mySize=0_pInt,phase,instance,source,sourceOffset,o
- integer(pInt) :: sizeState, sizeDotState, sizeDeltaState
- integer(pInt) :: NofMyPhase,p   
- character(len=pStringLen) :: &
-   extmsg = ''
+ integer(pInt) :: NofMyPhase,p,i   
  character(len=65536) :: &
    tag  = '', &
    line = ''
+ character(len=65536),   dimension(0), parameter :: emptyStringArray = [character(len=65536)::]
+ integer(kind(undefined_ID)) :: &
+   outputID
 
- write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_damage_isoBrittle_label//' init  -+>>>'
+ character(len=pStringLen) :: &
+   extmsg = ''
+ character(len=65536), dimension(:), allocatable :: &
+   outputs
+ 
+ write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_DAMAGE_ISOBRITTLE_LABEL//' init  -+>>>'
  write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"
 
@@ -156,6 +161,24 @@ subroutine source_damage_isoBrittle_init(fileUnit)
    if (prm%N                  <= 0.0_pReal) extmsg = trim(extmsg)//' isobrittle_n'
    if (prm%critStrainEnergy   <= 0.0_pReal) extmsg = trim(extmsg)//' isobrittle_criticalstrainenergy'
    
+!--------------------------------------------------------------------------------------------------
+!  exit if any parameter is out of range
+   if (extmsg /= '') &
+     call IO_error(211_pInt,ext_msg=trim(extmsg)//'('//SOURCE_DAMAGE_ISOBRITTLE_LABEL//')')
+
+!--------------------------------------------------------------------------------------------------
+!  output pararameters
+   outputs = config%getStrings('(output)',defaultVal=emptyStringArray)
+   allocate(prm%outputID(0))
+   do i=1_pInt, size(outputs)
+     outputID = undefined_ID
+     select case(outputs(i))
+       case ('isobrittle_drivingforce')
+       
+     end select
+
+   enddo
+
    end associate
    
  enddo
diff --git a/src/source_damage_isoDuctile.f90 b/src/source_damage_isoDuctile.f90
index 7186f8749..3b4b06727 100644
--- a/src/source_damage_isoDuctile.f90
+++ b/src/source_damage_isoDuctile.f90
@@ -40,6 +40,8 @@ module source_damage_isoDuctile
      critPlasticStrain, &
      N, &
      aTol
+   integer(kind(undefined_ID)), allocatable, dimension(:) :: &
+     outputID
  end type tParameters
 
  type(tParameters), dimension(:), allocatable, private :: param                                     !< containers of constitutive parameters (len Ninstance)
@@ -96,24 +98,26 @@ subroutine source_damage_isoDuctile_init(fileUnit)
    config_phase, &
    material_Nphase, &
    MATERIAL_partPhase
- 
- use numerics,only: &
-   numerics_integrator
 
  implicit none
  integer(pInt), intent(in) :: fileUnit
 
  integer(pInt), allocatable, dimension(:) :: chunkPos
  integer(pInt) :: Ninstance,mySize=0_pInt,phase,instance,source,sourceOffset,o
- integer(pInt) :: sizeState, sizeDotState, sizeDeltaState
- integer(pInt) :: NofMyPhase,p   
- character(len=pStringLen) :: &
-   extmsg = ''
+ integer(pInt) :: NofMyPhase,p,i
  character(len=65536) :: &
    tag  = '', &
    line = ''
+ character(len=65536),   dimension(0), parameter :: emptyStringArray = [character(len=65536)::]
+ integer(kind(undefined_ID)) :: &
+   outputID
 
- write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_damage_isoDuctile_label//' init  -+>>>'
+ character(len=pStringLen) :: &
+   extmsg = ''
+ character(len=65536), dimension(:), allocatable :: &
+   outputs
+
+ write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_DAMAGE_ISODUCTILE_LABEL//' init  -+>>>'
  write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"
 
@@ -157,6 +161,24 @@ subroutine source_damage_isoDuctile_init(fileUnit)
    
    if (prm%N                   <= 0.0_pReal) extmsg = trim(extmsg)//' isoductile_ratesensitivity'
    if (prm%critPlasticStrain   <= 0.0_pReal) extmsg = trim(extmsg)//' isoductile_criticalplasticstrain'
+   
+!--------------------------------------------------------------------------------------------------
+!  exit if any parameter is out of range
+   if (extmsg /= '') &
+     call IO_error(211_pInt,ext_msg=trim(extmsg)//'('//SOURCE_DAMAGE_ISODUCTILE_LABEL//')')
+
+!--------------------------------------------------------------------------------------------------
+!  output pararameters
+   outputs = config%getStrings('(output)',defaultVal=emptyStringArray)
+   allocate(prm%outputID(0))
+   do i=1_pInt, size(outputs)
+     outputID = undefined_ID
+     select case(outputs(i))
+       case ('isoductile_drivingforce')
+       
+     end select
+
+   enddo
 
    end associate