INCLUDE "ChemPlugin.f90"

PROGRAM LinksTree
    Use ChemPluginModule
    IMPLICIT NONE
    TYPE(ChemPlugin), dimension(:), allocatable :: cp
    TYPE(CpiLink) :: link0
    INTEGER :: nlevel = 4,  i, level, error, num_links

    WRITE(*,fmt='("Link ChemPlugin instances into a", I2, " bifurcating tree")') nlevel
    WRITE(*,fmt='("")')
    ALLOCATE(cp(2**nlevel - 1))
    DO i = 1, 2**nlevel - 1
       cp(i) = ChemPlugin()
    END DO

    DO level = 1, nlevel - 1
       Do i = 2**(level-1), 2**level - 1 
          link0 = cp(2*i)%Link(cp(i))
          link0 = cp(2*i + 1)%Link(cp(i))
       END DO
    END DO

   
    DO i = 1, 2**nlevel - 1
       num_links = cp(i)%nLinks()
       WRITE(*,'("Instance", I3, " is linked to", I2, " instance(s)")') i-1, num_Links
    END DO


END PROGRAM LinksTree