Index: /issm/branches/trunk-larour-SLPS2020/src/c/modules/QmuStatisticsx/QmuStatisticsx.cpp
===================================================================
--- /issm/branches/trunk-larour-SLPS2020/src/c/modules/QmuStatisticsx/QmuStatisticsx.cpp	(revision 25598)
+++ /issm/branches/trunk-larour-SLPS2020/src/c/modules/QmuStatisticsx/QmuStatisticsx.cpp	(revision 25599)
@@ -669,4 +669,5 @@
 	int nfields;
 	int range,lower_row,upper_row;
+	int nfilesperdirectory;
 	
 	/*intermediary:*/
@@ -688,5 +689,9 @@
 	int*         meansize=NULL;
 
+	/*only work on the statistical communicator: */
+	if (color==MPI_UNDEFINED)return 0;
+
 	/*Retrieve parameters:*/
+	parameters->FindParam(&nfilesperdirectory,QmuNfilesPerDirectoryEnum);
 	parameters->FindParam(&nsamples,QmuNsampleEnum);
 	parameters->FindParam(&directory,DirectoryNameEnum);
@@ -695,5 +700,6 @@
 	parameters->FindParam(&steps,&nsteps,StepsEnum);
 
-	/*Get rank:*/
+	/*Get rank from the stat comm communicator:*/
+	IssmComm::SetComm(statcomm);
 	int my_rank=IssmComm::GetRank();
 
@@ -983,4 +989,6 @@
 	}
 
+	_printf0_("Done with MeanVariance :\n");
+	IssmComm::SetComm(ISSM_MPI_COMM_WORLD);
 
 } /*}}}*/
@@ -1139,6 +1147,6 @@
 
 } /*}}}*/
-int OutputStatistics(Parameters* parameters,Results* results){ /*{{{*/
-	
+int OutputStatistics(Parameters* parameters,Results* results,int color,ISSM_MPI_Comm statcomm){ /*{{{*/
+
 	char   outputfilename[1000];
 	char* directory=NULL;
@@ -1149,4 +1157,7 @@
 	int nsteps;
 
+	/*only work on the statistical communicator: */
+	if (color==MPI_UNDEFINED)return 0;
+
 	FemModel* femmodel=new FemModel();
 	
Index: /issm/branches/trunk-larour-SLPS2020/src/c/modules/QmuStatisticsx/QmuStatisticsx.h
===================================================================
--- /issm/branches/trunk-larour-SLPS2020/src/c/modules/QmuStatisticsx/QmuStatisticsx.h	(revision 25598)
+++ /issm/branches/trunk-larour-SLPS2020/src/c/modules/QmuStatisticsx/QmuStatisticsx.h	(revision 25599)
@@ -9,5 +9,5 @@
 int ComputeMeanVariance(Parameters* parameters,Results* results,int color, ISSM_MPI_Comm statcomm);
 int ComputeSampleSeries(Parameters* parameters,Results* results,int color, ISSM_MPI_Comm statcomm);
-int OutputStatistics(Parameters* parameters,Results* results);
+int OutputStatistics(Parameters* parameters,Results* results,int color,ISSM_MPI_Comm statcomm);
 int ComputeHistogram(Parameters* parameters,Results* results,int color, ISSM_MPI_Comm statcomm);
 int readdata(IssmDouble** pdoublemat, int* pdoublematsize, IssmDouble* pdouble, FILE* fid,char* field,int step);