Changeset 1581 for XIOS


Ignore:
Timestamp:
09/19/18 14:30:22 (6 years ago)
Author:
ymipsl
Message:

XIOS remapper : bug fix for 2nd order conservative remapping with masked cells leading to crash in some cases.

YM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • XIOS/trunk/extern/remap/src/meshutil.cpp

    r877 r1581  
    133133  { 
    134134    for (int i = 0; i < elts[j]->n; i++)  
    135       if ( elts[j]->neighbour[i]== NOT_FOUND) neighbours[i]=NULL ; 
     135      if ( elts[j]->neighbour[i]== NOT_FOUND) neighbours[i]=NULL ; // no neighbour 
     136      else if (elts[elts[j]->neighbour[i]]->is.size() == 0) neighbours[i]=NULL ; // neighbour has none supermesh cell  
    136137      else  neighbours[i] = elts[elts[j]->neighbour[i]]; 
    137138 
Note: See TracChangeset for help on using the changeset viewer.