INCLUDE "ChemPlugin.f90"

PROGRAM mReact
    Use ChemPluginModule
    IMPLICIT NONE
    TYPE(ChemPlugin) :: cp
    INTEGER :: file_unit=20, ios=0, error
    REAL(8) :: deltat
    CHARACTER(LEN=255) :: line
    CHARACTER(LEN=255) :: input_file
    
    WRITE(*,'("mReact -- Use ChemPlugin to emulate React")')
    WRITE(*,'("")')

    ! Create ChemPlugin instance and capture output messsages
    cp = ChemPlugin("stdout")
    error = cp%Config("delxi = 0.01; step_increase = 1.5; pluses = banner")
    error = cp%Config("print = on; plot = on")  
    
    WRITE(*,'( "Enter React input script: ")', ADVANCE='no')
    READ(*,'(A)') input_file
    
    OPEN(unit=file_unit, file=input_file, action="read")
    Do WHILE(ios == 0) 
       READ(file_unit,'(A)',IOSTAT=ios) line
       IF ( ios == 0 ) THEN
          IF (line /= 'go') THEN
             error = cp%Config(line)
          ELSE
             error = cp%Initialize()
             DO WHILE(.TRUE.) 
                deltat = cp%ReportTimeStep()
                IF (cp%AdvanceTimeStep(deltat) /= 0) THEN
                   EXIT
                END IF
                IF (cp%AdvanceHeatTransport() /= 0) THEN
                   EXIT
                END IF
                IF (cp%AdvanceChemical() /= 0) THEN
                   EXIT
                END IF
             END DO
          END IF
       END IF
    END DO
     

END PROGRAM mReact