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 10425 for NEMO/trunk/src/OCE/TRA – NEMO

Ignore:
Timestamp:
2018-12-19T22:54:16+01:00 (5 years ago)
Author:
smasson
Message:

trunk: merge back dev_r10164_HPC09_ESIWACE_PREP_MERGE@10424 into the trunk

Location:
NEMO/trunk/src/OCE/TRA
Files:
17 edited

Legend:

Unmodified
Added
Removed

  • NEMO/trunk/src/OCE/TRA/eosbn2.F90

    r10068 r10425  
    533533         END DO 
    534534         ! 
    535          CALL lbc_lnk( prd, 'T', 1. )                    ! Lateral boundary conditions 
     535         CALL lbc_lnk( 'eosbn2', prd, 'T', 1. )                    ! Lateral boundary conditions 
    536536         ! 
    537537      CASE( np_seos )                !==  simplified EOS  ==! 
    … …  
    553553         END DO 
    554554         ! 
    555          CALL lbc_lnk( prd, 'T', 1. )                    ! Lateral boundary conditions 
     555         CALL lbc_lnk( 'eosbn2', prd, 'T', 1. )                    ! Lateral boundary conditions 
    556556         ! 
    557557      END SELECT 
    … …  
    664664         ! 
    665665      CASE DEFAULT 
    666          IF(lwp) WRITE(numout,cform_err) 
    667          IF(lwp) WRITE(numout,*) '          bad flag value for neos = ', neos 
    668          nstop = nstop + 1 
     666         WRITE(ctmp1,*) '          bad flag value for neos = ', neos 
     667         CALL ctl_stop( 'rab_3d:', ctmp1 ) 
    669668         ! 
    670669      END SELECT 
    … …  
    756755         END DO 
    757756         !                            ! Lateral boundary conditions 
    758          CALL lbc_lnk_multi( pab(:,:,jp_tem), 'T', 1. , pab(:,:,jp_sal), 'T', 1. )                     
     757         CALL lbc_lnk_multi( 'eosbn2', pab(:,:,jp_tem), 'T', 1. , pab(:,:,jp_sal), 'T', 1. )                     
    759758         ! 
    760759      CASE( np_seos )                  !==  simplified EOS  ==! 
    … …  
    776775         END DO 
    777776         !                            ! Lateral boundary conditions 
    778          CALL lbc_lnk_multi( pab(:,:,jp_tem), 'T', 1. , pab(:,:,jp_sal), 'T', 1. )                     
     777         CALL lbc_lnk_multi( 'eosbn2', pab(:,:,jp_tem), 'T', 1. , pab(:,:,jp_sal), 'T', 1. )                     
    779778         ! 
    780779      CASE DEFAULT 
    781          IF(lwp) WRITE(numout,cform_err) 
    782          IF(lwp) WRITE(numout,*) '          bad flag value for neos = ', neos 
    783          nstop = nstop + 1 
     780         WRITE(ctmp1,*) '          bad flag value for neos = ', neos 
     781         CALL ctl_stop( 'rab_2d:', ctmp1 ) 
    784782         ! 
    785783      END SELECT 
    … …  
    879877         ! 
    880878      CASE DEFAULT 
    881          IF(lwp) WRITE(numout,cform_err) 
    882          IF(lwp) WRITE(numout,*) '          bad flag value for neos = ', neos 
    883          nstop = nstop + 1 
     879         WRITE(ctmp1,*) '          bad flag value for neos = ', neos 
     880         CALL ctl_stop( 'rab_0d:', ctmp1 ) 
    884881         ! 
    885882      END SELECT 
    … …  
    10401037         ! 
    10411038      CASE DEFAULT 
    1042          IF(lwp) WRITE(numout,cform_err) 
    1043          IF(lwp) WRITE(numout,*) '          bad flag value for neos = ', neos 
    1044          nstop = nstop + 1 
     1039         WRITE(ctmp1,*) '          bad flag value for neos = ', neos 
     1040         CALL ctl_stop( 'eos_fzp_2d:', ctmp1 ) 
    10451041         ! 
    10461042      END SELECT       
    … …  
    10871083         ! 
    10881084      CASE DEFAULT 
    1089          IF(lwp) WRITE(numout,cform_err) 
    1090          IF(lwp) WRITE(numout,*) '          bad flag value for neos = ', neos 
    1091          nstop = nstop + 1 
     1085         WRITE(ctmp1,*) '          bad flag value for neos = ', neos 
     1086         CALL ctl_stop( 'eos_fzp_0d:', ctmp1 ) 
    10921087         ! 
    10931088      END SELECT 
    … …  
    12131208         ! 
    12141209      CASE DEFAULT 
    1215          IF(lwp) WRITE(numout,cform_err) 
    1216          IF(lwp) WRITE(numout,*) '          bad flag value for neos = ', neos 
    1217          nstop = nstop + 1 
     1210         WRITE(ctmp1,*) '          bad flag value for neos = ', neos 
     1211         CALL ctl_stop( 'eos_pen:', ctmp1 ) 
    12181212         ! 
    12191213      END SELECT 

  • NEMO/trunk/src/OCE/TRA/traadv_cen.F90

    r10069 r10425  
    123123               END DO 
    124124            END DO 
    125             CALL lbc_lnk_multi( ztu, 'U', -1. , ztv, 'V', -1. )   ! Lateral boundary cond. 
     125            CALL lbc_lnk_multi( 'traadv_cen', ztu, 'U', -1. , ztv, 'V', -1. )   ! Lateral boundary cond. 
    126126            ! 
    127127            DO jk = 1, jpkm1                       ! Horizontal advective fluxes 

  • NEMO/trunk/src/OCE/TRA/traadv_fct.F90

    r10068 r10425  
    169169            END DO 
    170170         END DO 
    171          CALL lbc_lnk( zwi, 'T', 1. )  ! Lateral boundary conditions on zwi  (unchanged sign) 
    172171         !                 
    173172         IF( l_trd .OR. l_hst )  THEN             ! trend diagnostics (contribution of upstream fluxes) 
    … …  
    208207               END DO 
    209208            END DO 
    210             CALL lbc_lnk_multi( zltu, 'T', 1. , zltv, 'T', 1. )   ! Lateral boundary cond. (unchanged sgn) 
     209            CALL lbc_lnk_multi( 'traadv_fct', zltu, 'T', 1. , zltv, 'T', 1. )   ! Lateral boundary cond. (unchanged sgn) 
    211210            ! 
    212211            DO jk = 1, jpkm1                 ! Horizontal advective fluxes 
    … …  
    233232               END DO 
    234233            END DO 
    235             CALL lbc_lnk_multi( ztu, 'U', -1. , ztv, 'V', -1. )   ! Lateral boundary cond. (unchanged sgn) 
     234            CALL lbc_lnk_multi( 'traadv_fct', ztu, 'U', -1. , ztv, 'V', -1. )   ! Lateral boundary cond. (unchanged sgn) 
    236235            ! 
    237236            DO jk = 1, jpkm1                 ! Horizontal advective fluxes 
    … …  
    279278         ENDIF 
    280279         ! 
    281          CALL lbc_lnk_multi( zwx, 'U', -1. , zwy, 'V', -1.,  zwz, 'W',  1. ) 
     280         CALL lbc_lnk_multi( 'traadv_fct', zwi, 'T', 1., zwx, 'U', -1. , zwy, 'V', -1.,  zwz, 'W',  1. ) 
    282281         ! 
    283282         !        !==  monotonicity algorithm  ==! 
    … …  
    399398         END DO 
    400399      END DO 
    401       CALL lbc_lnk_multi( zbetup, 'T', 1. , zbetdo, 'T', 1. )   ! lateral boundary cond. (unchanged sign) 
     400      CALL lbc_lnk_multi( 'traadv_fct', zbetup, 'T', 1. , zbetdo, 'T', 1. )   ! lateral boundary cond. (unchanged sign) 
    402401 
    403402      ! 3. monotonic flux in the i & j direction (paa & pbb) 
    … …  
    425424         END DO 
    426425      END DO 
    427       CALL lbc_lnk_multi( paa, 'U', -1. , pbb, 'V', -1. )   ! lateral boundary condition (changed sign) 
     426      CALL lbc_lnk_multi( 'traadv_fct', paa, 'U', -1. , pbb, 'V', -1. )   ! lateral boundary condition (changed sign) 
    428427      ! 
    429428   END SUBROUTINE nonosc 

  • NEMO/trunk/src/OCE/TRA/traadv_mus.F90

    r10068 r10425  
    140140         END DO 
    141141         ! lateral boundary conditions   (changed sign) 
    142          CALL lbc_lnk_multi( zwx, 'U', -1. , zwy, 'V', -1. ) 
     142         CALL lbc_lnk_multi( 'traadv_mus', zwx, 'U', -1. , zwy, 'V', -1. ) 
    143143         !                                !-- Slopes of tracer 
    144144         zslpx(:,:,jpk) = 0._wp                 ! bottom values 
    … …  
    188188            END DO 
    189189         END DO 
    190          CALL lbc_lnk_multi( zwx, 'U', -1. , zwy, 'V', -1. )   ! lateral boundary conditions   (changed sign) 
     190         CALL lbc_lnk_multi( 'traadv_mus', zwx, 'U', -1. , zwy, 'V', -1. )   ! lateral boundary conditions   (changed sign) 
    191191         ! 
    192192         DO jk = 1, jpkm1                 !-- Tracer advective trend 

  • NEMO/trunk/src/OCE/TRA/traadv_qck.F90

    r10068 r10425  
    150150            END DO 
    151151         END DO 
    152          CALL lbc_lnk_multi( zfc(:,:,:), 'T', 1. , zfd(:,:,:), 'T', 1. )   ! Lateral boundary conditions  
     152         CALL lbc_lnk_multi( 'traadv_qck', zfc(:,:,:), 'T', 1. , zfd(:,:,:), 'T', 1. )   ! Lateral boundary conditions  
    153153          
    154154         ! 
    … …  
    176176         END DO  
    177177         !--- Lateral boundary conditions  
    178          CALL lbc_lnk_multi( zfu(:,:,:), 'T', 1. , zfd(:,:,:), 'T', 1., zfc(:,:,:), 'T', 1.,  zwx(:,:,:), 'T', 1. ) 
     178         CALL lbc_lnk_multi( 'traadv_qck', zfu(:,:,:), 'T', 1. , zfd(:,:,:), 'T', 1., zfc(:,:,:), 'T', 1.,  zwx(:,:,:), 'T', 1. ) 
    179179 
    180180         !--- QUICKEST scheme 
    … …  
    189189            END DO 
    190190         END DO 
    191          CALL lbc_lnk( zfu(:,:,:), 'T', 1. )      ! Lateral boundary conditions  
     191         CALL lbc_lnk( 'traadv_qck', zfu(:,:,:), 'T', 1. )      ! Lateral boundary conditions  
    192192 
    193193         ! 
    … …  
    207207         END DO 
    208208         ! 
    209          CALL lbc_lnk( zwx(:,:,:), 'T', 1. ) ! Lateral boundary conditions 
     209         CALL lbc_lnk( 'traadv_qck', zwx(:,:,:), 'T', 1. ) ! Lateral boundary conditions 
    210210         ! 
    211211         ! Computation of the trend 
    … …  
    265265            END DO 
    266266         END DO 
    267          CALL lbc_lnk_multi( zfc(:,:,:), 'T', 1. , zfd(:,:,:), 'T', 1. )   ! Lateral boundary conditions  
     267         CALL lbc_lnk_multi( 'traadv_qck', zfc(:,:,:), 'T', 1. , zfd(:,:,:), 'T', 1. )   ! Lateral boundary conditions  
    268268 
    269269          
    … …  
    294294 
    295295         !--- Lateral boundary conditions  
    296          CALL lbc_lnk_multi( zfu(:,:,:), 'T', 1. , zfd(:,:,:), 'T', 1., zfc(:,:,:), 'T', 1., zwy(:,:,:), 'T', 1. ) 
     296         CALL lbc_lnk_multi( 'traadv_qck', zfu(:,:,:), 'T', 1. , zfd(:,:,:), 'T', 1., zfc(:,:,:), 'T', 1., zwy(:,:,:), 'T', 1. ) 
    297297 
    298298         !--- QUICKEST scheme 
    … …  
    307307            END DO 
    308308         END DO 
    309          CALL lbc_lnk( zfu(:,:,:), 'T', 1. )    !--- Lateral boundary conditions  
     309         CALL lbc_lnk( 'traadv_qck', zfu(:,:,:), 'T', 1. )    !--- Lateral boundary conditions  
    310310         ! 
    311311         ! Tracer flux on the x-direction 
    … …  
    324324         END DO 
    325325         ! 
    326          CALL lbc_lnk( zwy(:,:,:), 'T', 1. ) ! Lateral boundary conditions 
     326         CALL lbc_lnk( 'traadv_qck', zwy(:,:,:), 'T', 1. ) ! Lateral boundary conditions 
    327327         ! 
    328328         ! Computation of the trend 

  • NEMO/trunk/src/OCE/TRA/traadv_ubs.F90

    r10068 r10425  
    141141            !                                     
    142142         END DO          
    143          CALL lbc_lnk( zltu, 'T', 1. )   ;    CALL lbc_lnk( zltv, 'T', 1. )   ! Lateral boundary cond. (unchanged sgn) 
     143         CALL lbc_lnk( 'traadv_ubs', zltu, 'T', 1. )   ;    CALL lbc_lnk( 'traadv_ubs', zltv, 'T', 1. )   ! Lateral boundary cond. (unchanged sgn) 
    144144         !     
    145145         DO jk = 1, jpkm1        !==  Horizontal advective fluxes  ==!     (UBS) 
    … …  
    226226               END DO 
    227227            END DO 
    228             CALL lbc_lnk( zti, 'T', 1. )      ! Lateral boundary conditions on zti, zsi   (unchanged sign) 
     228            CALL lbc_lnk( 'traadv_ubs', zti, 'T', 1. )      ! Lateral boundary conditions on zti, zsi   (unchanged sign) 
    229229            ! 
    230230            !                          !*  anti-diffusive flux : high order minus low order 

  • NEMO/trunk/src/OCE/TRA/trabbc.F90

    r10068 r10425  
    9292      END DO 
    9393      ! 
    94       CALL lbc_lnk( tsa(:,:,:,jp_tem) , 'T', 1. ) 
     94      CALL lbc_lnk( 'trabbc', tsa(:,:,:,jp_tem) , 'T', 1. ) 
    9595      ! 
    9696      IF( l_trdtra ) THEN        ! Send the trend for diagnostics 

  • NEMO/trunk/src/OCE/TRA/trabbl.F90

    r10068 r10425  
    8484         &      e3u_bbl_0(jpi,jpj) , e3v_bbl_0(jpi,jpj) ,                                    STAT=tra_bbl_alloc ) 
    8585         ! 
    86       IF( lk_mpp            )   CALL mpp_sum ( tra_bbl_alloc ) 
     86      CALL mpp_sum ( 'trabbl', tra_bbl_alloc ) 
    8787      IF( tra_bbl_alloc > 0 )   CALL ctl_warn('tra_bbl_alloc: allocation of arrays failed.') 
    8888   END FUNCTION tra_bbl_alloc 
    … …  
    123123            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=           ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    124124         ! lateral boundary conditions ; just need for outputs 
    125          CALL lbc_lnk_multi( ahu_bbl, 'U', 1. , ahv_bbl, 'V', 1. ) 
     125         CALL lbc_lnk_multi( 'trabbl', 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_multi( utr_bbl, 'U', 1. , vtr_bbl, 'V', 1. ) 
     138         CALL lbc_lnk_multi( 'trabbl', 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 
    … …  
    525525      ! converte into REAL to use lbc_lnk ; impose a min value of 1 as a zero can be set in lbclnk 
    526526      zmbku(:,:) = REAL( mbku_d(:,:), wp )   ;     zmbkv(:,:) = REAL( mbkv_d(:,:), wp )   
    527       CALL lbc_lnk_multi( zmbku,'U',1., zmbkv,'V',1.)  
     527      CALL lbc_lnk_multi( 'trabbl', zmbku,'U',1., zmbkv,'V',1.)  
    528528      mbku_d(:,:) = MAX( INT( zmbku(:,:) ), 1 ) ;  mbkv_d(:,:) = MAX( NINT( zmbkv(:,:) ), 1 ) 
    529529      ! 
    … …  
    548548         END DO 
    549549      END DO 
    550       CALL lbc_lnk_multi( e3u_bbl_0, 'U', 1. , e3v_bbl_0, 'V', 1. )      ! lateral boundary conditions 
     550      CALL lbc_lnk_multi( 'trabbl', e3u_bbl_0, 'U', 1. , e3v_bbl_0, 'V', 1. )      ! lateral boundary conditions 
    551551      ! 
    552552      !                             !* masked diffusive flux coefficients 

  • NEMO/trunk/src/OCE/TRA/tradmp.F90

    r10351 r10425  
    6666      ALLOCATE( resto(jpi,jpj,jpk), STAT= tra_dmp_alloc ) 
    6767      ! 
    68       IF( lk_mpp            )   CALL mpp_sum ( tra_dmp_alloc ) 
     68      CALL mpp_sum ( 'tradmp', tra_dmp_alloc ) 
    6969      IF( tra_dmp_alloc > 0 )   CALL ctl_warn('tra_dmp_alloc: allocation of arrays failed') 
    7070      ! 

  • NEMO/trunk/src/OCE/TRA/traldf_lap_blp.F90

    r10068 r10425  
    213213      END SELECT 
    214214      ! 
    215       CALL lbc_lnk( zlap(:,:,:,:) , 'T', 1. )     ! Lateral boundary conditions (unchanged sign) 
     215      CALL lbc_lnk( 'traldf_lap_blp', zlap(:,:,:,:) , 'T', 1. )     ! Lateral boundary conditions (unchanged sign) 
    216216      !                                               ! Partial top/bottom cell: GRADh( zlap )   
    217217      IF( ln_isfcav .AND. ln_zps ) THEN   ;   CALL zps_hde_isf( kt, kjpt, zlap, zglu, zglv, zgui, zgvi )  ! both top & bottom 

  • NEMO/trunk/src/OCE/TRA/traldf_triad.F90

    r10068 r10425  
    9898      IF( .NOT.ALLOCATED(zdkt3d) )  THEN 
    9999         ALLOCATE( zdkt3d(jpi,jpj,0:1) , STAT=ierr ) 
    100          IF( lk_mpp   )   CALL mpp_sum ( ierr ) 
     100         CALL mpp_sum ( 'traldf_triad', ierr ) 
    101101         IF( ierr > 0 )   CALL ctl_stop('STOP', 'tra_ldf_triad: unable to allocate arrays') 
    102102      ENDIF 

  • NEMO/trunk/src/OCE/TRA/tramle.F90

    r10105 r10425  
    321321               END DO 
    322322            END DO 
    323             CALL lbc_lnk_multi( rfu, 'U', 1. , rfv, 'V', 1. ) 
     323            CALL lbc_lnk_multi( 'tramle', rfu, 'U', 1. , rfv, 'V', 1. ) 
    324324            ! 
    325325         ELSEIF( nn_mle == 1 ) THEN           ! MLE array allocation & initialisation 

  • NEMO/trunk/src/OCE/TRA/tranpc.F90

    r10068 r10425  
    309309         ENDIF 
    310310         ! 
    311          CALL lbc_lnk_multi( tsa(:,:,:,jp_tem), 'T', 1., tsa(:,:,:,jp_sal), 'T', 1. ) 
     311         CALL lbc_lnk_multi( 'tranpc', tsa(:,:,:,jp_tem), 'T', 1., tsa(:,:,:,jp_sal), 'T', 1. ) 
    312312         ! 
    313313         IF( lwp .AND. l_LB_debug ) THEN 

  • NEMO/trunk/src/OCE/TRA/tranxt.F90

    r10095 r10425  
    107107#endif 
    108108      !                                              ! local domain boundaries  (T-point, unchanged sign) 
    109       CALL lbc_lnk_multi( tsa(:,:,:,jp_tem), 'T', 1., tsa(:,:,:,jp_sal), 'T', 1. ) 
     109      CALL lbc_lnk_multi( 'tranxt', tsa(:,:,:,jp_tem), 'T', 1., tsa(:,:,:,jp_sal), 'T', 1. ) 
    110110      ! 
    111111      IF( ln_bdy )   CALL bdy_tra( kt )  ! BDY open boundaries 
    … …  
    163163         ENDIF 
    164164         ! 
    165          CALL lbc_lnk_multi( tsb(:,:,:,jp_tem), 'T', 1., tsb(:,:,:,jp_sal), 'T', 1., & 
     165         CALL lbc_lnk_multi( 'tranxt', tsb(:,:,:,jp_tem), 'T', 1., tsb(:,:,:,jp_sal), 'T', 1., & 
    166166                  &          tsn(:,:,:,jp_tem), 'T', 1., tsn(:,:,:,jp_sal), 'T', 1., & 
    167167                  &          tsa(:,:,:,jp_tem), 'T', 1., tsa(:,:,:,jp_sal), 'T', 1.  ) 

  • NEMO/trunk/src/OCE/TRA/traqsr.F90

    r10068 r10425  
    275275         END DO 
    276276      END DO 
    277       CALL lbc_lnk( fraqsr_1lev(:,:), 'T', 1._wp ) 
     277      CALL lbc_lnk( 'traqsr', fraqsr_1lev(:,:), 'T', 1._wp ) 
    278278      ! 
    279279      IF( iom_use('qsr3d') ) THEN      ! output the shortwave Radiation distribution 

  • NEMO/trunk/src/OCE/TRA/trazdf.F90

    r10364 r10425  
    9292         END DO 
    9393!!gm this should be moved in trdtra.F90 and done on all trends 
    94          CALL lbc_lnk_multi( ztrdt, 'T', 1. , ztrds, 'T', 1. ) 
     94         CALL lbc_lnk_multi( 'trazdf', ztrdt, 'T', 1. , ztrds, 'T', 1. ) 
    9595!!gm 
    9696         CALL trd_tra( kt, 'TRA', jp_tem, jptra_zdf, ztrdt ) 

  • NEMO/trunk/src/OCE/TRA/zpshde.F90

    r10068 r10425  
    144144            END DO 
    145145         END DO 
    146          CALL lbc_lnk_multi( pgtu(:,:,jn), 'U', -1. , pgtv(:,:,jn), 'V', -1. )   ! Lateral boundary cond. 
    147          ! 
    148146      END DO 
     147      ! 
     148      CALL lbc_lnk_multi( 'zpshde', pgtu(:,:,:), 'U', -1. , pgtv(:,:,:), 'V', -1. )   ! Lateral boundary cond. 
    149149      !                 
    150150      IF( PRESENT( prd ) ) THEN    !==  horizontal derivative of density anomalies (rd)  ==!    (optional part) 
    … …  
    183183            END DO 
    184184         END DO 
    185          CALL lbc_lnk_multi( pgru , 'U', -1. , pgrv , 'V', -1. )   ! Lateral boundary conditions 
     185         CALL lbc_lnk_multi( 'zpshde', pgru , 'U', -1. , pgrv , 'V', -1. )   ! Lateral boundary conditions 
    186186         ! 
    187187      END IF 
    … …  
    305305            END DO 
    306306         END DO 
    307          CALL lbc_lnk_multi( pgtu(:,:,jn), 'U', -1. , pgtv(:,:,jn), 'V', -1. )   ! Lateral boundary cond. 
    308          ! 
    309307      END DO 
     308      ! 
     309      CALL lbc_lnk_multi( 'zpshde', pgtu(:,:,:), 'U', -1. , pgtv(:,:,:), 'V', -1. )   ! Lateral boundary cond. 
    310310 
    311311      ! horizontal derivative of density anomalies (rd) 
    … …  
    353353         END DO 
    354354 
    355          CALL lbc_lnk_multi( pgru , 'U', -1. , pgrv , 'V', -1. )   ! Lateral boundary conditions 
     355         CALL lbc_lnk_multi( 'zpshde', pgru , 'U', -1. , pgrv , 'V', -1. )   ! Lateral boundary conditions 
    356356         ! 
    357357      END IF 
    … …  
    406406         ! 
    407407      END DO 
    408       CALL lbc_lnk_multi( pgtui(:,:,:), 'U', -1. , pgtvi(:,:,:), 'V', -1. )   ! Lateral boundary cond. 
     408      CALL lbc_lnk_multi( 'zpshde', 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_multi( pgrui, 'U', -1. , pgrvi, 'V', -1. )   ! Lateral boundary conditions 
     451         CALL lbc_lnk_multi( 'zpshde', pgrui, 'U', -1. , pgrvi, 'V', -1. )   ! Lateral boundary conditions 
    452452         ! 
    453453      END IF   
Note: See TracChangeset for help on using the changeset viewer.