#include <iostream>
#include <iomanip>
#include "ChemPlugin.h"


int main(int argc, char** argv) {
   std::cout << "ChemPlugin example -- pH titration" << std::endl << std::endl;

   // Create a ChemPlugin instance.
   ChemPlugin cp("stdout");

   // Configure the instance.
   cp.Config("Ca++ = 1 mmol/kg; Na+ = 1 mmol/kg");
   cp.Config("Cl- = 3 mmol/kg; HCO3- = 2 mmol/kg; pH = 4");
   cp.Config("react 3 mmol/kg NaOH; delxi = 0.1");

   // Initialize the instance.
   cp.Initialize();
   cp.PrintOutput("Initial.txt");

   // Time marching loop.
   while (true) {
      double deltat = cp.ReportTimeStep();
      if (cp.AdvanceTimeStep(deltat)) break;
      if (cp.AdvanceChemical()) break;
      cp.PrintOutput("Path.txt");
   }

   // Any keystroke closes the console.
   std::cin.get();
   return 0;
}