added thermo_solve.py