void read_blottner_data_ascii( MixtureViscosity<BlottnerViscosity<NumericType>,NumericType >& mu,
const std::string &filename )
{
std::ifstream in(filename.c_str());
if(!in.is_open())
{
std::cerr << "ERROR: unable to load file " << filename << std::endl;
antioch_error();
}
// skip the header
skip_comment_lines(in, '#');
std::string name;
NumericType a, b, c;
while (in.good())
{
in >> name; // Species Name
in >> a; //
in >> b; //
in >> c; //
// If we are still good, we have a valid set of transport
// data for this species. Otherwise, we read past end-of-file
// in the section above
if (in.good())
{
const ChemicalMixture<NumericType>& chem_mixture = mu.chemical_mixture();
// Check if this is a species we want.
if( chem_mixture.species_name_map().find(name) !=
chem_mixture.species_name_map().end() )
{
// Pack up coefficients
std::vector<NumericType> coeffs(3);
coeffs[0] = a;
coeffs[1] = b;
coeffs[2] = c;
mu.add(name, coeffs);
}
}
}
in.close();
return;
}
void read_blottner_data_ascii( MixtureViscosity<BlottnerViscosity<NumericType>,NumericType >& mu,
const std::string &filename )
{
std::ifstream in(filename.c_str());
if(!in.is_open())
{
std::cerr << "ERROR: unable to load file " << filename << std::endl;
antioch_error();
}
// skip the header
skip_comment_lines(in, '#');
std::string name;
NumericType a, b, c;
while (in.good())
{
in >> name; // Species Name
in >> a; //
in >> b; //
in >> c; //
// If we are still good, we have a valid set of transport
// data for this species. Otherwise, we read past end-of-file
// in the section above
if (in.good())
{
const ChemicalMixture<NumericType>& chem_mixture = mu.chemical_mixture();
// Check if this is a species we want.
if( chem_mixture.species_name_map().find(name) !=
chem_mixture.species_name_map().end() )
{
// Pack up coefficients
std::vector<NumericType> coeffs(3);
coeffs[0] = a;
coeffs[1] = b;
coeffs[2] = c;
mu.add(name, coeffs);
}
}
}
in.close();
// If we requested Blottner viscosity for our mixture, we'd better
// have Blottner viscosity data for every species in our mixture.
const TransportMixture<NumericType>& trans_mixture = mu.chemical_mixture();
const unsigned int n_species = trans_mixture.n_species();
if (mu.species_viscosities().size() < n_species)
antioch_error_msg
("Could not find Blottner viscosity data for more than " <<
mu.species_viscosities().size() << " of " << n_species <<
" requested species in '" << filename << "'.");
for (unsigned int s=0; s != n_species; ++s)
if (!mu.species_viscosities()[s])
{
const Species& species = trans_mixture.species_list()[s];
const std::string& name = trans_mixture.species_inverse_name_map().find(species)->second;
antioch_error_msg
("Could not find Blottner viscosity data for species '" << name <<
"' in '" << filename << "'.");
}
return;
}
