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.
dynbfr.F90 in branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DYN – NEMO

source: branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DYN/dynbfr.F90 @ 9124

Last change on this file since 9124 was 9124, checked in by gm, 6 years ago

dev_merge_2017: ln_timing instead of nn_timing + restricted timing to nemo_init and routine called by step in OPA_SRC

  • Property svn:keywords set to Id
File size: 5.1 KB
Line 
1MODULE dynbfr
2   !!==============================================================================
3   !!                 ***  MODULE  dynbfr  ***
4   !! Ocean dynamics :  bottom friction component of the momentum mixing trend
5   !!==============================================================================
6   !! History :  3.2  ! 2008-11  (A. C. Coward)  Original code
7   !!            3.4  ! 2011-09  (H. Liu) Make it consistent with semi-implicit Bottom friction (ln_drgimp =T)
8   !!            4.0  ! 2017-05  (G. Madec)  drag coef. defined at t-point (zdfdrg.F90)
9   !!----------------------------------------------------------------------
10
11   !!----------------------------------------------------------------------
12   !!   dyn_bfr       : Update the momentum trend with the bottom friction contribution
13   !!----------------------------------------------------------------------
14   USE oce            ! ocean dynamics and tracers variables
15   USE dom_oce        ! ocean space and time domain variables
16   USE zdf_oce        ! vertical physics: variables
17   USE zdfdrg         ! vertical physics: top/bottom drag coef.
18   USE trd_oce        ! trends: ocean variables
19   USE trddyn         ! trend manager: dynamics
20   !
21   USE in_out_manager ! I/O manager
22   USE prtctl         ! Print control
23
24   IMPLICIT NONE
25   PRIVATE
26
27   PUBLIC   dyn_bfr   !  routine called by step.F90
28
29   !! * Substitutions
30#  include "vectopt_loop_substitute.h90"
31   !!----------------------------------------------------------------------
32   !! NEMO/OPA 4.0 , NEMO Consortium (2017)
33   !! $Id$
34   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
35   !!----------------------------------------------------------------------
36CONTAINS
37   
38   SUBROUTINE dyn_bfr( kt )
39      !!----------------------------------------------------------------------
40      !!                  ***  ROUTINE dyn_bfr  ***
41      !!
42      !! ** Purpose :   compute the bottom friction ocean dynamics physics.
43      !!
44      !!              only for explicit bottom friction form
45      !!              implicit bfr is implemented in dynzdf_imp
46      !!
47      !! ** Action  :   (ua,va)   momentum trend increased by bottom friction trend
48      !!---------------------------------------------------------------------
49      INTEGER, INTENT(in) ::   kt   ! ocean time-step index
50      !!
51      INTEGER  ::   ji, jj       ! dummy loop indexes
52      INTEGER  ::   ikbu, ikbv   ! local integers
53      REAL(wp) ::   zm1_2dt      ! local scalar
54      REAL(wp) ::   zCdu, zCdv   !   -      -
55      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::   ztrdu, ztrdv
56      !!---------------------------------------------------------------------
57      !
58!!gm bug : time step is only rdt (not 2 rdt if euler start !)
59         zm1_2dt = - 1._wp / ( 2._wp * rdt )
60
61      IF( l_trddyn ) THEN      ! trends: store the input trends
62         ALLOCATE( ztrdu(jpi,jpj,jpk) , ztrdv(jpi,jpj,jpk) )
63         ztrdu(:,:,:) = ua(:,:,:)
64            ztrdv(:,:,:) = va(:,:,:)
65      ENDIF
66
67
68      DO jj = 2, jpjm1
69         DO ji = 2, jpim1
70            ikbu = mbku(ji,jj)          ! deepest wet ocean u- & v-levels
71            ikbv = mbkv(ji,jj)
72            !
73            ! Apply stability criteria on absolute value  : abs(bfr/e3) < 1/(2dt) => bfr/e3 > -1/(2dt)
74            zCdu = 0.5*( rCdU_bot(ji+1,jj)+rCdU_bot(ji,jj) ) / e3u_n(ji,jj,ikbu)
75            zCdv = 0.5*( rCdU_bot(ji,jj+1)+rCdU_bot(ji,jj) ) / e3v_n(ji,jj,ikbv)
76            !
77            ua(ji,jj,ikbu) = ua(ji,jj,ikbu) + MAX(  zCdu , zm1_2dt  ) * ub(ji,jj,ikbu)
78            va(ji,jj,ikbv) = va(ji,jj,ikbv) + MAX(  zCdv , zm1_2dt  ) * vb(ji,jj,ikbv)
79         END DO
80      END DO
81      !
82      IF( ln_isfcav ) THEN        ! ocean cavities
83         DO jj = 2, jpjm1
84            DO ji = 2, jpim1
85               ikbu = miku(ji,jj)          ! first wet ocean u- & v-levels
86               ikbv = mikv(ji,jj)
87               !
88               ! Apply stability criteria on absolute value  : abs(bfr/e3) < 1/(2dt) => bfr/e3 > -1/(2dt)
89               zCdu = 0.5*( rCdU_top(ji+1,jj)+rCdU_top(ji,jj) ) / e3u_n(ji,jj,ikbu)    ! NB: Cdtop masked
90               zCdv = 0.5*( rCdU_top(ji,jj+1)+rCdU_top(ji,jj) ) / e3v_n(ji,jj,ikbv)
91               !
92               ua(ji,jj,ikbu) = ua(ji,jj,ikbu) + MAX(  zCdu , zm1_2dt  ) * ub(ji,jj,ikbu)
93               va(ji,jj,ikbv) = va(ji,jj,ikbv) + MAX(  zCdv , zm1_2dt  ) * vb(ji,jj,ikbv)
94           END DO
95         END DO
96      ENDIF
97      !
98      IF( l_trddyn ) THEN      ! trends: send trends to trddyn for further diagnostics
99         ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:)
100         ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:)
101         CALL trd_dyn( ztrdu(:,:,:), ztrdv(:,:,:), jpdyn_bfr, kt )
102         DEALLOCATE( ztrdu, ztrdv )
103      ENDIF
104      !                                          ! print mean trends (used for debugging)
105      IF(ln_ctl)   CALL prt_ctl( tab3d_1=ua, clinfo1=' bfr  - Ua: ', mask1=umask,               &
106         &                       tab3d_2=va, clinfo2=       ' Va: ', mask2=vmask, clinfo3='dyn' )
107      !
108   END SUBROUTINE dyn_bfr
109
110   !!==============================================================================
111END MODULE dynbfr
Note: See TracBrowser for help on using the repository browser.