Changing young's modulus

After each step, I want to change the young's modulus of a material. I have written a subroutine in fortran to calculate the new elastic modulus for each time step but I am not sure how or which subroutine to use to transfer the new elastic modulus to Abaqus for each time step.