Changeset 2119

Timestamp:

03/30/21 11:06:23 (3 years ago)

Author:

jderouillat

Message:

Fix in reorder_domain

File:

• XIOS/dev/dev_ym/XIOS_COUPLING/src/transformation/domain_algorithm/domain_algorithm_reorder.cpp (modified) (3 diffs)

XIOS/dev/dev_ym/XIOS_COUPLING/src/transformation/domain_algorithm/domain_algorithm_reorder.cpp

@@ -44 +44 @@
 TRY
 {
+  domainDestination->type.setValue( CDomain::type_attr::rectilinear );
+  domainDestination->ni_glo = domainSource->ni_glo;
+  domainDestination->nj_glo = domainSource->nj_glo;
+  domainDestination->ni = domainSource->ni;
+  domainDestination->nj = domainSource->nj;
+  domainDestination->ibegin = domainSource->ibegin;
+  domainDestination->jbegin = domainSource->jbegin;
+
   reorderDomain->checkValid(domainSource);
   domainDestination->checkAttributes() ; // for now but maybe use domainSource as template for domain destination
 
-  if (domainDestination->type !=  CDomain::type_attr::rectilinear)
+  if (domainSource->type !=  CDomain::type_attr::rectilinear)
   {
       ERROR("CDomainAlgorithmReorder::CDomainAlgorithmReorder(CDomain* domainDestination, CDomain* domainSource, CReorderDomain* reorderDomain)",
@@ -62 +70 @@
   }
   
-  if (!reorderDomain->invert_lat.isEmpty())
+  if (!reorderDomain->invert_lat.isEmpty() && reorderDomain->invert_lat.getValue() )
   {
     CArray<int,1>& j_index=domainDestination->j_index ;
     int nglo = j_index.numElements() ;
     int nj_glo =domainDestination->nj_glo ;
-    
     for (size_t i = 0; i < nglo ; ++i)
     {
@@ -80 +87 @@
     CArray<int,1>& i_index=domainDestination->i_index ;
     int nglo = i_index.numElements() ;
-    
     for (size_t i = 0; i < nglo ; ++i)
     {