Clustered solutions around harmonic centers to a coupled elliptic system