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

Timestamp:

11/16/15 16:21:45 (9 years ago)

Author:

rlacroix

Message:

Add a new field attribute grid_path to allow chaining spatial transformations.

Intermediate grids must be separated by comma and the source/destination grids must not be listed.

File:

• XIOS/trunk/src/node/field.cpp (modified) (1 diff)

XIOS/trunk/src/node/field.cpp

@@ -819 +819 @@
    {
      if (grid && !grid->isTransformed() && hasDirectFieldReference() && grid != getDirectFieldReference()->grid)
-       grid->transformGrid(getDirectFieldReference()->grid);
+     {
+       std::vector<CGrid*> grids;
+       // Source grid
+       grids.push_back(getDirectFieldReference()->grid);
+       // Intermediate grids
+       if (!grid_path.isEmpty())
+       {
+         std::string gridId;
+         size_t start = 0, end;
+
+         do
+         {
+           end = grid_path.getValue().find(',', start);
+           if (end != std::string::npos)
+           {
+             gridId = grid_path.getValue().substr(start, end - start);
+             start = end + 1;
+           }
+           else
+             gridId = grid_path.getValue().substr(start);
+
+           if (!CGrid::has(gridId))
+             ERROR("void CField::solveTransformedGrid()",
+                   << "Invalid grid_path, the grid '" << gridId << "' does not exist.");
+
+           grids.push_back(CGrid::get(gridId));
+         }
+         while (end != std::string::npos);
+       }
+       // Destination grid
+       grids.push_back(grid);
+
+       for (size_t i = 0, count = grids.size() - 1; i < count; ++i)
+       {
+         CGrid *gridSrc  = grids[i];
+         CGrid *gridDest = grids[i + 1];
+         if (!gridDest->isTransformed())
+           gridDest->transformGrid(gridSrc);
+       }
+     }
    }