Changeset 753 for XIOS/trunk/src/transformation/domain_algorithm_interpolate.cpp

Timestamp:

10/21/15 16:56:23 (9 years ago)

Author:

mhnguyen

Message:

Updating pole processing for rectangular grid source

Test
+) On Curie
+) test_remap passes

File:

• XIOS/trunk/src/transformation/domain_algorithm_interpolate.cpp (modified) (2 diffs)

XIOS/trunk/src/transformation/domain_algorithm_interpolate.cpp

@@ -76 +76 @@
   else // if domain source is rectilinear, not do anything now
   {
+    bool isNorthPole = false;
+    bool isSouthPole = false;
+    if (poleValue == std::abs(domainSrc_->lat_start)) isNorthPole = true;
+    if (poleValue == std::abs(domainSrc_->lat_end)) isSouthPole = true;
+
     nVertexSrc = constNVertex;
-    domainSrc_->fillInRectilinearBoundLonLat(boundsLonSrc, boundsLatSrc);
-  }
-
-  bool isNorthPole = false;
-  bool isSouthPole = false;
+    domainSrc_->fillInRectilinearBoundLonLat(boundsLonSrc, boundsLatSrc, isNorthPole, isSouthPole);
+  }
+
   std::map<int,std::vector<std::pair<int,double> > > interpMapValueNorthPole;
   std::map<int,std::vector<std::pair<int,double> > > interpMapValueSouthPole;
@@ -116 +119 @@
   else
   {
+    bool isNorthPole = false;
+    bool isSouthPole = false;
     if (poleValue == std::abs(domainDest_->lat_start)) isNorthPole = true;
     if (poleValue == std::abs(domainDest_->lat_end)) isSouthPole = true;