from ChemPlugin import *

print("I made a titration model!\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.01")

for cmd in cmds:
    cp.Config(cmd)

# Initialize the instance
cp.Initialize()
# >>> Write print output to "Titration.txt" <<<
# >>> Write a header to "Titration.gtp" <<<
# >>> Append a data block to "Titration.gtp" <<<

# Time Marching loop
while True:
    deltat = cp.ReportTimeStep()
    if cp.AdvanceTimeStep(deltat): break
    if cp.AdvanceChemical(): break
    # >>> Append a print block <<<
    # >>> Append a plot block <<<

# >>> Do we need a trailer? <<<
input("Done!")