Ignore:
Timestamp:
2018-06-15T17:03:42+02:00 (2 years ago)
Author:
smasson
Message:

trunk: set proper NP folding corners, see #2101

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/LBC/mpp_nfd_generic.h90

    r9690 r9805  
    9595         ALLOCATE( ztabl(jpimax   ,4,ipk,ipl,ipf) , ztabr(jpimax*jpmaxngh,4,ipk,ipl,ipf) )  
    9696         ! 
    97          ztabr(:,:,:,:,:) = 0._wp 
    98          ztabl(:,:,:,:,:) = 0._wp 
    99          ! 
    10097         DO jf = 1, ipf 
    10198            DO jl = 1, ipl 
     
    119116            iproc = nfipproc(isendto(jr),jpnj) 
    120117            IF(iproc /= -1) THEN 
     118               iilb = nimppt(iproc+1) 
     119               ildi = nldit (iproc+1) 
    121120               ilei = nleit (iproc+1) 
    122                ildi = nldit (iproc+1) 
     121               IF( iilb           ==      1 )   ildi = 1     ! e-w boundary already done -> force to take 1st column 
     122               IF( iilb + jpi - 1 == jpiglo )   ilei = jpi   ! e-w boundary already done -> force to take last column 
    123123               iilb = nfiimpp(isendto(jr),jpnj) - nfiimpp(isendto(1),jpnj) 
    124124            ENDIF 
     
    182182            &                znorthgloio, itaille, MPI_DOUBLE_PRECISION, ncomm_north, ierr ) 
    183183         ! 
    184          ztab(:,:,:,:,:) = 0._wp 
    185184         DO jr = 1, ndim_rank_north         ! recover the global north array 
    186185            iproc = nrank_north(jr) + 1 
     186            iilb  = nimppt(iproc) 
    187187            ildi  = nldit (iproc) 
    188188            ilei  = nleit (iproc) 
    189             iilb  = nimppt(iproc) 
     189            IF( iilb           ==      1 )   ildi = 1     ! e-w boundary already done -> force to take 1st column 
     190            IF( iilb + jpi - 1 == jpiglo )   ilei = jpi   ! e-w boundary already done -> force to take last column 
    190191            DO jf = 1, ipf 
    191192               DO jl = 1, ipl 
Note: See TracChangeset for help on using the changeset viewer.