Changeset 13275


Ignore:
Timestamp:
2020-07-08T19:36:41+02:00 (4 weeks ago)
Author:
smasson
Message:

Extra_Halo: final bugfix on bestpartition when nn_hls > 1, see #2366

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE/LBC/mppini.F90

    r13269 r13275  
    766766         DO jj = 1, knbj 
    767767            DO ji = 2, knbi 
    768                kimppt(ji,jj) = kimppt(ji-1,jj) + klci(ji-1,jj) - 2 * nn_hls 
     768               kimppt(ji,jj) = kimppt(ji-1,jj) + klci(ji-1,jj) - i2hls 
    769769            END DO 
    770770         END DO 
     
    774774         DO jj = 2, knbj 
    775775            DO ji = 1, knbi 
    776                kjmppt(ji,jj) = kjmppt(ji,jj-1) + klcj(ji,jj-1) - 2 * nn_hls 
     776               kjmppt(ji,jj) = kjmppt(ji,jj-1) + klcj(ji,jj-1) - i2hls 
    777777            END DO 
    778778         END DO 
     
    10961096            IF( iarea == inbj ) THEN                                   ! the last line was not read 
    10971097               IF( jperio == 2 .OR. jperio == 7 ) THEN                 !   north-south periodocity 
    1098                   CALL readbot_strip( 1, 1, lloce(2:inx-1,iny) )       !   read the first line -> last line of lloce 
    1099                ELSE 
    1100                   lloce(2:inx-1,iny) = .FALSE.                         !   closed boundary 
     1098                  CALL readbot_strip( 1, 1, lloce(2:inx-1,iny) )       !      read the first line -> last line of lloce 
     1099               ELSEIF( jperio == 3 .OR. jperio == 4 ) THEN             !   north-pole folding T-pivot, T-point  
     1100                  lloce(2,iny) = lloce(2,iny-2)                        !      here we have 1 halo (even if nn_hls>1) 
     1101                  DO ji = 3,inx-1 
     1102                     lloce(ji,iny  ) = lloce(inx-ji+2,iny-2)           !      ok, we have at least 3 lines 
     1103                  END DO 
     1104                  DO ji = inx/2+2,inx-1 
     1105                     lloce(ji,iny-1) = lloce(inx-ji+2,iny-1) 
     1106                  END DO 
     1107               ELSEIF( jperio == 5 .OR. jperio == 6 ) THEN             !   north-pole folding F-pivot, T-point, 1 halo 
     1108                  lloce(inx/2+1,iny-1) = lloce(inx/2,iny-1)            !      here we have 1 halo (even if nn_hls>1) 
     1109                  lloce(inx  -1,iny-1) = lloce(2    ,iny-1) 
     1110                  DO ji = 2,inx-1 
     1111                     lloce(ji,iny) = lloce(inx-ji+1,iny-1) 
     1112                  END DO 
     1113               ELSE                                                    !   closed boundary 
     1114                  lloce(2:inx-1,iny) = .FALSE. 
    11011115               ENDIF 
    11021116            ENDIF 
Note: See TracChangeset for help on using the changeset viewer.