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.
bdydyn.F90 in NEMO/trunk/src/OCE/BDY – NEMO

source: NEMO/trunk/src/OCE/BDY/bdydyn.F90 @ 14325

Last change on this file since 14325 was 13237, checked in by smasson, 4 years ago

trunk: Mid-year merge, merge back KERNEL-06_techene_e3

  • Property svn:keywords set to Id
File size: 5.7 KB
RevLine 
[911]1MODULE bdydyn
[1125]2   !!======================================================================
[911]3   !!                       ***  MODULE  bdydyn  ***
[3294]4   !! Unstructured Open Boundary Cond. :   Apply boundary conditions to velocities
[1125]5   !!======================================================================
6   !! History :  1.0  !  2005-02  (J. Chanut, A. Sellar)  Original code
7   !!             -   !  2007-07  (D. Storkey) Move Flather implementation to separate routine.
8   !!            3.0  !  2008-04  (NEMO team)  add in the reference version
[1502]9   !!            3.2  !  2008-04  (R. Benshila) consider velocity instead of transport
[2528]10   !!            3.3  !  2010-09  (E.O'Dea) modifications for Shelf configurations
11   !!            3.3  !  2010-09  (D.Storkey) add ice boundary conditions
[3294]12   !!            3.4  !  2011     (D. Storkey) rewrite in preparation for OBC-BDY merge
[1125]13   !!----------------------------------------------------------------------
[3294]14   !!   bdy_dyn        : split velocities into barotropic and baroclinic parts
15   !!                    and call bdy_dyn2d and bdy_dyn3d to apply boundary
16   !!                    conditions
[1125]17   !!----------------------------------------------------------------------
[911]18   USE oce             ! ocean dynamics and tracers
19   USE dom_oce         ! ocean space and time domain
20   USE bdy_oce         ! ocean open boundary conditions
[3294]21   USE bdydyn2d        ! open boundary conditions for barotropic solution
22   USE bdydyn3d        ! open boundary conditions for baroclinic velocities
[911]23   USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
[1125]24   USE in_out_manager  !
[4689]25   USE domvvl          ! variable volume
[911]26
27   IMPLICIT NONE
28   PRIVATE
29
[5930]30   PUBLIC   bdy_dyn    ! routine called in dyn_nxt
[911]31
[13237]32   !! * Substitutions
33#  include "domzgr_substitute.h90"
[1125]34   !!----------------------------------------------------------------------
[9598]35   !! NEMO/OCE 4.0 , NEMO Consortium (2018)
[1146]36   !! $Id$
[10068]37   !! Software governed by the CeCILL license (see ./LICENSE)
[1125]38   !!----------------------------------------------------------------------
[911]39CONTAINS
40
[12377]41   SUBROUTINE bdy_dyn( kt, Kbb, puu, pvv, Kaa, dyn3d_only )
[1125]42      !!----------------------------------------------------------------------
[3294]43      !!                  ***  SUBROUTINE bdy_dyn  ***
[1125]44      !!
[3294]45      !! ** Purpose : - Wrapper routine for bdy_dyn2d and bdy_dyn3d.
[911]46      !!
[1125]47      !!----------------------------------------------------------------------
[12377]48      INTEGER                             , INTENT(in)    ::   kt           ! Main time step counter
49      INTEGER                             , INTENT(in)    ::   Kbb, Kaa     ! Ocean time level indices
50      REAL(wp), DIMENSION(jpi,jpj,jpk,jpt), INTENT(inout) ::   puu, pvv     ! Ocean velocities (to be updated at open boundaries)
51      LOGICAL, OPTIONAL                   , INTENT(in)    ::   dyn3d_only   ! T => only update baroclinic velocities
[6140]52      !
53      INTEGER ::   jk, ii, ij, ib_bdy, ib, igrd     ! Loop counter
54      LOGICAL ::   ll_dyn2d, ll_dyn3d, ll_orlanski
[12377]55      REAL(wp), DIMENSION(jpi,jpj) :: zua2d, zva2d     ! after barotropic velocities
[6140]56      !!----------------------------------------------------------------------
57      !
[3294]58      ll_dyn2d = .true.
59      ll_dyn3d = .true.
[6140]60      !
[3294]61      IF( PRESENT(dyn3d_only) ) THEN
[6140]62         IF( dyn3d_only )   ll_dyn2d = .false.
[3294]63      ENDIF
[6140]64      !
[4292]65      ll_orlanski = .false.
66      DO ib_bdy = 1, nb_bdy
[6140]67         IF ( cn_dyn2d(ib_bdy) == 'orlanski' .OR. cn_dyn2d(ib_bdy) == 'orlanski_npo' &
68     &   .OR. cn_dyn3d(ib_bdy) == 'orlanski' .OR. cn_dyn3d(ib_bdy) == 'orlanski_npo')   ll_orlanski = .true.
69      END DO
[4292]70
[3294]71      !-------------------------------------------------------
72      ! Split velocities into barotropic and baroclinic parts
73      !-------------------------------------------------------
[911]74
[6140]75      !                          ! "After" velocities:
[12377]76      zua2d(:,:) = 0._wp
77      zva2d(:,:) = 0._wp     
[4370]78      DO jk = 1, jpkm1
[12377]79         zua2d(:,:) = zua2d(:,:) + e3u(:,:,jk,Kaa) * puu(:,:,jk,Kaa) * umask(:,:,jk)
80         zva2d(:,:) = zva2d(:,:) + e3v(:,:,jk,Kaa) * pvv(:,:,jk,Kaa) * vmask(:,:,jk)
[4370]81      END DO
[12377]82      zua2d(:,:) = zua2d(:,:) * r1_hu(:,:,Kaa)
83      zva2d(:,:) = zva2d(:,:) * r1_hv(:,:,Kaa)
[4292]84
[3294]85      DO jk = 1 , jpkm1
[12377]86         puu(:,:,jk,Kaa) = ( puu(:,:,jk,Kaa) - zua2d(:,:) ) * umask(:,:,jk)
87         pvv(:,:,jk,Kaa) = ( pvv(:,:,jk,Kaa) - zva2d(:,:) ) * vmask(:,:,jk)
[3294]88      END DO
89
[4292]90
[6140]91      IF( ll_orlanski ) THEN     ! "Before" velocities (Orlanski condition only)
[4292]92         DO jk = 1 , jpkm1
[12377]93            puu(:,:,jk,Kbb) = ( puu(:,:,jk,Kbb) - uu_b(:,:,Kbb) ) * umask(:,:,jk)
94            pvv(:,:,jk,Kbb) = ( pvv(:,:,jk,Kbb) - vv_b(:,:,Kbb) ) * vmask(:,:,jk)
[4292]95         END DO
[6140]96      ENDIF
[4292]97
[3294]98      !-------------------------------------------------------
99      ! Apply boundary conditions to barotropic and baroclinic
100      ! parts separately
101      !-------------------------------------------------------
102
[12377]103      IF( ll_dyn2d )   CALL bdy_dyn2d( kt, zua2d, zva2d, uu_b(:,:,Kbb), vv_b(:,:,Kbb), r1_hu(:,:,Kaa), r1_hv(:,:,Kaa), ssh(:,:,Kaa) )
[3294]104
[12377]105      IF( ll_dyn3d )   CALL bdy_dyn3d( kt, Kbb, puu, pvv, Kaa )
[3294]106
107      !-------------------------------------------------------
108      ! Recombine velocities
109      !-------------------------------------------------------
[6140]110      !
[3294]111      DO jk = 1 , jpkm1
[12377]112         puu(:,:,jk,Kaa) = ( puu(:,:,jk,Kaa) + zua2d(:,:) ) * umask(:,:,jk)
113         pvv(:,:,jk,Kaa) = ( pvv(:,:,jk,Kaa) + zva2d(:,:) ) * vmask(:,:,jk)
[3294]114      END DO
[6140]115      !
[4292]116      IF ( ll_orlanski ) THEN
117         DO jk = 1 , jpkm1
[12377]118            puu(:,:,jk,Kbb) = ( puu(:,:,jk,Kbb) + uu_b(:,:,Kbb) ) * umask(:,:,jk)
119            pvv(:,:,jk,Kbb) = ( pvv(:,:,jk,Kbb) + vv_b(:,:,Kbb) ) * vmask(:,:,jk)
[4292]120         END DO
121      END IF
[6140]122      !
[3294]123   END SUBROUTINE bdy_dyn
124
[1125]125   !!======================================================================
[911]126END MODULE bdydyn
Note: See TracBrowser for help on using the repository browser.