Changeset 2442 for XIOS2/trunk/src/transformation/domain_algorithm_interpolate.cpp

Timestamp:

12/13/22 11:32:58 (19 months ago)

Author:

jderouillat

Message:

Fix bounds management in transformations. Disabled temporarily bounds settings in the rectilinear case of the generic_testcase

File:

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

XIOS2/trunk/src/transformation/domain_algorithm_interpolate.cpp

@@ -125 +125 @@
   int niSrc = domainSrc_->ni.getValue(), njSrc = domainSrc_->nj.getValue();
   bool hasBoundSrc = domainSrc_->hasBounds;
-  if (hasBoundSrc) nVertexSrc = domainSrc_->nvertex.getValue();
+  if ((hasBoundSrc) && (!domainSrc_->nvertex.isEmpty())) nVertexSrc = domainSrc_->nvertex.getValue(); // default is constNVertex = 4
   CArray<double,2> boundsLonSrc(nVertexSrc,localDomainSrcSize);
   CArray<double,2> boundsLatSrc(nVertexSrc,localDomainSrcSize);
@@ -197 +197 @@
   int niDest = domainDest_->ni.getValue(), njDest = domainDest_->nj.getValue();
   bool hasBoundDest = domainDest_->hasBounds;
-  if (hasBoundDest) nVertexDest = domainDest_->nvertex.getValue();
+  if ((hasBoundDest) && (!domainDest_->nvertex.isEmpty())) nVertexDest = domainDest_->nvertex.getValue(); // default is constNVertex = 4
   CArray<double,2> boundsLonDest(nVertexDest,localDomainDestSize);
   CArray<double,2> boundsLatDest(nVertexDest,localDomainDestSize);