Changeset 2011 for XIOS/dev/dev_ym/XIOS_COUPLING/src/node/interpolate_axis.cpp

Timestamp:

01/12/21 23:05:02 (3 years ago)

Author:

ymipsl

Message:

• bug fix when createing mask on server side when overlapping grid
• implement axis interpolation on pressure coordinate
• big cleaning in transformation

YM

File:

• XIOS/dev/dev_ym/XIOS_COUPLING/src/node/interpolate_axis.cpp (modified) (2 diffs)

XIOS/dev/dev_ym/XIOS_COUPLING/src/node/interpolate_axis.cpp

@@ -1 +1 @@
 #include "interpolate_axis.hpp"
+#include "axis_algorithm_interpolate_coordinate.hpp"
+#include "axis_algorithm_interpolate.hpp"
 #include "type.hpp"
 #include "field.hpp"
@@ -84 +86 @@
     return auxInputs;
   }
+
+  CGenericAlgorithmTransformation* CInterpolateAxis::createAlgorithm(bool isSource,
+                                                        CGrid* gridDst, CGrid* gridSrc,
+                                                        int elementPositionInGrid,
+                                                        std::map<int, int>& elementPositionInGridSrc2ScalarPosition,
+                                                        std::map<int, int>& elementPositionInGridSrc2AxisPosition,
+                                                        std::map<int, int>& elementPositionInGridSrc2DomainPosition,
+                                                        std::map<int, int>& elementPositionInGridDst2ScalarPosition,
+                                                        std::map<int, int>& elementPositionInGridDst2AxisPosition,
+                                                        std::map<int, int>& elementPositionInGridDst2DomainPosition)
+  {
+    if (!coordinate.isEmpty())  return CAxisAlgorithmInterpolateCoordinate::create(isSource, gridDst,  gridSrc, this, elementPositionInGrid, 
+                                      elementPositionInGridSrc2ScalarPosition, elementPositionInGridSrc2AxisPosition, elementPositionInGridSrc2DomainPosition,
+                                      elementPositionInGridDst2ScalarPosition, elementPositionInGridDst2AxisPosition, elementPositionInGridDst2DomainPosition);
+    else return CAxisAlgorithmInterpolate::create(isSource, gridDst,  gridSrc, this, elementPositionInGrid, 
+                elementPositionInGridSrc2ScalarPosition, elementPositionInGridSrc2AxisPosition, elementPositionInGridSrc2DomainPosition,
+                elementPositionInGridDst2ScalarPosition, elementPositionInGridDst2AxisPosition, elementPositionInGridDst2DomainPosition);
+  }
 }