Changeset 2119


Ignore:
Timestamp:
03/30/21 11:06:23 (18 months ago)
Author:
jderouillat
Message:

Fix in reorder_domain

File:
1 edited

Legend:

Unmodified
Added
Removed
  • XIOS/dev/dev_ym/XIOS_COUPLING/src/transformation/domain_algorithm/domain_algorithm_reorder.cpp

    r2011 r2119  
    4444TRY 
    4545{ 
     46  domainDestination->type.setValue( CDomain::type_attr::rectilinear ); 
     47  domainDestination->ni_glo = domainSource->ni_glo; 
     48  domainDestination->nj_glo = domainSource->nj_glo; 
     49  domainDestination->ni = domainSource->ni; 
     50  domainDestination->nj = domainSource->nj; 
     51  domainDestination->ibegin = domainSource->ibegin; 
     52  domainDestination->jbegin = domainSource->jbegin; 
     53 
    4654  reorderDomain->checkValid(domainSource); 
    4755  domainDestination->checkAttributes() ; // for now but maybe use domainSource as template for domain destination 
    4856 
    49   if (domainDestination->type !=  CDomain::type_attr::rectilinear) 
     57  if (domainSource->type !=  CDomain::type_attr::rectilinear) 
    5058  { 
    5159      ERROR("CDomainAlgorithmReorder::CDomainAlgorithmReorder(CDomain* domainDestination, CDomain* domainSource, CReorderDomain* reorderDomain)", 
     
    6270  } 
    6371   
    64   if (!reorderDomain->invert_lat.isEmpty()) 
     72  if (!reorderDomain->invert_lat.isEmpty() && reorderDomain->invert_lat.getValue() ) 
    6573  { 
    6674    CArray<int,1>& j_index=domainDestination->j_index ; 
    6775    int nglo = j_index.numElements() ; 
    6876    int nj_glo =domainDestination->nj_glo ; 
    69      
    7077    for (size_t i = 0; i < nglo ; ++i) 
    7178    { 
     
    8087    CArray<int,1>& i_index=domainDestination->i_index ; 
    8188    int nglo = i_index.numElements() ; 
    82      
    8389    for (size_t i = 0; i < nglo ; ++i) 
    8490    { 
Note: See TracChangeset for help on using the changeset viewer.