import ChemPlugin.ChemPlugin;

public class Titration5 {

   public static void main(String[] args) {

      System.out.println("ChemPlugin example -- pH titration\n");

      // create ChemPlugin instance
      ChemPlugin cp = new ChemPlugin("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.PlotHeader("myPlot.gtp", "char");
      cp.PlotBlock();

      // Time marching loop
      while (true) {
         double deltat = cp.ReportTimeStep();
         if (cp.AdvanceTimeStep(deltat) != 0){
            System.exit(1);
         }
         if (cp.AdvanceChemical() != 0){
            System.exit(0);
         }
         cp.PlotBlock();
      }

   }

}