New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 415 – NEMO

Changeset 415


Ignore:
Timestamp:
2006-03-20T18:35:34+01:00 (18 years ago)
Author:
opalod
Message:

nemo_v1_bugfix_033 : CT : correct the index used for sea-ice North-fold (jpiglo instead of jpi)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/lib_mpp.F90

    r392 r415  
    5353   IMPLICIT NONE 
    5454 
     55   PRIVATE 
     56   PUBLIC  mynode, mpparent, mpp_isl, mpp_min, mpp_max, mpp_sum,  mpp_lbc_north 
     57   PUBLIC  mpp_lbc_north_e, mpp_minloc, mpp_maxloc, mpp_lnk_3d, mpp_lnk_2d, mpp_lnk_2d_e, mpplnks 
     58   PUBLIC  mpprecv, mppsend, mppscatter, mppgather, mppobc, mpp_ini_north, mppstop, mppsync 
     59 
    5560   !! * Interfaces 
    5661   !! define generic interface for these routine as they are called sometimes 
     
    285290      WRITE(numout,*) '           mpi send type            c_mpi_send = ', c_mpi_send 
    286291 
    287       SELECT CASE ( c_mpi_send ) 
    288       CASE ( 'S' )                ! Standard mpi send (blocking) 
    289          WRITE(numout,*) '           Standard blocking mpi send (send)' 
    290292#if defined key_agrif 
    291          IF ( Agrif_Root() ) THEN 
    292 #endif       
    293          CALL mpi_init( ierr ) 
     293      IF( Agrif_Root() ) THEN 
     294#endif 
     295         SELECT CASE ( c_mpi_send ) 
     296         CASE ( 'S' )                ! Standard mpi send (blocking) 
     297            WRITE(numout,*) '           Standard blocking mpi send (send)' 
     298            CALL mpi_init( ierr ) 
     299         CASE ( 'B' )                ! Buffer mpi send (blocking) 
     300            WRITE(numout,*) '           Buffer blocking mpi send (bsend)' 
     301            CALL mpi_init_opa( ierr ) 
     302         CASE ( 'I' )                ! Immediate mpi send (non-blocking send) 
     303            WRITE(numout,*) '           Immediate non-blocking send (isend)' 
     304            l_isend = .TRUE. 
     305            CALL mpi_init( ierr ) 
     306         CASE DEFAULT 
     307            WRITE(numout,cform_err) 
     308            WRITE(numout,*) '           bad value for c_mpi_send = ', c_mpi_send 
     309            nstop = nstop + 1 
     310         END SELECT 
     311 
    294312#if defined key_agrif 
    295          ENDIF 
    296 #endif     
    297       CASE ( 'B' )                ! Buffer mpi send (blocking) 
    298          WRITE(numout,*) '           Buffer blocking mpi send (bsend)' 
    299 #if defined key_agrif 
    300          IF ( Agrif_Root() ) THEN 
    301 #endif       
    302          CALL mpi_init_opa( ierr ) 
    303 #if defined key_agrif 
    304          ENDIF 
    305 #endif     
    306       CASE ( 'I' )                ! Immediate mpi send (non-blocking send) 
    307          WRITE(numout,*) '           Immediate non-blocking send (isend)' 
    308          l_isend = .TRUE. 
    309 #if defined key_agrif 
    310          IF ( Agrif_Root() ) THEN 
    311 #endif       
    312          CALL mpi_init( ierr ) 
    313 #if defined key_agrif 
    314          ENDIF 
    315 #endif     
    316       CASE DEFAULT 
    317          WRITE(numout,cform_err) 
    318          WRITE(numout,*) '           bad value for c_mpi_send = ', c_mpi_send 
    319          nstop = nstop + 1 
    320       END SELECT 
     313      ENDIF 
     314#endif 
    321315 
    322316      CALL mpi_comm_rank( mpi_comm_world, rank, ierr ) 
     
    40734067                ztab( 2 ,ijpj) = 0.e0 
    40744068                DO ji = 2 , jpiglo-1 
    4075                    ijt = jpi - ji + 2 
     4069                   ijt = jpiglo - ji + 2 
    40764070                   ztab(ji,ijpj)= 0.5 * ( ztab(ji,ijpj-1) + psgn * ztab(ijt,ijpj-1) ) 
    40774071                END DO 
     
    43374331                DO jl = 0, jpr2dj 
    43384332                   DO ji = 2 , jpiglo-1 
    4339                       ijt = jpi - ji + 2 
     4333                      ijt = jpiglo - ji + 2 
    43404334                      ztab(ji,ijpj+jl)= 0.5 * ( ztab(ji,ijpj-1-jl) + psgn * ztab(ijt,ijpj-1-jl) ) 
    43414335                   END DO 
Note: See TracChangeset for help on using the changeset viewer.