Changeset 27275

Timestamp:

09/09/22 11:48:02 (3 years ago)

Author:

Eric.Larour

Message:

CHG: indices were not being picked up!

File:

• issm/branches/trunk-larour-SLPS2022/src/c/modules/QmuStatisticsx/QmuStatisticsx.cpp (modified) (6 diffs)

issm/branches/trunk-larour-SLPS2022/src/c/modules/QmuStatisticsx/QmuStatisticsx.cpp

@@ -106 +106 @@
         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 +134 @@
                         /*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 +1174 @@
         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 +1217 @@
                 _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 +1256 @@
                 }
                 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 +1293 @@
                         }
                 }
-        }
+        } //}}}
         _printf0_("Done with SampleSeries:\n");
         IssmComm::SetComm(ISSM_MPI_COMM_WORLD);