Time dependent amplitude in subroutines?

I need to apply non uniform temperature boundary condition to a surface. I will be given an excel file with time and temperature data for n locations on the surface. I was planning on using the disp subroutine to apply this boundary condition. I'm not sure how to incorporate the time/temperature data. I'm having a hard time finding resources for this. So far, for subroutines, I have been able to follow an example using dload to apply pressure to a region in the center of a plate. This is steady state analysis with single magnitude pressure. How can I vary this with time using a table? Can I reference amplitudes defined in the input file or can I reference a table somewhere else? I'm still looking but asking in case I cant find an answer online. Thanks