from ChemPlugin import *

print("ChemPlugin example -- pH titration\n")

# Create ChemPlugin instance
cp = ChemPlugin("stdout")

# Configure the instance
cmds =  ("Ca++ = 1 mmol/kg",
         "Na+ = 1 mmol/kg",
         "Cl- = 3 mmol/kg",
         "HCO3- = 2 mmol/kg",
         "pH = 4")
        
cmds += ("react 3 mmol/kg NaOH", "delxi = 0.1")

for cmd in cmds:
    cp.Config(cmd)

# Initialize the instance
cp.Initialize()
cp.PlotHeader("myPlot.gtp", "char")
cp.PlotBlock()

# Time Marching loop
while True:
    deltat = cp.ReportTimeStep()
    if cp.AdvanceTimeStep(deltat): break
    if cp.AdvanceChemical(): break
    cp.PlotBlock()
    cp.PlotTrailer()
    print("Pausing at Xi = %4.2f" % cp.Report1("Xi"))
    input()

input()