Ignore:
Timestamp:
2019-05-23T17:13:38+02:00 (2 years ago)
Author:
girrmann
Message:

dev_r10984_HPC-13 : missing part of [11024], see #2285

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/BDY/bdyini.F90

    r11024 r11044  
    12461246      END DO 
    12471247       
    1248       ! detect corners and their orientation index 1 to 4 depending on the orientation 
    1249       ! detect geometries with 3 neighbours  index 5 to 8 depending on the orientation 
    1250       ! else                                 index 0 
     1248      ! detect corner interior and its orientation index 1 to 4  depending on the orientation 
     1249      ! detect corner exterior and its orientation index 5 to 8  depending on the orientation 
     1250      ! detect geometries with 3 neighbours        index 9 to 12 depending on the orientation 
     1251      ! else                                       index 0 
    12511252      DO ib_bdy = 1, nb_bdy 
    12521253         DO igrd = 1, jpbgrd 
     
    12591260               ii        =  idx_bdy(ib_bdy)%nbi(ib,igrd) 
    12601261               ij        =  idx_bdy(ib_bdy)%nbj(ib,igrd) 
    1261                !IF( ii == 1 .OR. ii == jpi .OR. ij == 1 .OR. ij == jpj )   CYCLE 
    12621262               llnobdy = pmask(ii  ,ij+1) == 1.   
    12631263               llsobdy = pmask(ii  ,ij-1) == 1.  
     
    12781278               END IF 
    12791279               IF( inbdy == 2 )   THEN   ! exterior of a corner 
    1280                   idx_bdy(ib_bdy)%ntreat(ib,igrd) = 0   ! regular treatment with flags 
     1280                  !        o      !        o      !    _____|       !       |_____   
     1281                  !  5 ____x o    !  6   o x___   ! 7      x o      !  8   o x       
     1282                  !         |     !       |       !        o        !        o  
     1283                  IF( llnobdy .AND. lleabdy )   idx_bdy(ib_bdy)%ntreat(ib,igrd) = 5 
     1284                  IF( llnobdy .AND. llwebdy )   idx_bdy(ib_bdy)%ntreat(ib,igrd) = 6 
     1285                  IF( llsobdy .AND. lleabdy )   idx_bdy(ib_bdy)%ntreat(ib,igrd) = 7 
     1286                  IF( llsobdy .AND. llwebdy )   idx_bdy(ib_bdy)%ntreat(ib,igrd) = 8 
    12811287               END IF 
    12821288               IF( inbdy == 3 )   THEN   ! 3 neighbours __   __ 
    12831289                  !    |_  o      !        o  _|  !       |_|     !       o          
    1284                   !  5  _| x o    !  6   o x |_   !  7   o x o    ! 8   o x o        
    1285                   !    |   o      !        o   |  !        o      !    __|¨|__       
    1286                   IF( llnobdy .AND. lleabdy .AND. llsobdy )   idx_bdy(ib_bdy)%ntreat(ib,igrd) = 5 
    1287                   IF( llnobdy .AND. llwebdy .AND. llsobdy )   idx_bdy(ib_bdy)%ntreat(ib,igrd) = 6 
    1288                   IF( llwebdy .AND. llsobdy .AND. lleabdy )   idx_bdy(ib_bdy)%ntreat(ib,igrd) = 7 
    1289                   IF( llwebdy .AND. llnobdy .AND. lleabdy )   idx_bdy(ib_bdy)%ntreat(ib,igrd) = 8 
     1290                  !  9  _| x o    ! 10   o x |_   ! 11   o x o    ! 12  o x o        
     1291                  !    |   o      !        o   |  !        o      !    __|¨|__     
     1292                  IF( llnobdy .AND. lleabdy .AND. llsobdy )   idx_bdy(ib_bdy)%ntreat(ib,igrd) = 9 
     1293                  IF( llnobdy .AND. llwebdy .AND. llsobdy )   idx_bdy(ib_bdy)%ntreat(ib,igrd) = 10 
     1294                  IF( llwebdy .AND. llsobdy .AND. lleabdy )   idx_bdy(ib_bdy)%ntreat(ib,igrd) = 11 
     1295                  IF( llwebdy .AND. llnobdy .AND. lleabdy )   idx_bdy(ib_bdy)%ntreat(ib,igrd) = 12 
    12901296               END IF 
    12911297               IF( inbdy == 4 )   THEN 
     
    12971303            END DO 
    12981304         END DO 
    1299          !CALL lbc_lnk( 'bdyini', )  
    13001305      END DO 
    13011306      ! 
Note: See TracChangeset for help on using the changeset viewer.