Often in engineering calculations we want to be able to work with material properties that vary with extensive properties like Temperature and Pressure. For example the specific heat of air varies appreciably with temperature. As in that link, these property variations are determined from experiments and tabulated for specific temperature intervals.
The table represents a numerical function, i.e. it gives
What I'm looking for is a module to code this in python. Is there a scipy or numpy function which takes a pair of arrays, and returns a function which, when called with a Temperature argument, will return the corresponding heat capacity (linearly interpolating where necessary).
To be clear, I don't want to curve fit the data. It's cumbersome and some accuracy is lost.