diff --git a/code/DAMASK_spectral_driver.f90 b/code/DAMASK_spectral_driver.f90
index 85c856076..2e1fd51cf 100644
--- a/code/DAMASK_spectral_driver.f90
+++ b/code/DAMASK_spectral_driver.f90
@@ -536,6 +536,7 @@ program DAMASK_spectral_Driver
          if (.not. cutBack) &
            write(statUnit,*) totalIncsCounter, time, cutBackLevel, &
                              solres%converged, solres%iterationsNeeded                              ! write statistics about accepted solution
+           flush(statUnit)
        enddo subIncLooping
        cutBackLevel = max(0_pInt, cutBackLevel - 1_pInt)                                            ! try half number of subincs next inc
        if(solres%converged) then                                                                    ! report converged inc
@@ -550,6 +551,7 @@ program DAMASK_spectral_Driver
        if (mod(inc,loadCases(currentLoadCase)%outputFrequency) == 0_pInt) then                      ! at output frequency
          write(6,'(1/,a)') ' ... writing results to file ......................................'
          write(resUnit)  materialpoint_results                                                      ! write result to file
+         flush(resUnit)
        endif
        if( loadCases(currentLoadCase)%restartFrequency > 0_pInt .and. &                             ! at frequency of writing restart information set restart parameter for FEsolving 
                       mod(inc,loadCases(currentLoadCase)%restartFrequency) == 0_pInt) then          ! ToDo first call to CPFEM_general will write? 
diff --git a/code/constitutive.f90 b/code/constitutive.f90
index 6d25802dc..34366e0be 100644
--- a/code/constitutive.f90
+++ b/code/constitutive.f90
@@ -154,6 +154,8 @@ subroutine constitutive_init
  character(len=32) :: outputName                                                                    !< name of output, intermediate fix until HDF5 output is ready
  logical :: knownPlasticity, nonlocalConstitutionPresent
 #ifdef HDF
+ integer(pInt) :: &
+  phase                                                                                             !< phase
  integer(pInt), dimension(:,:,:), allocatable :: mappingConstitutive
  integer(pInt), dimension(:,:,:), allocatable :: mappingCrystallite
  integer(pInt), dimension(:), allocatable :: ConstitutivePosition
@@ -265,8 +267,8 @@ subroutine constitutive_init
        end select
        instance = phase_plasticityInstance(material_phase(g,i,e))
 #ifdef HDF
-       InstancePosition(instance) = InstancePosition(instance)+1_pInt
-       mapping(g,e,1:2) = [instancePosition(instance),instance]
+       ConstitutivePosition(instance) = ConstitutivePosition(instance)+1_pInt
+       mappingConstitutive(g,e,1:2)   = [ConstitutivePosition(instance),instance]
 #endif
        select case(phase_plasticity(material_phase(g,i,e)))
          case (PLASTICITY_NONE_ID)
@@ -446,7 +448,7 @@ subroutine constitutive_init
    endforall
  enddo
 #ifdef HDF
- call  HDF5_mappingConstitutive(mapping)
+ call  HDF5_mappingConstitutive(mappingConstitutive)
 #endif
 !--------------------------------------------------------------------------------------------------
 ! write out state size file