Is there a way to incorporate time series data into coupled ODE system using DifferentialEquations.jl?

Have a look at How to linearize a model I created the function cp from numerical data and used it in Modelingtoolkit equations…

Yes, you might need another equation…