Ignore:
Timestamp:
01/29/14 11:49:29 (10 years ago)
Author:
ymipsl
Message:

Add functionnalities for reduced gaussian grid support.

YM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • XIOS/trunk/src/node/grid.cpp

    r435 r467  
    425425    { 
    426426       rank=domain->connectedServer[ns] ; 
    427        int ib=domain->ib_srv[ns] ; 
    428        int ie=domain->ie_srv[ns] ; 
    429        int jb=domain->jb_srv[ns] ; 
    430        int je=domain->je_srv[ns] ; 
    431427        
    432428       int i,j ; 
     
    434430       for(int k=0;k<storeIndex_client.numElements();k++) 
    435431       { 
    436          i=out_i_client(k) ; 
    437          j=out_j_client(k) ; 
    438          if (i>=ib-1 && i<=ie-1 && j>=jb-1 && j<=je-1) nb++ ;  
     432         i=out_i_client(k)- domain->ibegin +1; 
     433         j=out_j_client(k)- domain->jbegin +1; 
     434         if (domain->mapConnectedServer(i,j)==ns)  nb++ ; 
    439435       } 
    440436        
     
    448444       for(int k=0;k<storeIndex_client.numElements();k++) 
    449445       { 
    450          i=out_i_client(k) ; 
    451          j=out_j_client(k) ; 
    452          if (i>=ib-1 && i<=ie-1 && j>=jb-1 && j<=je-1)  
     446         i=out_i_client(k)- domain->ibegin +1 ; 
     447         j=out_j_client(k)- domain->jbegin +1 ; 
     448         if (domain->mapConnectedServer(i,j)==ns) 
    453449         { 
    454450            storeIndex(nb)=k ; 
    455             out_i(nb)=out_i_client(k) ; 
    456             out_j(nb)=out_j_client(k) ; 
     451            out_i(nb)=domain->i_index(i,j) + domain->ibegin-1; 
     452            out_j(nb)=domain->j_index(i,j) + domain->jbegin-1;  
    457453            out_l(nb)=out_l_client(k) ; 
    458454            nb++ ; 
Note: See TracChangeset for help on using the changeset viewer.