Changeset 10537


Ignore:
Timestamp:
2019-01-16T21:41:21+01:00 (15 months ago)
Author:
smasson
Message:

trunk: bugfix in bdyini, see #2213

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/BDY/bdyini.F90

    r10481 r10537  
    941941                     ii = idx_bdy(ib_bdy)%nbi(icount,igrd) 
    942942                     ij = idx_bdy(ib_bdy)%nbj(icount,igrd) 
    943                      if((com_east .ne. 1) .and. (ii == (nlci-1)) .and. (nbondi .le. 0)) then 
     943                     if((com_east .ne. 1) .and. (ii == nlci) .and. (nbondi .le. 0)) then 
    944944                        com_east = 1 
    945945                     elseif((com_west .ne. 1) .and. (ii == 2) .and. (nbondi .ge. 0) .and. (nbondi .ne. 2)) then 
     
    948948                     if((com_south .ne. 1) .and. (ij == 2) .and. (nbondj .ge. 0) .and. (nbondj .ne. 2)) then 
    949949                        com_south = 1 
    950                      elseif((com_north .ne. 1) .and. (ij == (nlcj-1)) .and. (nbondj .le. 0)) then 
     950                     elseif((com_north .ne. 1) .and. (ij == nlcj) .and. (nbondj .le. 0)) then 
    951951                        com_north = 1 
    952952                     endif  
     
    960960                       & nbrdta(ib,igrd,ib_bdy) == ir  ) THEN 
    961961                       ii = nbidta(ib,igrd,ib_bdy)- iw_b(1)+2 
    962                        if((com_west_b .ne. 1) .and. (ii == (nlcit(nowe+1)-1))) then 
     962                       if((com_west_b .ne. 1) .and. (ii == (nlcit(nowe+1)))) then 
    963963                          ij = nbjdta(ib,igrd,ib_bdy) - is_b(1)+2 
    964964                          if((ij == 2) .and. (nbondj == 0 .or. nbondj == 1)) then 
     
    989989                       & nbrdta(ib,igrd,ib_bdy) == ir  ) THEN 
    990990                       ii = nbidta(ib,igrd,ib_bdy)- iw_b(1)+2 
    991                        if((com_west_b .ne. 1) .and. (ii == (nlcit(nowe+1)-1))) then 
     991                       if((com_west_b .ne. 1) .and. (ii == (nlcit(nowe+1)))) then 
    992992                          ij = nbjdta(ib,igrd,ib_bdy) - is_b(1)+2 
    993993                          if((ij == 2) .and. (nbondj == 0 .or. nbondj == 1)) then 
     
    10301030                       & nbrdta(ib,igrd,ib_bdy) == ir  ) THEN 
    10311031                       ij = nbjdta(ib,igrd,ib_bdy)- is_b(3)+2 
    1032                        if((com_south_b .ne. 1) .and. (ij == (nlcjt(noso+1)-1))) then 
     1032                       if((com_south_b .ne. 1) .and. (ij == (nlcjt(noso+1)))) then 
    10331033                          com_south_b = 1 
    10341034                       endif  
     
    10521052                       & nbrdta(ib,igrd,ib_bdy) == ir  ) THEN 
    10531053                       ij = nbjdta(ib,igrd,ib_bdy)- is_b(3)+2 
    1054                        if((com_south_b .ne. 1) .and. (ij == (nlcjt(noso+1)-1))) then 
     1054                       if((com_south_b .ne. 1) .and. (ij == (nlcjt(noso+1)))) then 
    10551055                          com_south_b = 1 
    10561056                       endif  
Note: See TracChangeset for help on using the changeset viewer.