INCLUDE "ChemPlugin.f90"

PROGRAM LinksChain
    Use ChemPluginModule
    IMPLICIT NONE
    TYPE(ChemPlugin) :: cp_inlet
    TYPE(ChemPlugin), dimension(:), allocatable :: cp
    TYPE(CpiLink) :: link0
    INTEGER :: nchain = 8,  i, error, num_links

    WRITE(*,fmt=('("Link ChemPlugin instances into a one-dimensional chain")'))
    WRITE(*,fmt=('("")'))

    nChain = 8
    cp_inlet = ChemPlugin()
    ALLOCATE(cp(nChain))
    DO i = 1, nChain
       cp(i) = ChemPlugin()
    END DO

    ! Link the instances into a chain
    link0 = cp(1)%Link(cp_inlet)
    Do i = 2, nChain
       link0 = cp(i)%Link(cp(i-1))
    END DO 
    link0 = cp(nChain)%Link() ! Outlet
    num_links = cp_inlet%nLinks()
    WRITE(*,'("Inlet is linked to ", I1 , " instance")') num_links
    Do i = 1, nChain
        num_links = cp(i)%nLinks()
       WRITE(*,'("Instance", I2, " is linked to ", I1, " instances")') i - 1, num_Links
    END DO

END PROGRAM LinksChain