count_trial Subroutine

    subroutine count_trial()
        use Parallel_neci, only: MPISumAll
        implicit none
        integer ::  ntrialtot, ncontot, ntrial, ncon

        call count_trial_this_proc(ntrial, ncon)
        call MPISumAll(ntrial, ntrialtot)
        call MPISumAll(ncon, ncontot)
        write(stdout, *) "Trial states ", ntrialtot
        write(stdout, *) "Connected states ", ncontot
    end subroutine count_trial