diff --git a/src/plastic_kinematichardening.f90 b/src/plastic_kinematichardening.f90
index a82d9066a..3bd35aaea 100644
--- a/src/plastic_kinematichardening.f90
+++ b/src/plastic_kinematichardening.f90
@@ -1,7 +1,7 @@
 !--------------------------------------------------------------------------------------------------
 !> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH
 !> @author Zhuowen Zhao, Michigan State University
-!> @brief Introducing Voce-type kinematic hardening rule into crystal phenopowerlaw plasticity  
+!> @brief Introducing Voce-type kinematic hardening rule into crystal plasticity  
 !! formulation using a power law fitting
 !--------------------------------------------------------------------------------------------------
 module plastic_kinehardening
@@ -318,7 +318,7 @@ subroutine plastic_kinehardening_init(fileUnit)
          do j = 1_pInt, Nchunks_SlipSlip
            param(instance)%interaction_slipslip(j) = IO_floatValue(line,chunkPos,1_pInt+j)
          enddo
-       case ('nonSchmidCoeff')
+       case ('nonschmidcoeff')
          if (chunkPos(1) < 1_pInt + Nchunks_nonSchmid) &
            call IO_warning(52_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_KINEHARDENING_label//')')
          do j = 1_pInt,Nchunks_nonSchmid
@@ -350,10 +350,10 @@ subroutine plastic_kinehardening_init(fileUnit)
        case ('n_slip')
          param(instance)%n_slip                   = IO_floatValue(line,chunkPos,2_pInt)
              
-       case ('aTolResistance')
+       case ('atol_resistance')
          param(instance)%aTolResistance           = IO_floatValue(line,chunkPos,2_pInt)
        
-       case ('aTolShear')
+       case ('atol_shear')
          param(instance)%aTolShear                = IO_floatValue(line,chunkPos,2_pInt)
                
        case default
diff --git a/src/plastic_phenopowerlaw.f90 b/src/plastic_phenopowerlaw.f90
index 319f1b585..facc76ff1 100644
--- a/src/plastic_phenopowerlaw.f90
+++ b/src/plastic_phenopowerlaw.f90
@@ -987,13 +987,15 @@ subroutine plastic_phenopowerlaw_dotState(Tstar_v,ipc,ip,el)
    index_myFamily = sum(lattice_NslipSystem(1:f-1_pInt,ph))                                         ! at which index starts my family
    slipSystems1: do i = 1_pInt,plastic_phenopowerlaw_Nslip(f,instance)
      j = j+1_pInt
-     left_SlipSlip(j) = 1.0_pReal + plastic_phenopowerlaw_H_int(f,instance)                         ! modified no system-dependent left part
-     left_SlipTwin(j) = 1.0_pReal                                                                   ! no system-dependent left part
-     right_SlipSlip(j) = abs(1.0_pReal-plasticState(ph)%state(j,of) / &
-                                    (plastic_phenopowerlaw_tausat_slip(f,instance)+ssat_offset)) &
-                         **plastic_phenopowerlaw_a_slip(instance)&
-                         *sign(1.0_pReal,1.0_pReal-plasticState(ph)%state(j,of) / &
-                                    (plastic_phenopowerlaw_tausat_slip(f,instance)+ssat_offset))
+     left_SlipSlip(j) = (1.0_pReal + plastic_phenopowerlaw_H_int(f,instance)) &
+                        *abs(1.0_pReal-plasticState(ph)%state(j,of) / &                             ! no system-dependent left part
+                                      (plastic_phenopowerlaw_tausat_slip(f,instance)+ssat_offset)) &
+                        **plastic_phenopowerlaw_a_slip(instance)&
+                        *sign(1.0_pReal,1.0_pReal-plasticState(ph)%state(j,of) / &
+                                      (plastic_phenopowerlaw_tausat_slip(f,instance)+ssat_offset))              
+     left_SlipTwin(j) = 1.0_pReal                                                                   
+     right_SlipSlip(j) = 1.0_pReal                                                                  ! system-dependent part (beta summation)
+                 
      right_TwinSlip(j) = 1.0_pReal                                                                  ! no system-dependent part
 
 !--------------------------------------------------------------------------------------------------