The goal is to obtain shrinkage curves from the macroscopic approach, I am using CREEP subroutine. The base code from the picture is not having any variable to give relative density as an input. In my case, the swelling strain is the function of relative density. Please give me any suggestions regarding this problem.