Changeset 676 for XIOS/trunk/src/node/field.cpp

Timestamp:

08/25/15 16:52:45 (9 years ago)

Author:

rlacroix

Message:

Add support for indexed output.

If the new field attribute "indexed_output" is set to true and a mask is defined (either at grid, domain or axis level), the indexed data will be outputed instead of the full data with missing values.

See ​http://cfconventions.org/Data/cf-conventions/cf-conventions-1.5/build/cf-conventions.html#compression-by-gathering for more information.

File:

• XIOS/trunk/src/node/field.cpp (modified) (6 diffs)

XIOS/trunk/src/node/field.cpp

@@ -35 +35 @@
       , hasOutputFile(false)
       , domAxisIds_("", ""), areAllReferenceSolved(false)
+      , useCompressedOutput(false)
       , isReadDataRequestPending(false)
    { setVirtualVariableGroup(); }
@@ -45 +46 @@
       , hasOutputFile(false)
       , domAxisIds_("", ""), areAllReferenceSolved(false)
+      , useCompressedOutput(false)
       , isReadDataRequestPending(false)
    { setVirtualVariableGroup(); }
@@ -455 +457 @@
    {
       return !this->refObject.empty();
+   }
+
+   //----------------------------------------------------------------
+
+   bool CField::getUseCompressedOutput() const
+   {
+     return useCompressedOutput;
+   }
+
+   void CField::setUseCompressedOutput()
+   {
+     useCompressedOutput = true;
    }
 
@@ -814 +828 @@
       for(it=data_srv.begin();it!=data_srv.end();it++)
       {
-         grid->outputField(it->first, it->second, fieldOut.dataFirst()) ;
+         grid->outputField(it->first, it->second, fieldOut.dataFirst());
       }
    }
@@ -824 +838 @@
       for (it = data_srv.begin(); it != data_srv.end(); it++)
       {
-         grid->outputField(it->first, it->second, fieldOut.dataFirst()) ;
+         grid->outputField(it->first, it->second, fieldOut.dataFirst());
       }
    }
@@ -852 +866 @@
       {
          grid->inputField(it->first, fieldOut.dataFirst(), it->second);
+      }
+   }
+
+   void CField::outputCompressedField(CArray<double,1>& fieldOut)
+   {
+      map<int, CArray<double,1> >::iterator it;
+
+      for (it = data_srv.begin(); it != data_srv.end(); it++)
+      {
+         grid->outputCompressedField(it->first, it->second, fieldOut.dataFirst());
       }
    }