diff --git a/code/constitutive_dislotwin.f90 b/code/constitutive_dislotwin.f90
index 94ed789eb..8bbf2adba 100644
--- a/code/constitutive_dislotwin.f90
+++ b/code/constitutive_dislotwin.f90
@@ -424,7 +424,7 @@ subroutine constitutive_dislotwin_init(fileUnit)
 ! parameters depending on slip number of twin families
        case ('ntwin')
          if (positions(1) < Nchunks_TwinFamilies + 1_pInt) &
-           call IO_warning(50_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_DISLOTWIN_label//')')
+           call IO_warning(51_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_DISLOTWIN_label//')')
          if (positions(1) > Nchunks_TwinFamilies + 1_pInt) &
            call IO_error(150_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_DISLOTWIN_label//')')
          Nchunks_TwinFamilies = positions(1) - 1_pInt
diff --git a/code/constitutive_phenopowerlaw.f90 b/code/constitutive_phenopowerlaw.f90
index 08cce95b6..efd211c7c 100644
--- a/code/constitutive_phenopowerlaw.f90
+++ b/code/constitutive_phenopowerlaw.f90
@@ -317,7 +317,7 @@ allocate(constitutive_phenopowerlaw_sizePostResults(maxNinstance),
 ! parameters depending on number of twin families
        case ('ntwin')
          if (positions(1) < Nchunks_TwinFamilies + 1_pInt) &
-           call IO_warning(50_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_PHENOPOWERLAW_label//')')
+           call IO_warning(51_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_PHENOPOWERLAW_label//')')
          if (positions(1) > Nchunks_TwinFamilies + 1_pInt) &
            call IO_error(150_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_PHENOPOWERLAW_label//')')
          Nchunks_TwinFamilies = positions(1) - 1_pInt
diff --git a/code/crystallite.f90 b/code/crystallite.f90
index 633c7b5b6..ffdb2b3a3 100644
--- a/code/crystallite.f90
+++ b/code/crystallite.f90
@@ -3389,9 +3389,9 @@ logical function crystallite_stateJump(g,i,e)
      crystallite_stateJump = .false.
      return
    endif
-   plasticState(mappingConstitutive(2,g,i,e))%state(1:mySizeDotState,mappingConstitutive(1,g,i,e)) =                &
-                plasticState(mappingConstitutive(2,g,i,e))%state(1:mySizeDotState,mappingConstitutive(1,g,i,e))     &
-              + plasticState(mappingConstitutive(2,g,i,e))%deltaState(1:mySizeDotState,mappingConstitutive(1,g,i,e))
+   plasticState(  mappingConstitutive(2,g,i,e))%state     (1:mySizeDotState,mappingConstitutive(1,g,i,e)) =  &
+     plasticState(mappingConstitutive(2,g,i,e))%state     (1:mySizeDotState,mappingConstitutive(1,g,i,e)) +  &
+     plasticState(mappingConstitutive(2,g,i,e))%deltaState(1:mySizeDotState,mappingConstitutive(1,g,i,e))
 #else 
    mySizeDotState = constitutive_sizeDotState(g,i,e)
    if (any(constitutive_deltaState(g,i,e)%p(1:mySizeDotState) &