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 9094 – NEMO

Changeset 9094


Ignore:
Timestamp:
2017-12-15T18:26:58+01:00 (5 years ago)
Author:
cetlod
Message:

Use of lbclnk_multi in subdir LDF & TRA

Location:
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC
Files:
16 edited

Legend:

Unmodified
Added
Removed

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/LDF/ldfc1d_c2d.F90

    r7753 r9094  
    101101            END DO 
    102102         END DO 
    103          CALL lbc_lnk( pah1, 'U', 1. )   ! Lateral boundary conditions 
    104          CALL lbc_lnk( pah2, 'V', 1. )    
     103         ! Lateral boundary conditions 
     104         CALL lbc_lnk_multi( pah1, 'U', 1. , pah2, 'V', 1. )    
    105105         ! 
    106106      END SELECT 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90

    r9019 r9094  
    321321         ENDIF 
    322322         ! 
    323          CALL lbc_lnk( ahmt, 'T', 1. )   ;   CALL lbc_lnk( ahmf, 'F', 1. ) 
     323         CALL lbc_lnk_multi( ahmt, 'T', 1., ahmf, 'F', 1. ) 
    324324         ! 
    325325         ! 
    … …  
    403403         ENDIF 
    404404         ! 
    405          CALL lbc_lnk( ahmt, 'T', 1. )   ;   CALL lbc_lnk( ahmf, 'F', 1. ) 
     405         CALL lbc_lnk_multi( ahmt, 'T', 1. , ahmf, 'F', 1. ) 
    406406         ! 
    407407      END SELECT 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/LDF/ldfslp.F90

    r9019 r9094  
    231231         END DO 
    232232      END DO 
    233       CALL lbc_lnk( zwz, 'U', -1. )   ;   CALL lbc_lnk( zww, 'V', -1. )      ! lateral boundary conditions 
     233      CALL lbc_lnk_multi( zwz, 'U', -1., zww, 'V', -1. )      ! lateral boundary conditions 
    234234      ! 
    235235      !                                            !* horizontal Shapiro filter 
    … …  
    313313         END DO 
    314314      END DO 
    315       CALL lbc_lnk( zwz, 'T', -1. )   ;    CALL lbc_lnk( zww, 'T', -1. )      ! lateral boundary conditions 
     315      CALL lbc_lnk_multi( zwz, 'T', -1., zww, 'T', -1. )      ! lateral boundary conditions 
    316316      ! 
    317317      !                                           !* horizontal Shapiro filter 
    … …  
    362362      ! IV. Lateral boundary conditions 
    363363      ! =============================== 
    364       CALL lbc_lnk( uslp , 'U', -1. )      ;      CALL lbc_lnk( vslp , 'V', -1. ) 
    365       CALL lbc_lnk( wslpi, 'W', -1. )      ;      CALL lbc_lnk( wslpj, 'W', -1. ) 
     364      CALL lbc_lnk_multi( uslp , 'U', -1. , vslp , 'V', -1. , wslpi, 'W', -1., wslpj, 'W', -1. ) 
    366365 
    367366      IF(ln_ctl) THEN 
    … …  
    726725      END DO 
    727726      !!gm this lbc_lnk should be useless.... 
    728       CALL lbc_lnk( uslpml , 'U', -1. )   ;   CALL lbc_lnk( vslpml , 'V', -1. )   ! lateral boundary cond. (sign change) 
    729       CALL lbc_lnk( wslpiml, 'W', -1. )   ;   CALL lbc_lnk( wslpjml, 'W', -1. )   ! lateral boundary conditions 
     727      CALL lbc_lnk_multi( uslpml , 'U', -1. , vslpml , 'V', -1. , wslpiml, 'W', -1. , wslpjml, 'W', -1. )  
    730728      ! 
    731729      IF( ln_timing )   CALL timing_stop('ldf_slp_mxl') 
    … …  
    799797!               END DO 
    800798!            END DO 
    801 !            CALL lbc_lnk( uslp , 'U', -1. )   ;   CALL lbc_lnk( vslp , 'V', -1. )      ! Lateral boundary conditions 
    802 !            CALL lbc_lnk( wslpi, 'W', -1. )   ;   CALL lbc_lnk( wslpj, 'W', -1. ) 
     799!            CALL lbc_lnk_multi( uslp , 'U', -1. ; CALL lbc_lnk( vslp , 'V', -1.,  wslpi, 'W', -1.,  wslpj, 'W', -1. ) 
    803800!!gm         ENDIF 
    804801      ENDIF 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90

    r9019 r9094  
    564564         END DO  
    565565      END DO  
    566       CALL lbc_lnk( paeiu(:,:,1), 'U', 1. )   ;   CALL lbc_lnk( paeiv(:,:,1), 'V', 1. )      ! lateral boundary condition 
     566      CALL lbc_lnk_multi( paeiu(:,:,1), 'U', 1. , paeiv(:,:,1), 'V', 1. )      ! lateral boundary condition 
    567567 
    568568      DO jk = 2, jpkm1                          !==  deeper values equal the surface one  ==! 
    … …  
    678678      ! 
    679679      !                                                  !==  eiv stream function: output  ==! 
    680       CALL lbc_lnk( psi_uw, 'U', -1. )                         ! lateral boundary condition 
    681       CALL lbc_lnk( psi_vw, 'V', -1. ) 
     680      CALL lbc_lnk_multi( psi_uw, 'U', -1. , psi_vw, 'V', -1. ) 
    682681      ! 
    683682!!gm      CALL iom_put( "psi_eiv_uw", psi_uw )                 ! output 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90

    r9019 r9094  
    755755            END DO 
    756756         END DO 
    757          ! 
    758          CALL lbc_lnk( pab(:,:,jp_tem), 'T', 1. )                    ! Lateral boundary conditions 
    759          CALL lbc_lnk( pab(:,:,jp_sal), 'T', 1. )                     
     757         !                            ! Lateral boundary conditions 
     758         CALL lbc_lnk_multi( pab(:,:,jp_tem), 'T', 1. , pab(:,:,jp_sal), 'T', 1. )                     
    760759         ! 
    761760      CASE( np_seos )                  !==  simplified EOS  ==! 
    … …  
    776775            END DO 
    777776         END DO 
    778          ! 
    779          CALL lbc_lnk( pab(:,:,jp_tem), 'T', 1. )                    ! Lateral boundary conditions 
    780          CALL lbc_lnk( pab(:,:,jp_sal), 'T', 1. )                     
     777         !                            ! Lateral boundary conditions 
     778         CALL lbc_lnk_multi( pab(:,:,jp_tem), 'T', 1. , pab(:,:,jp_sal), 'T', 1. )                     
    781779         ! 
    782780      CASE DEFAULT 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_cen.F90

    r9019 r9094  
    126126               END DO 
    127127            END DO 
    128             CALL lbc_lnk( ztu, 'U', -1. )   ;    CALL lbc_lnk( ztv, 'V', -1. )   ! Lateral boundary cond. 
     128            CALL lbc_lnk_multi( ztu, 'U', -1. , ztv, 'V', -1. )   ! Lateral boundary cond. 
    129129            ! 
    130130            DO jk = 1, jpkm1                       ! Horizontal advective fluxes 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_fct.F90

    r9019 r9094  
    211211               END DO 
    212212            END DO 
    213             CALL lbc_lnk( zltu, 'T', 1. )   ;    CALL lbc_lnk( zltv, 'T', 1. )   ! Lateral boundary cond. (unchanged sgn) 
     213            CALL lbc_lnk_multi( zltu, 'T', 1. , zltv, 'T', 1. )   ! Lateral boundary cond. (unchanged sgn) 
    214214            ! 
    215215            DO jk = 1, jpkm1                 ! Horizontal advective fluxes 
    … …  
    236236               END DO 
    237237            END DO 
    238             CALL lbc_lnk( ztu, 'U', -1. )   ;    CALL lbc_lnk( ztv, 'V', -1. )   ! Lateral boundary cond. (unchanged sgn) 
     238            CALL lbc_lnk_multi( ztu, 'U', -1. , ztv, 'V', -1. )   ! Lateral boundary cond. (unchanged sgn) 
    239239            ! 
    240240            DO jk = 1, jpkm1                 ! Horizontal advective fluxes 
    … …  
    282282         ENDIF 
    283283         ! 
    284          CALL lbc_lnk( zwx, 'U', -1. )   ;   CALL lbc_lnk( zwy, 'V', -1. )         ! Lateral bondary conditions 
    285          CALL lbc_lnk( zwz, 'W',  1. ) 
     284         CALL lbc_lnk_multi( zwx, 'U', -1. , zwy, 'V', -1.,  zwz, 'W',  1. ) 
    286285         ! 
    287286         !        !==  monotonicity algorithm  ==! 
    … …  
    402401         END DO 
    403402      END DO 
    404       CALL lbc_lnk( zbetup, 'T', 1. )   ;   CALL lbc_lnk( zbetdo, 'T', 1. )   ! lateral boundary cond. (unchanged sign) 
     403      CALL lbc_lnk_multi( zbetup, 'T', 1. , zbetdo, 'T', 1. )   ! lateral boundary cond. (unchanged sign) 
    405404 
    406405      ! 3. monotonic flux in the i & j direction (paa & pbb) 
    … …  
    428427         END DO 
    429428      END DO 
    430       CALL lbc_lnk( paa, 'U', -1. )   ;   CALL lbc_lnk( pbb, 'V', -1. )   ! lateral boundary condition (changed sign) 
     429      CALL lbc_lnk_multi( paa, 'U', -1. , pbb, 'V', -1. )   ! lateral boundary condition (changed sign) 
    431430      ! 
    432431      IF( ln_timing )   CALL timing_stop('nonosc') 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_mle.F90

    r9019 r9094  
    324324               END DO 
    325325            END DO 
    326             CALL lbc_lnk( rfu, 'U', 1. )   ;   CALL lbc_lnk( rfv, 'V', 1. ) 
     326            CALL lbc_lnk_multi( rfu, 'U', 1. , rfv, 'V', 1. ) 
    327327            ! 
    328328         ELSEIF( nn_mle == 1 ) THEN           ! MLE array allocation & initialisation 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_mus.F90

    r9019 r9094  
    142142           END DO 
    143143         END DO 
    144          CALL lbc_lnk( zwx, 'U', -1. )          ! lateral boundary conditions   (changed sign) 
    145          CALL lbc_lnk( zwy, 'V', -1. ) 
     144         ! lateral boundary conditions   (changed sign) 
     145         CALL lbc_lnk_multi( zwx, 'U', -1. , zwy, 'V', -1. ) 
    146146         !                                !-- Slopes of tracer 
    147147         zslpx(:,:,jpk) = 0._wp                 ! bottom values 
    … …  
    191191            END DO 
    192192         END DO 
    193          CALL lbc_lnk( zwx, 'U', -1. )   ;   CALL lbc_lnk( zwy, 'V', -1. )   ! lateral boundary conditions   (changed sign) 
     193         CALL lbc_lnk_multi( zwx, 'U', -1. , zwy, 'V', -1. )   ! lateral boundary conditions   (changed sign) 
    194194         ! 
    195195         DO jk = 1, jpkm1                 !-- Tracer advective trend 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_qck.F90

    r9019 r9094  
    155155            END DO 
    156156         END DO 
    157          CALL lbc_lnk( zfc(:,:,:), 'T', 1. )   ;   CALL lbc_lnk( zfd(:,:,:), 'T', 1. )   ! Lateral boundary conditions  
     157         CALL lbc_lnk_multi( zfc(:,:,:), 'T', 1. , zfd(:,:,:), 'T', 1. )   ! Lateral boundary conditions  
    158158          
    159159         ! 
    … …  
    181181         END DO  
    182182         !--- Lateral boundary conditions  
    183          CALL lbc_lnk( zfu(:,:,:), 'T', 1. )   ;   CALL lbc_lnk( zfd(:,:,:), 'T', 1. ) 
    184          CALL lbc_lnk( zfc(:,:,:), 'T', 1. )   ;   CALL lbc_lnk( zwx(:,:,:), 'T', 1. ) 
     183         CALL lbc_lnk_multi( zfu(:,:,:), 'T', 1. , zfd(:,:,:), 'T', 1., zfc(:,:,:), 'T', 1.,  zwx(:,:,:), 'T', 1. ) 
    185184 
    186185         !--- QUICKEST scheme 
    … …  
    271270            END DO 
    272271         END DO 
    273          CALL lbc_lnk( zfc(:,:,:), 'T', 1. )   ;   CALL lbc_lnk( zfd(:,:,:), 'T', 1. )   ! Lateral boundary conditions  
     272         CALL lbc_lnk_multi( zfc(:,:,:), 'T', 1. , zfd(:,:,:), 'T', 1. )   ! Lateral boundary conditions  
    274273 
    275274          
    … …  
    300299 
    301300         !--- Lateral boundary conditions  
    302          CALL lbc_lnk( zfu(:,:,:), 'T', 1. )      ;     CALL lbc_lnk( zfd(:,:,:), 'T', 1. ) 
    303          CALL lbc_lnk( zfc(:,:,:), 'T', 1. )      ;     CALL lbc_lnk( zwy(:,:,:), 'T', 1. ) 
     301         CALL lbc_lnk_multi( zfu(:,:,:), 'T', 1. , zfd(:,:,:), 'T', 1., zfc(:,:,:), 'T', 1., zwy(:,:,:), 'T', 1. ) 
    304302 
    305303         !--- QUICKEST scheme 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA/trabbl.F90

    r9019 r9094  
    123123            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=           ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    124124         ! lateral boundary conditions ; just need for outputs 
    125          CALL lbc_lnk( ahu_bbl, 'U', 1. )     ;     CALL lbc_lnk( ahv_bbl, 'V', 1. ) 
     125         CALL lbc_lnk_multi( ahu_bbl, 'U', 1. , ahv_bbl, 'V', 1. ) 
    126126         CALL iom_put( "ahu_bbl", ahu_bbl )   ! bbl diffusive flux i-coef 
    127127         CALL iom_put( "ahv_bbl", ahv_bbl )   ! bbl diffusive flux j-coef 
    … …  
    136136            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=           ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    137137         ! lateral boundary conditions ; just need for outputs 
    138          CALL lbc_lnk( utr_bbl, 'U', 1. )     ;   CALL lbc_lnk( vtr_bbl, 'V', 1. ) 
     138         CALL lbc_lnk_multi( utr_bbl, 'U', 1. , vtr_bbl, 'V', 1. ) 
    139139         CALL iom_put( "uoce_bbl", utr_bbl )  ! bbl i-transport 
    140140         CALL iom_put( "voce_bbl", vtr_bbl )  ! bbl j-transport 
    … …  
    491491      INTEGER ::   ji, jj                      ! dummy loop indices 
    492492      INTEGER ::   ii0, ii1, ij0, ij1, ios     ! local integer 
    493       REAL(wp), DIMENSION(jpi,jpj) ::   zmbk   ! workspace 
     493      REAL(wp), DIMENSION(jpi,jpj) ::   zmbku, zmbkv   ! workspace 
    494494      !! 
    495495      NAMELIST/nambbl/ ln_trabbl, nn_bbl_ldf, nn_bbl_adv, rn_ahtbbl, rn_gambbl 
    … …  
    539539      END DO 
    540540      ! converte into REAL to use lbc_lnk ; impose a min value of 1 as a zero can be set in lbclnk 
    541       zmbk(:,:) = REAL( mbku_d(:,:), wp )   ;   CALL lbc_lnk(zmbk,'U',1.)   ;   mbku_d(:,:) = MAX( INT( zmbk(:,:) ), 1 ) 
    542       zmbk(:,:) = REAL( mbkv_d(:,:), wp )   ;   CALL lbc_lnk(zmbk,'V',1.)   ;   mbkv_d(:,:) = MAX( INT( zmbk(:,:) ), 1 ) 
     541      zmbku(:,:) = REAL( mbku_d(:,:), wp )   ;     zmbkv(:,:) = REAL( mbkv_d(:,:), wp )   
     542      CALL lbc_lnk_multi( zmbku,'U',1., zmbk,'V',1.)  
     543      mbku_d(:,:) = MAX( INT( zmbku(:,:) ), 1 ) ;  mbkv_d(:,:) = MAX( INT( zmbkv(:,:) ), 1 ) 
    543544      ! 
    544545                                        !* sign of grad(H) at u- and v-points; zero if grad(H) = 0 
    … …  
    562563         END DO 
    563564      END DO 
    564       CALL lbc_lnk( e3u_bbl_0, 'U', 1. )   ;   CALL lbc_lnk( e3v_bbl_0, 'V', 1. )      ! lateral boundary conditions 
     565      CALL lbc_lnk_multi( e3u_bbl_0, 'U', 1. , e3v_bbl_0, 'V', 1. )      ! lateral boundary conditions 
    565566      ! 
    566567      !                             !* masked diffusive flux coefficients 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_lap_blp.F90

    r9019 r9094  
    227227      END SELECT 
    228228      ! 
    229       DO jn = 1, kjpt 
    230          CALL lbc_lnk( zlap(:,:,:,jn) , 'T', 1. )     ! Lateral boundary conditions (unchanged sign) 
    231       END DO 
     229      CALL lbc_lnk( zlap(:,:,:,:) , 'T', 1. )     ! Lateral boundary conditions (unchanged sign) 
    232230      !                                               ! Partial top/bottom cell: GRADh( zlap )   
    233231      IF( ln_isfcav .AND. ln_zps ) THEN   ;   CALL zps_hde_isf( kt, kjpt, zlap, zglu, zglv, zgui, zgvi )  ! both top & bottom 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA/tranpc.F90

    r9019 r9094  
    309309         ENDIF 
    310310         ! 
    311          CALL lbc_lnk( tsa(:,:,:,jp_tem), 'T', 1. )   ;   CALL lbc_lnk( tsa(:,:,:,jp_sal), 'T', 1. ) 
     311         CALL lbc_lnk_multi( tsa(:,:,:,jp_tem), 'T', 1., tsa(:,:,:,jp_sal), 'T', 1. ) 
    312312         ! 
    313313         IF( lwp .AND. l_LB_debug ) THEN 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

    r9023 r9094  
    106106      CALL Agrif_tra                     ! AGRIF zoom boundaries 
    107107#endif 
    108       ! 
    109       CALL lbc_lnk( tsa(:,:,:,jp_tem), 'T', 1._wp )      ! local domain boundaries  (T-point, unchanged sign) 
    110       CALL lbc_lnk( tsa(:,:,:,jp_sal), 'T', 1._wp ) 
     108      !                                              ! local domain boundaries  (T-point, unchanged sign) 
     109      CALL lbc_lnk_multi( tsa(:,:,:,jp_tem), 'T', 1., tsa(:,:,:,jp_sal), 'T', 1. ) 
    111110      ! 
    112111      IF( ln_bdy )   CALL bdy_tra( kt )  ! BDY open boundaries 
    … …  
    164163         ENDIF 
    165164         ! 
    166          DO jn = 1, jpts 
    167             CALL lbc_lnk( tsb(:,:,:,jn), 'T', 1._wp )  
    168             CALL lbc_lnk( tsn(:,:,:,jn), 'T', 1._wp ) 
    169             CALL lbc_lnk( tsa(:,:,:,jn), 'T', 1._wp ) 
    170          END DO 
     165         CALL lbc_lnk_multi( tsb(:,:,:,jp_tem), 'T', 1., tsb(:,:,:,jp_sal), 'T', 1., & 
     166                  &          tsn(:,:,:,jp_tem), 'T', 1., tsn(:,:,:,jp_sal), 'T', 1., & 
     167                  &          tsa(:,:,:,jp_tem), 'T', 1., tsa(:,:,:,jp_sal), 'T', 1.  ) 
     168         ! 
    171169      ENDIF      
    172170      ! 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf.F90

    r9019 r9094  
    8686         END DO 
    8787!!gm this should be moved in trdtra.F90 and done on all trends 
    88          CALL lbc_lnk( ztrdt, 'T', 1. ) 
    89          CALL lbc_lnk( ztrds, 'T', 1. ) 
     88         CALL lbc_lnk_multi( ztrdt, 'T', 1. , ztrds, 'T', 1. ) 
    9089!!gm 
    9190         CALL trd_tra( kt, 'TRA', jp_tem, jptra_zdf, ztrdt ) 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA/zpshde.F90

    r9019 r9094  
    144144            END DO 
    145145         END DO 
    146          CALL lbc_lnk( pgtu(:,:,jn), 'U', -1. )   ;   CALL lbc_lnk( pgtv(:,:,jn), 'V', -1. )   ! Lateral boundary cond. 
     146         CALL lbc_lnk_multi( pgtu(:,:,jn), 'U', -1. , pgtv(:,:,jn), 'V', -1. )   ! Lateral boundary cond. 
    147147         ! 
    148148      END DO 
    … …  
    183183            END DO 
    184184         END DO 
    185          CALL lbc_lnk( pgru , 'U', -1. )   ;   CALL lbc_lnk( pgrv , 'V', -1. )   ! Lateral boundary conditions 
     185         CALL lbc_lnk_multi( pgru , 'U', -1. , pgrv , 'V', -1. )   ! Lateral boundary conditions 
    186186         ! 
    187187      END IF 
    … …  
    305305            END DO 
    306306         END DO 
    307          CALL lbc_lnk( pgtu(:,:,jn), 'U', -1. )   ;   CALL lbc_lnk( pgtv(:,:,jn), 'V', -1. )   ! Lateral boundary cond. 
     307         CALL lbc_lnk_multi( pgtu(:,:,jn), 'U', -1. , pgtv(:,:,jn), 'V', -1. )   ! Lateral boundary cond. 
    308308         ! 
    309309      END DO 
    … …  
    353353         END DO 
    354354 
    355          CALL lbc_lnk( pgru , 'U', -1. )   ;   CALL lbc_lnk( pgrv , 'V', -1. )   ! Lateral boundary conditions 
     355         CALL lbc_lnk_multi( pgru , 'U', -1. , pgrv , 'V', -1. )   ! Lateral boundary conditions 
    356356         ! 
    357357      END IF 
    … …  
    404404            END DO 
    405405         END DO 
    406          CALL lbc_lnk( pgtui(:,:,jn), 'U', -1. ); CALL lbc_lnk( pgtvi(:,:,jn), 'V', -1. )   ! Lateral boundary cond. 
    407406         ! 
    408407      END DO 
     408      CALL lbc_lnk_multi( pgtui(:,:,:), 'U', -1. , pgtvi(:,:,:), 'V', -1. )   ! Lateral boundary cond. 
    409409 
    410410      IF( PRESENT( prd ) ) THEN    !==  horizontal derivative of density anomalies (rd)  ==!    (optional part) 
    … …  
    449449            END DO 
    450450         END DO 
    451          CALL lbc_lnk( pgrui   , 'U', -1. ); CALL lbc_lnk( pgrvi   , 'V', -1. )   ! Lateral boundary conditions 
     451         CALL lbc_lnk_multi( pgrui, 'U', -1. , pgrvi, 'V', -1. )   ! Lateral boundary conditions 
    452452         ! 
    453453      END IF   
Note: See TracChangeset for help on using the changeset viewer.