26 #include "grins_config.h"
28 #ifdef GRINS_HAVE_ANTIOCH
34 #include "libmesh/getpot.h"
37 #include "antioch/read_reaction_set_data_xml.h"
38 #include "antioch/cea_mixture_ascii_parsing.h"
39 #include "antioch/stat_mech_thermo.h"
45 _reaction_set( new Antioch::ReactionSet<
libMesh::Real>( (*_antioch_gas.get()) ) ),
46 _cea_mixture( new Antioch::CEAThermoMixture<
libMesh::Real>( (*_antioch_gas.get()) ) )
48 if( !input.have_variable(
"Physics/Chemistry/chem_file") )
50 std::cerr <<
"Error: Must specify XML chemistry file to use Antioch." << std::endl;
54 std::string xml_filename = input(
"Physics/Chemistry/chem_file",
"DIE!");
55 bool verbose_read = input(
"screen-options/verbose_kinetics_read",
false );
57 Antioch::read_reaction_set_data_xml<libMesh::Real>( xml_filename, verbose_read, *
_reaction_set.get() );
59 Antioch::read_cea_mixture_data_ascii_default( *
_cea_mixture.get() );
73 Antioch::StatMechThermodynamics<libMesh::Real> thermo( *(this->
_antioch_gas.get()) );
77 for(
unsigned int s = 0; s < this->
n_species(); s++ )
87 #endif // GRINS_HAVE_ANTIOCH
unsigned int n_species() const
boost::scoped_ptr< Antioch::CEAThermoMixture< libMesh::Real > > _cea_mixture
Wrapper class for Antioch::ChemicalMixture.
boost::scoped_ptr< Antioch::ReactionSet< libMesh::Real > > _reaction_set
virtual ~AntiochMixture()
void build_stat_mech_ref_correction()
std::vector< libMesh::Real > _h_stat_mech_ref_correction
boost::scoped_ptr< Antioch::ChemicalMixture< libMesh::Real > > _antioch_gas