diff --git a/code/plastic_isotropic.f90 b/code/plastic_isotropic.f90
index 81d3055cf..d8dc0f384 100644
--- a/code/plastic_isotropic.f90
+++ b/code/plastic_isotropic.f90
@@ -45,7 +45,7 @@ module plastic_isotropic
      gdot0, &
      n, &
      h0, &
-     h0_slopeLnRate, &
+     h0_slopeLnRate=0.0_pReal, &
      tausat, &
      a, &
      aTolFlowstress, &
@@ -466,7 +466,7 @@ subroutine plastic_isotropic_LiAndItsTangent(Li,dLi_dTstar_3333,Tstar_v,ipc,ip,e
    math_spherical33, &
    math_mul33xx33
  use material, only: &
-   phaseAt, phasememberAt, &
+   phasememberAt, &
    plasticState, &
    material_phase, &
    phase_plasticityInstance
@@ -495,7 +495,7 @@ subroutine plastic_isotropic_LiAndItsTangent(Li,dLi_dTstar_3333,Tstar_v,ipc,ip,e
    k, l, m, n
 
  of = phasememberAt(ipc,ip,el)                                                                      ! phasememberAt should be tackled by material and be renamed to material_phasemember
- instance = phase_plasticityInstance(phaseAt(ipc,ip,el))                                            ! "phaseAt" equivalent to "material_phase" !!
+ instance = phase_plasticityInstance(material_phase(ipc,ip,el))
 
  Tstar_sph_33 = math_spherical33(math_Mandel6to33(Tstar_v))                                         ! spherical part of 2nd Piola-Kirchhoff stress
  squarenorm_Tstar_sph = math_mul33xx33(Tstar_sph_33,Tstar_sph_33)
@@ -535,7 +535,7 @@ subroutine plastic_isotropic_dotState(Tstar_v,ipc,ip,el)
  use math, only: &
    math_mul6x6
  use material, only: &
-   phaseAt, phasememberAt, &
+   phasememberAt, &
    plasticState, &
    material_phase, &
    phase_plasticityInstance
@@ -559,7 +559,7 @@ subroutine plastic_isotropic_dotState(Tstar_v,ipc,ip,el)
    of                                                                                               !< shortcut notation for offset position in state array
 
  of = phasememberAt(ipc,ip,el)                                                                      ! phasememberAt should be tackled by material and be renamed to material_phasemember
- instance = phase_plasticityInstance(phaseAt(ipc,ip,el))                                            ! "phaseAt" equivalent to "material_phase" !!
+ instance = phase_plasticityInstance(material_phase(ipc,ip,el))
 
 !--------------------------------------------------------------------------------------------------
 ! norm of (deviatoric) 2nd Piola-Kirchhoff stress
@@ -616,7 +616,7 @@ function plastic_isotropic_postResults(Tstar_v,ipc,ip,el)
  use material, only: &
    material_phase, &
    plasticState, &
-   phaseAt, phasememberAt, &
+   phasememberAt, &
    phase_plasticityInstance
 
  implicit none
@@ -640,7 +640,7 @@ function plastic_isotropic_postResults(Tstar_v,ipc,ip,el)
    o
 
  of = phasememberAt(ipc,ip,el)                                                                      ! phasememberAt should be tackled by material and be renamed to material_phasemember
- instance = phase_plasticityInstance(phaseAt(ipc,ip,el))                                            ! "phaseAt" equivalent to "material_phase" !!
+ instance = phase_plasticityInstance(material_phase(ipc,ip,el))
  
 !--------------------------------------------------------------------------------------------------
 ! norm of (deviatoric) 2nd Piola-Kirchhoff stress