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