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 3583 for branches/2012/dev_MERCATOR_2012_rev3555/NEMOGCM/NEMO/OPA_SRC/BDY/bdydyn3d.F90 – NEMO

Ignore:
Timestamp:
2012-11-16T17:18:17+01:00 (10 years ago)
Author:
cbricaud
Message:

add modification from dev_r3327_MERCATOR1_BDY branch in dev_MERCATOR_2012_rev3555 branch

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_MERCATOR_2012_rev3555/NEMOGCM/NEMO/OPA_SRC/BDY/bdydyn3d.F90

    r3294 r3583  
    1414   !!---------------------------------------------------------------------- 
    1515   USE timing          ! Timing 
     16   USE wrk_nemo        ! Memory Allocation 
    1617   USE oce             ! ocean dynamics and tracers  
    1718   USE dom_oce         ! ocean space and time domain 
     
    1920   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    2021   USE in_out_manager  ! 
     22   Use phycst 
    2123 
    2224   IMPLICIT NONE 
     
    2426 
    2527   PUBLIC   bdy_dyn3d     ! routine called by bdy_dyn 
    26  
     28   PUBLIC   bdy_dyn3d_dmp ! routine called by step 
     29 
     30   !! * Substitutions 
     31#  include "domzgr_substitute.h90" 
    2732   !!---------------------------------------------------------------------- 
    2833   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
     
    5560         CASE(jp_frs) 
    5661            CALL bdy_dyn3d_frs( idx_bdy(ib_bdy), dta_bdy(ib_bdy), kt ) 
     62         CASE(2) 
     63            CALL bdy_dyn3d_spe( idx_bdy(ib_bdy), dta_bdy(ib_bdy), kt ) 
     64         CASE(3) 
     65            CALL bdy_dyn3d_zro( idx_bdy(ib_bdy), dta_bdy(ib_bdy), kt ) 
    5766         CASE DEFAULT 
    5867            CALL ctl_stop( 'bdy_dyn3d : unrecognised option for open boundaries for baroclinic velocities' ) 
     
    6271   END SUBROUTINE bdy_dyn3d 
    6372 
     73   SUBROUTINE bdy_dyn3d_spe( idx, dta, kt ) 
     74      !!---------------------------------------------------------------------- 
     75      !!                  ***  SUBROUTINE bdy_dyn3d_spe  *** 
     76      !! 
     77      !! ** Purpose : - Apply a specified value for baroclinic velocities 
     78      !!                at open boundaries. 
     79      !! 
     80      !!---------------------------------------------------------------------- 
     81      INTEGER                     ::   kt 
     82      TYPE(OBC_INDEX), INTENT(in) ::   idx  ! OBC indices 
     83      TYPE(OBC_DATA),  INTENT(in) ::   dta  ! OBC external data 
     84      !! 
     85      INTEGER  ::   jb, jk         ! dummy loop indices 
     86      INTEGER  ::   ii, ij, igrd   ! local integers 
     87      REAL(wp) ::   zwgt           ! boundary weight 
     88      !!---------------------------------------------------------------------- 
     89      ! 
     90      IF( nn_timing == 1 ) CALL timing_start('bdy_dyn3d_spe') 
     91      ! 
     92      igrd = 2                      ! Relaxation of zonal velocity 
     93      DO jb = 1, idx%nblenrim(igrd) 
     94         DO jk = 1, jpkm1 
     95            ii   = idx%nbi(jb,igrd) 
     96            ij   = idx%nbj(jb,igrd) 
     97            ua(ii,ij,jk) = dta%u3d(jb,jk) * umask(ii,ij,jk) 
     98         END DO 
     99      END DO 
     100      ! 
     101      igrd = 3                      ! Relaxation of meridional velocity 
     102      DO jb = 1, idx%nblenrim(igrd) 
     103         DO jk = 1, jpkm1 
     104            ii   = idx%nbi(jb,igrd) 
     105            ij   = idx%nbj(jb,igrd) 
     106            va(ii,ij,jk) = dta%v3d(jb,jk) * vmask(ii,ij,jk) 
     107         END DO 
     108      END DO 
     109      CALL lbc_lnk( ua, 'U', -1. )   ;   CALL lbc_lnk( va, 'V', -1. )   ! Boundary points should be updated 
     110      ! 
     111      IF( kt .eq. nit000 ) CLOSE( unit = 102 ) 
     112 
     113      IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn3d_spe') 
     114 
     115   END SUBROUTINE bdy_dyn3d_spe 
     116 
     117   SUBROUTINE bdy_dyn3d_zro( idx, dta, kt ) 
     118      !!---------------------------------------------------------------------- 
     119      !!                  ***  SUBROUTINE bdy_dyn3d_zro  *** 
     120      !! 
     121      !! ** Purpose : - baroclinic velocities = 0. at open boundaries. 
     122      !! 
     123      !!---------------------------------------------------------------------- 
     124      INTEGER                     ::   kt 
     125      TYPE(OBC_INDEX), INTENT(in) ::   idx  ! OBC indices 
     126      TYPE(OBC_DATA),  INTENT(in) ::   dta  ! OBC external data 
     127      !! 
     128      INTEGER  ::   ib, ik         ! dummy loop indices 
     129      INTEGER  ::   ii, ij, igrd, zcoef   ! local integers 
     130      REAL(wp) ::   zwgt           ! boundary weight 
     131      !!---------------------------------------------------------------------- 
     132      ! 
     133      IF( nn_timing == 1 ) CALL timing_start('bdy_dyn3d_zro') 
     134      ! 
     135      igrd = 2                       ! Everything is at T-points here 
     136      DO ib = 1, idx%nblenrim(igrd) 
     137         ii = idx%nbi(ib,igrd) 
     138         ij = idx%nbj(ib,igrd) 
     139         DO ik = 1, jpkm1 
     140            ua(ii,ij,ik) = 0._wp 
     141         END DO 
     142      END DO 
     143 
     144      igrd = 3                       ! Everything is at T-points here 
     145      DO ib = 1, idx%nblenrim(igrd) 
     146         ii = idx%nbi(ib,igrd) 
     147         ij = idx%nbj(ib,igrd) 
     148         DO ik = 1, jpkm1 
     149            va(ii,ij,ik) = 0._wp 
     150         END DO 
     151      END DO 
     152      ! 
     153      CALL lbc_lnk( ua, 'U', -1. )   ;   CALL lbc_lnk( va, 'V', -1. )   ! Boundary points should be updated 
     154      ! 
     155      IF( kt .eq. nit000 ) CLOSE( unit = 102 ) 
     156 
     157      IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn3d_zro') 
     158 
     159   END SUBROUTINE bdy_dyn3d_zro 
     160 
    64161   SUBROUTINE bdy_dyn3d_frs( idx, dta, kt ) 
    65162      !!---------------------------------------------------------------------- 
     
    111208   END SUBROUTINE bdy_dyn3d_frs 
    112209 
     210   SUBROUTINE bdy_dyn3d_dmp( kt ) 
     211      !!---------------------------------------------------------------------- 
     212      !!                  ***  SUBROUTINE bdy_dyn3d_dmp  *** 
     213      !! 
     214      !! ** Purpose : Apply damping for baroclinic velocities at open boundaries. 
     215      !! 
     216      !!---------------------------------------------------------------------- 
     217      INTEGER                     ::   kt 
     218      !! 
     219      INTEGER  ::   jb, jk         ! dummy loop indices 
     220      INTEGER  ::   ii, ij, igrd   ! local integers 
     221      REAL(wp) ::   zwgt           ! boundary weight 
     222      INTEGER  ::  ib_bdy          ! loop index 
     223      !!---------------------------------------------------------------------- 
     224      ! 
     225      IF( nn_timing == 1 ) CALL timing_start('bdy_dyn3d_dmp') 
     226      ! 
     227      !------------------------------------------------------- 
     228      ! Remove barotropic part from before velocity 
     229      !------------------------------------------------------- 
     230      CALL wrk_alloc(jpi,jpj,pu2d,pv2d)  
     231 
     232      pu2d(:,:) = 0.e0 
     233      pv2d(:,:) = 0.e0 
     234 
     235      DO jk = 1, jpkm1 
     236#if defined key_vvl 
     237         pu2d(:,:) = pu2d(:,:) + fse3u_b(:,:,jk)* ub(:,:,jk)   *umask(:,:,jk)  
     238         pv2d(:,:) = pv2d(:,:) + fse3v_b(:,:,jk)* vb(:,:,jk)   *vmask(:,:,jk) 
     239#else 
     240         pu2d(:,:) = pu2d(:,:) + fse3u_0(:,:,jk) * ub(:,:,jk)  * umask(:,:,jk) 
     241         pv2d(:,:) = pv2d(:,:) + fse3v_0(:,:,jk) * vb(:,:,jk)  * vmask(:,:,jk) 
     242#endif 
     243      END DO 
     244 
     245      IF( lk_vvl ) THEN 
     246         pu2d(:,:) = pu2d(:,:) * umask(:,:,1) / ( hu_0(:,:) + sshu_b(:,:) + 1._wp - umask(:,:,1) ) 
     247         pv2d(:,:) = pv2d(:,:) * vmask(:,:,1) / ( hv_0(:,:) + sshv_b(:,:) + 1._wp - vmask(:,:,1) ) 
     248      ELSE 
     249         pu2d(:,:) = pv2d(:,:) * hur(:,:) 
     250         pv2d(:,:) = pu2d(:,:) * hvr(:,:) 
     251      ENDIF 
     252 
     253      DO ib_bdy=1, nb_bdy 
     254         IF ( ln_dyn3d_dmp(ib_bdy).and.nn_dyn3d(ib_bdy).gt.0 ) THEN 
     255            igrd = 2                      ! Relaxation of zonal velocity 
     256            DO jb = 1, idx_bdy(ib_bdy)%nblen(igrd) 
     257               ii   = idx_bdy(ib_bdy)%nbi(jb,igrd) 
     258               ij   = idx_bdy(ib_bdy)%nbj(jb,igrd) 
     259               zwgt = idx_bdy(ib_bdy)%nbd(jb,igrd) 
     260               DO jk = 1, jpkm1 
     261                  ua(ii,ij,jk) = ( ua(ii,ij,jk) + zwgt * ( dta_bdy(ib_bdy)%u3d(jb,jk) - ub(ii,ij,jk) + pu2d(ii,ij)) ) * umask(ii,ij,jk) 
     262               END DO 
     263            END DO 
     264            ! 
     265            igrd = 3                      ! Relaxation of meridional velocity 
     266            DO jb = 1, idx_bdy(ib_bdy)%nblen(igrd) 
     267               ii   = idx_bdy(ib_bdy)%nbi(jb,igrd) 
     268               ij   = idx_bdy(ib_bdy)%nbj(jb,igrd) 
     269               zwgt = idx_bdy(ib_bdy)%nbd(jb,igrd) 
     270               DO jk = 1, jpkm1 
     271                  va(ii,ij,jk) = ( va(ii,ij,jk) + zwgt * ( dta_bdy(ib_bdy)%v3d(jb,jk) - vb(ii,ij,jk) + pv2d(ii,ij)) ) * vmask(ii,ij,jk) 
     272               END DO 
     273            END DO 
     274         ENDIF 
     275      ENDDO 
     276      ! 
     277      CALL wrk_dealloc(jpi,jpj,pu2d,pv2d)  
     278      ! 
     279      CALL lbc_lnk( ua, 'U', -1. )   ;   CALL lbc_lnk( va, 'V', -1. )   ! Boundary points should be updated 
     280      ! 
     281      IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn3d_dmp') 
     282 
     283   END SUBROUTINE bdy_dyn3d_dmp 
    113284 
    114285#else 
Note: See TracChangeset for help on using the changeset viewer.