Return thermodynamic properties at given temperature, density and composition, by EOS
pure function thermo_dTX input Density d; input Temperature T; input Real z[subsNum]; output Real p, h, s, Cv, Cp, Dvp, DvT; end thermo_dTX;