# - computeFlowVars (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._computeFlowVars(t, ['centers:Mach', 'centers:Etot'])
h._computeFlowVars(t, ['Mach', 'Etot'], 'SolutionFlow')

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