use warnings;
use strict;

use ChemPlugin;

$|=1;

print "ChemPlugin example -- pH titration\n\n";

# create the ChemPlugin instance
my $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();

while(1){
   my $deltat = $cp->ReportTimeStep();
   if($cp->AdvanceTimeStep($deltat)){
      last; 
   }
   if($cp->AdvanceChemical()){
      last;
   }
   $cp->PlotBlock();
}

$cp->destroy();

