Changeset 5487

Timestamp:

11/02/20 09:55:31 (4 years ago)

Author:

ymipsl

Message:

XIOS modification to implement working ensemble by pool of members. Only slice of pool memeber are written in the netcdf file.

YM

File:

• CONFIG_DEVT/LMDZOR_V6.2_work_ENSEMBLES/modeles/XIOS/src/io/nc4_data_output.cpp (modified) (1 diff)

CONFIG_DEVT/LMDZOR_V6.2_work_ENSEMBLES/modeles/XIOS/src/io/nc4_data_output.cpp

@@ -2634 +2634 @@
                 }
 
-
+                // suppose first dim is dedicated for ensemble
+                double defaultValue=field->default_value ;
+                size_t slice=1 ;
+                for(int i=1 ; i<count.size() ; i++) slice=slice*count[i] ;
+                
+                size_t nEnsemble=count[0] ;
+                size_t pos=0 ;
+                start[0] = 0 ;
+                count[0] = 0 ;
+                for(size_t i=0 ; i<nEnsemble ; i++)
+                {
+                  bool allMissing=true ;
+                  for (size_t j=0 ; j<slice ; j++,pos++)
+                    if (fieldData(pos)!=defaultValue) { allMissing=false ; break;}
+                  if (allMissing) 
+                  {
+                    if (count[0]==0) start[0]++ ;
+                    else break ;
+                  }
+                  else count[0]++ ;  
+                }  
+
+                CArray<double,1> newData(fieldData.dataFirst()+start[0]*slice,shape(count[0]*slice),neverDeleteData) ;
+                info(0)<<"write data parameter  "<<fieldid<<"  "<<field->getNStep()<<"  start : " ;
+                for(auto it : start) info(0)<<" "<<it ;
+                info(0)<<"  count : " ; 
+                for(auto it : count) info(0)<<"  "<<it ;
+                info(0)<<endl ; 
+                 
                 CTimer::get("Files : writing data").resume();
-                SuperClassWriter::writeData(fieldData, fieldid, isCollective, field->getNStep() - 1, &start, &count);
+                SuperClassWriter::writeData(newData, fieldid, isCollective, field->getNStep() - 1, &start, &count);
                 CTimer::get("Files : writing data").suspend();