-c Corrected AL 8/19/2014: each processor needs whole iset array not only its
-c own element
-c call mpi_scatter(i2set,1,mpi_integer,
-c & iset,1,mpi_integer,king,
-c & CG_COMM,ierr)
- call mpi_bcast(i2set(0),nodes,mpi_integer,king,
- & CG_COMM,ierr)
- iset=i2set(me)
+ call mpi_scatter(i2set,1,mpi_integer,
+ & iset,1,mpi_integer,king,
+ & CG_COMM,ierr)
+