Defining a history function for Delay Differential Equations using past data points

You would just use an interpolation like DataInterpolations.jl for the h(p,t) that you give to the solver. That will then by used for any t < tspan[1].