Index: /issm/branches/trunk-larour-SLPS2022/src/c/modules/QmuStatisticsx/QmuStatisticsx.cpp
===================================================================
--- /issm/branches/trunk-larour-SLPS2022/src/c/modules/QmuStatisticsx/QmuStatisticsx.cpp	(revision 27274)
+++ /issm/branches/trunk-larour-SLPS2022/src/c/modules/QmuStatisticsx/QmuStatisticsx.cpp	(revision 27275)
@@ -106,4 +106,5 @@
 	iomodel->FindConstant(&numstatistics,"md.qmu.statistics.numstatistics");
 	for (int i=1;i<=numstatistics;i++){
+	//for (int i=9;i<=9;i++){
 
 		char* directory=NULL;
@@ -133,5 +134,5 @@
 			/*fetch indices: */
 			sprintf(string,"md.qmu.statistics.method(%i).indices",i);
-			iomodel->FetchData(&indices,&dummy,&nindices,string);
+			iomodel->FetchData(&indices,&nindices,&dummy,string);
 			parameters->AddObject(new IntVecParam(IndicesEnum,indices,nindices));
 
@@ -1173,7 +1174,22 @@
 	xsize=xNew<int>(nfields*nsteps);
 
+	/*_printf0_("nindices: " << nindices << "\n");
+	for (int i=0;i<nindices;i++){
+		_printf0_(indices[i] << " " );
+	}
+	_printf0_("\n");
+
+	_printf0_("nsteps: " << nsteps << "\n");
+	for (int i=0;i<nsteps;i++){
+		_printf0_(steps[i] << " " );
+	}
+	_printf0_("\n");*/
+
+
+
 	/*Start opening files:*/
 	for (int i=(lower_row+1);i<=upper_row;i++){
 		_printf0_("reading file #: " << i << "\n");
+		/*Create memory buffer for file, to speed things up: {{{*/
 		char file[1000];
 		long int  length;
@@ -1201,6 +1217,6 @@
 		_printf0_("    processing file:\n");
 		fid=fmemopen(buffer, length, "rb");
-
-		/*start reading data from the buffer directly:*/
+		/*}}}*/
+		/*Retrieve the values for all fields and time steps:{{{*/
 		for (int f=0;f<nfields;f++){
 			fseek(fid,0,SEEK_SET);
@@ -1240,10 +1256,9 @@
 		}
 		fclose(fid);
-
-		/*delete buffer:*/
 		xDelete<char>(buffer);
+		/*}}}*/
 	}
 	ISSM_MPI_Barrier(IssmComm::GetComm());
-	_printf0_("Done reading files, now assembling time series.\n");
+	_printf0_("Done reading files, now assembling time series.\n"); //{{{
 
 	for (int f=0;f<nfields;f++){
@@ -1278,5 +1293,5 @@
 			}
 		}
-	}
+	} //}}}
 	_printf0_("Done with SampleSeries:\n");
 	IssmComm::SetComm(ISSM_MPI_COMM_WORLD);