Good answers above, but there is an additional point you might need to note if you are crossing a phase boundary. Some parameters, eg density, can change very rapidly across a phase boundary and a simple interpolation may not always give a good result. This problem can be minimised by using a very detailed matrix of values to be interpolated, or by using a more complicated interpolation algorithm than just fitting splines and using the interpolation function.
Alan