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 911 for trunk/NEMO/OPA_SRC/DYN/wzvmod.F90 – NEMO

Ignore:
Timestamp:
2008-04-28T11:31:32+02:00 (16 years ago)
Author:
ctlod
Message:

Implementation of the BDY package, see ticket: #126

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/DYN/wzvmod.F90

    r888 r911  
    77   !!             7.0  !  96-01  (G. Madec)  Statement function for e3 
    88   !!             8.5  !  02-07  (G. Madec)  Free form, F90 
     9   !!              "   !  07-07  (D. Storkey) Zero zhdiv at open boundary (BDY)  
    910   !!---------------------------------------------------------------------- 
    1011   !!   wzv        : Compute the vertical velocity 
     
    1819   USE prtctl          ! Print control 
    1920   USE phycst 
     21   USE bdy_oce         ! unstructured open boundaries 
    2022   USE lbclnk          ! ocean lateral boundary condition (or mpp link) 
    2123 
     
    5456 
    5557      !! * Local declarations 
     58      INTEGER  ::     jgrd, jb           ! temporary scalars 
    5659      INTEGER  ::           jk           ! dummy loop indices 
    5760      !! Variable volume 
     
    8891         ! Horizontal divergence of barotropic transports 
    8992         !-------------------------------------------------- 
     93         zhdiv(:,:) = 0.e0 
    9094         DO jj = 2, jpjm1 
    9195            DO ji = 2, jpim1   ! vector opt. 
     
    105109         IF( lp_obc_north )   zhdiv(nin0  :nin1  ,njn0p1:njn1p1) = 0.e0    ! north 
    106110         IF( lp_obc_south )   zhdiv(nis0  :nis1  ,njs0  :njs1)   = 0.e0    ! south 
     111#endif 
     112 
     113#if defined key_bdy || defined key_bdy_tides 
     114         jgrd=1 !: tracer grid. 
     115         DO jb = 1, nblenrim(jgrd) 
     116           ji = nbi(jb,jgrd) 
     117           jj = nbj(jb,jgrd) 
     118           zhdiv(ji, jj) = 0.e0 
     119           zhdiv(ji, jj) = 0.e0 
     120           zhdiv(ji, jj) = 0.e0 
     121         END DO 
    107122#endif 
    108123 
Note: See TracChangeset for help on using the changeset viewer.