# - computeFluidVsp (PyTree) -
import Converter.PyTree as C
import Converter.Internal as Internal
import Thermolib.PyTree as Thermolib

Internal.__FlowSolutionCenters__ = 'SolutionFlow'

t = C.convertFile2PyTree("Simple/archive_CHARME.hdf")

h = Thermolib.Thermo('epicea.xml', workdir='Simple')

#h._computeFluidVsp(t, ['centers:H', 'centers:Cp'])
h._computeFluidVsp(t, ['H', 'Cp'], 'SolutionFlow')

C.convertPyTree2File(t, 'out.hdf')