diff --git a/trunk/constitutive.f90 b/trunk/constitutive.f90 index 16e32b1c6..fbd862003 100644 --- a/trunk/constitutive.f90 +++ b/trunk/constitutive.f90 @@ -933,7 +933,7 @@ enddo dLp_dTstar=0.0_pReal do i=1,material_Nslip(matID) dgdot_dtauslip(i) = material_gdot0_slip(matID)*(abs(tau_slip(i))/state(i))**& - (material_n_slip(matID)-1.0_pReal)*material_n_slip(matID)/constitutive_state_new(i,ipc,ip,el) + (material_n_slip(matID)-1.0_pReal)*material_n_slip(matID)/state(i) forall (k=1:3,l=1:3,m=1:3,n=1:3) dLp_dTstar(k,l,m,n) = dLp_dTstar(k,l,m,n)+ & dgdot_dtauslip(i)*constitutive_Sslip(k,l,i,material_CrystalStructure(matID))* &