# - computeFluidVars (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._computeFluidVars(t, ['centers:Rho', 'centers:H'])
h._computeFluidVars(t, ['Rho', 'H'], 'SolutionFlow')

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