26{
29
30 std::string line;
31 std::string aqueous_solution_prev_;
32 std::size_t chemical_system_id = 0;
33 while (std::getline(in, line))
34 {
35 if (line.find("USE reaction_pressure none") != std::string::npos)
36 {
37 break;
38 }
39
40 if (line.find("SOLUTION_RAW") != std::string::npos)
41 {
42 aqueous_solution_prev_ =
43 "SOLUTION_RAW " +
44 std::to_string(num_chemical_systems + chemical_system_id + 1) +
45 "\n";
46 continue;
47 }
48
49 aqueous_solution_prev_ += line + "\n";
50
51 if (line.find("-gammas") != std::string::npos)
52 {
54 aqueous_solution_prev_.clear();
55 ++chemical_system_id;
56 }
57 }
58}
std::vector< std::string > aqueous_solutions_prev