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 3952 for branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/LBC/lib_mpp.F90 – NEMO

Ignore:
Timestamp:
2013-07-02T15:20:13+02:00 (11 years ago)
Author:
flavoni
Message:

ok for LIM2, still reproducibility error on LIM3 case

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/LBC/lib_mpp.F90

    r3938 r3952  
    362362 
    363363 
    364    SUBROUTINE mpp_lnk_3d( ptab, cd_type, psgn, cd_mpp, pval, ldup ) 
     364   SUBROUTINE mpp_lnk_3d( ptab, cd_type, psgn, cd_mpp, pval ) 
    365365      !!---------------------------------------------------------------------- 
    366366      !!                  ***  routine mpp_lnk_3d  *** 
     
    390390      CHARACTER(len=3), OPTIONAL      , INTENT(in   ) ::   cd_mpp   ! fill the overlap area only  
    391391      REAL(wp)        , OPTIONAL      , INTENT(in   ) ::   pval     ! background value (used at closed boundaries) 
    392       LOGICAL         , OPTIONAL      , INTENT(in   ) ::   ldup     ! duplicate value (used at closed boundaries) 
    393392      !! 
    394393      INTEGER  ::   ji, jj, jk, jl             ! dummy loop indices 
     
    428427            ptab( 1 ,:,:) = ptab(jpim1,:,:) 
    429428            ptab(jpi,:,:) = ptab(  2  ,:,:) 
    430          ELSEIF ( nperio == 0 .AND. PRESENT( ldup ) ) THEN 
    431             ptab(3,:,:) = ptab(2,:,:) 
    432             ptab(:,3,:) = ptab(:,2,:) 
    433             ptab(nlci-2,:,:) = ptab(nlci-1,:,:) 
    434             ptab(:,nlcj-2,:) = ptab(:,nlcj-1,:) 
    435             
    436429         ELSE                                     !* closed 
    437430            IF( .NOT. cd_type == 'F' )   ptab(     1       :jpreci,:,:) = zland    ! south except F-point 
     
    566559 
    567560 
    568    SUBROUTINE mpp_lnk_2d( pt2d, cd_type, psgn, cd_mpp, pval, ldup ) 
     561   SUBROUTINE mpp_lnk_2d( pt2d, cd_type, psgn, cd_mpp, pval ) 
    569562      !!---------------------------------------------------------------------- 
    570563      !!                  ***  routine mpp_lnk_2d  *** 
     
    592585      CHARACTER(len=3), OPTIONAL  , INTENT(in   ) ::   cd_mpp   ! fill the overlap area only  
    593586      REAL(wp)        , OPTIONAL  , INTENT(in   ) ::   pval     ! background value (used at closed boundaries) 
    594       LOGICAL         , OPTIONAL  , INTENT(in   ) ::   ldup     ! duplicate value (used at closed boundaries) 
    595587      !! 
    596588      INTEGER  ::   ji, jj, jl   ! dummy loop indices 
     
    629621            pt2d( 1 ,:) = pt2d(jpim1,:)                                    ! west 
    630622            pt2d(jpi,:) = pt2d(  2  ,:)                                    ! east 
    631          ELSEIF ( nperio == 0 .AND. PRESENT( ldup ) ) THEN 
    632             pt2d(3,:) = pt2d(2,:) 
    633             pt2d(:,3) = pt2d(:,2) 
    634             pt2d(nlci-2,:) = pt2d(nlci-1,:) 
    635             pt2d(:,nlcj-2) = pt2d(:,nlcj-1) 
    636623         ELSE                                     ! closed 
    637624            IF( .NOT. cd_type == 'F' )   pt2d(     1       :jpreci,:) = zland    ! south except F-point 
     
    20242011      CALL MPI_ALLREDUCE( kice, zwork, jpnij, mpi_integer, mpi_sum, mpi_comm_opa, ierr ) 
    20252012      ndim_rank_ice = SUM( zwork )           
    2026   
    2027       !SF ! if there is no ice in the domain, return to the main program (clem modif) 
    2028      !SF  IF ( ndim_rank_ice == 0 ) THEN 
    2029       !SF    DEALLOCATE(kice, zwork) 
    2030      !SF     RETURN 
    2031      !SF  ENDIF 
    20322013 
    20332014      ! Allocate the right size to nrank_north 
     
    20462027      ! Create the world group 
    20472028      CALL MPI_COMM_GROUP( mpi_comm_opa, ngrp_world, ierr ) 
    2048   
     2029 
    20492030      ! Create the ice group from the world group 
    20502031      CALL MPI_GROUP_INCL( ngrp_world, ndim_rank_ice, nrank_ice, ngrp_ice, ierr ) 
Note: See TracChangeset for help on using the changeset viewer.