INCLUDE "ChemPlugin.f90"

PROGRAM simple
    USE ChemPluginModule
    IMPLICIT NONE
    REAL(8) :: deltat
    REAL(8)::  pH(1)
    TYPE(ChemPlugin) :: cp
    INTEGER :: error, i
     CHARACTER(LEN=255), dimension(2), PARAMETER :: &
        cmds = [character(len=255) :: &
        "Na+ = 1 mmol/kg; Cl- = 1 mmol/kg; pH = 3","react 2 mmol NaOH; delxi = 0.1"]

    cp = ChemPlugin()
   
    Do i = 1, 2
       error = cp%Config(cmds(i))
    end do
    
    error = cp%Initialize()
    error = cp%Report(pH,"pH")
     WRITE(*,'("pH = ",F8.4)') pH
     do while(.TRUE.)
        deltat = cp%ReportTimeStep() 
        if (cp%AdvanceTimeStep(deltat) /= 0) then
           EXIT
        end if
        if (cp%AdvanceChemical() /= 0) then
            EXIT
        end if 
        error = cp%Report(pH, "pH")
        WRITE(*,'("pH = ",F8.4)') pH
    end do
END PROGRAM simple