Changeset 9850


Ignore:
Timestamp:
2018-06-27T17:12:25+02:00 (2 years ago)
Author:
smasson
Message:

dev_r9759_HPC09_ESIWACE: report bugfix from NEMOtrunk on not-masked NP, see #2101

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r9759_HPC09_ESIWACE/src/OCE/LBC/mpp_nfd_generic.h90

    r9814 r9850  
    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 
     
    123120            iproc = nfipproc(isendto(jr),jpnj) 
    124121            IF(iproc /= -1) THEN 
     122               iilb = nimppt(iproc+1) 
     123               ildi = nldit (iproc+1) 
    125124               ilei = nleit (iproc+1) 
    126                ildi = nldit (iproc+1) 
     125               IF( iilb           ==      1 )   ildi = 1     ! e-w boundary already done -> force to take 1st column 
     126               IF( iilb + jpi - 1 == jpiglo )   ilei = jpi   ! e-w boundary already done -> force to take last column 
    127127               iilb = nfiimpp(isendto(jr),jpnj) - nfiimpp(isendto(1),jpnj) 
    128128            ENDIF 
     
    199199         ENDIF 
    200200         ! 
    201          ztab(:,:,:,:,:) = 0._wp 
    202201         DO jr = 1, ndim_rank_north         ! recover the global north array 
    203202            iproc = nrank_north(jr) + 1 
     203            iilb  = nimppt(iproc) 
    204204            ildi  = nldit (iproc) 
    205205            ilei  = nleit (iproc) 
    206             iilb  = nimppt(iproc) 
     206            IF( iilb           ==      1 )   ildi = 1     ! e-w boundary already done -> force to take 1st column 
     207            IF( iilb + jpi - 1 == jpiglo )   ilei = jpi   ! e-w boundary already done -> force to take last column 
    207208            DO jf = 1, ipf 
    208209               DO jl = 1, ipl 
Note: See TracChangeset for help on using the changeset viewer.