Index: /issm/trunk-jpl/src/c/toolkits/mumps/MpiDenseMumpsSolve.cpp
===================================================================
--- /issm/trunk-jpl/src/c/toolkits/mumps/MpiDenseMumpsSolve.cpp	(revision 16067)
+++ /issm/trunk-jpl/src/c/toolkits/mumps/MpiDenseMumpsSolve.cpp	(revision 16068)
@@ -20,8 +20,14 @@
 /*}}}*/
 
-void MumpsSettings(DMUMPS_STRUC_C &theMumpsStruc,ISSM_MPI_Comm &comm) { 
+void MumpsInit(DMUMPS_STRUC_C &theMumpsStruc) { 
 	theMumpsStruc.par          = 1;  
 	theMumpsStruc.sym          = 0;
-	theMumpsStruc.comm_fortran = MPI_Comm_c2f(comm);
+	theMumpsStruc.comm_fortran = MPI_Comm_c2f(IssmComm::GetComm());
+	theMumpsStruc.job          = -1;
+	dmumps_c(&theMumpsStruc);
+}
+
+// must be preceded by a call to MumpsInit
+void MumpsSettings(DMUMPS_STRUC_C &theMumpsStruc) { 
 	/*Control statements:{{{ */
 	theMumpsStruc.icntl[1-1] = 6; //error verbose
@@ -35,10 +41,4 @@
 	theMumpsStruc.icntl[30-1] = 0;
 	/*}}}*/
-}
-
-// must be preceded by a call to MumpsSettings
-void MumpsInit(DMUMPS_STRUC_C &theMumpsStruc) { 
-	theMumpsStruc.job          = -1;
-	dmumps_c(&theMumpsStruc);
 }
 
@@ -75,8 +75,7 @@
 		IssmPDouble *rhs) { 
 	/*Initialize mumps: {{{*/
-	ISSM_MPI_Comm   comm=IssmComm::GetComm();
 	DMUMPS_STRUC_C theMumpsStruc;
-	MumpsSettings(theMumpsStruc,comm);
 	MumpsInit(theMumpsStruc);
+	MumpsSettings(theMumpsStruc);
 	/*}}}*/
 	// now setup the rest of theMumpsStruc 
@@ -195,4 +194,5 @@
 		   rhs);
 
+if (my_rank==0) for (int i=0;i<Kff_M;++i) std::cout << i << " : " << rhs[i] << std::endl;
 	/*Now scatter from cpu 0 to all other cpus: {{{*/
 	ISSM_MPI_Scatterv( rhs, recvcounts, displs, ISSM_MPI_DOUBLE, uf, uf_m, ISSM_MPI_DOUBLE, 0, comm);