Changeset 1571 for XIOS/dev/dev_olga/src/node/axis.hpp

Timestamp:

09/06/18 18:55:44 (6 years ago)

Author:

oabramkina

Message:

Taking care of cases when a process doesn't possess a grid (its domain or axe size is zero).

Tested with CMIP6 toy models and IPSL model.

To do: remove grid mask from transformations.

File:

• XIOS/dev/dev_olga/src/node/axis.hpp (modified) (3 diffs)

XIOS/dev/dev_olga/src/node/axis.hpp

@@ -114 +114 @@
          void computeWrittenIndex();
          void computeWrittenCompressedIndex(MPI_Comm);
+
+         template <class T> void outputAxis(CArray<T,1>& valueIn, CArray<T,1>& valueOut);
+
          bool hasTransformation();
          void solveInheritanceTransformation();
@@ -126 +129 @@
         bool hasLabel;
 
-        CArray<size_t,1> localIndexToWriteOnServer;
+//        CArray<size_t,1> localIndexToWriteOnServer;
+        CArray<int,1> localIndexToWriteOnServer;
+        CArray<size_t,1> localIndexToWriteOnClient;
 
       private:
@@ -185 +190 @@
    ///--------------------------------------------------------------
 
+   template <class T> void CAxis::outputAxis(CArray<T,1>& valueIn, CArray<T,1>& valueOut)
+   {
+     int sizeServer = localIndexToWriteOnServer.numElements();
+     int sizeClient = localIndexToWriteOnClient.numElements();
+     valueOut.resize(sizeClient);
+
+     if (sizeServer == sizeClient)
+     {
+       valueOut = valueIn;
+     }
+     else
+     // a hole
+     {
+       valueOut = valueIn;
+       //valueOut = -1; initialize here to some value ????
+     }
+//      CArray<size_t,1>& outIndexClient = grid->localIndexToWriteOnClient;
+//      CArray<size_t,1>& outIndexServer = grid->localIndexToWriteOnServer;
+//      for (size_t idx = 0; idx < outIndexServer.numElements(); ++idx)
+//      {
+//        fieldOut(outIndexServer(idx)) = recvDataSrv(outIndexClient(idx));
+//      }
+   }
+
    // Declare/Define CAxisGroup and CAxisDefinition
    DECLARE_GROUP(CAxis);