source: trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynbfr.F90 @ 3294

Last change on this file since 3294 was 3294, checked in by rblod, 9 years ago

Merge of 3.4beta into the trunk

  • Property svn:keywords set to Id
File size: 4.6 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
8   !!                            Bottom friction (ln_bfrimp = .true.)
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         ! ocean vertical physics variables
17   USE zdfbfr          ! ocean bottom friction variables
18   USE trdmod          ! ocean active dynamics and tracers trends
19   USE trdmod_oce      ! ocean variables trends
20   USE in_out_manager  ! I/O manager
21   USE prtctl          ! Print control
22   USE timing          ! Timing
23   USE wrk_nemo        ! Memory Allocation
24
25   IMPLICIT NONE
26   PRIVATE
27
28   PUBLIC   dyn_bfr    !  routine called by step.F90
29
30   !! * Substitutions
31#  include "domzgr_substitute.h90"
32#  include "zdfddm_substitute.h90"
33#  include "vectopt_loop_substitute.h90"
34   !!----------------------------------------------------------------------
35   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
36   !! $Id$
37   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
38   !!----------------------------------------------------------------------
39CONTAINS
40   
41   SUBROUTINE dyn_bfr( kt )
42      !!----------------------------------------------------------------------
43      !!                  ***  ROUTINE dyn_bfr  ***
44      !!
45      !! ** Purpose :   compute the bottom friction ocean dynamics physics.
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), POINTER, DIMENSION(:,:,:) ::  ztrdu, ztrdv
55      !!---------------------------------------------------------------------
56      !
57      IF( nn_timing == 1 )  CALL timing_start('dyn_bfr')
58      !
59      IF( .NOT.ln_bfrimp) THEN     ! only for explicit bottom friction form
60                                    ! implicit bfr is implemented in dynzdf_imp
61
62        zm1_2dt = - 1._wp / ( 2._wp * rdt )
63
64        IF( l_trddyn )   THEN                      ! temporary save of ua and va trends
65           CALL wrk_alloc( jpi,jpj,jpk, ztrdu, ztrdv )
66           ztrdu(:,:,:) = ua(:,:,:)
67           ztrdv(:,:,:) = va(:,:,:)
68        ENDIF
69
70
71# if defined key_vectopt_loop
72        DO jj = 1, 1
73           DO ji = jpi+2, jpij-jpi-1   ! vector opt. (forced unrolling)
74# else
75        DO jj = 2, jpjm1
76           DO ji = 2, jpim1
77# endif
78              ikbu = mbku(ji,jj)          ! deepest ocean u- & v-levels
79              ikbv = mbkv(ji,jj)
80              !
81              ! Apply stability criteria on absolute value  : abs(bfr/e3) < 1/(2dt) => bfr/e3 > -1/(2dt)
82              ua(ji,jj,ikbu) = ua(ji,jj,ikbu) + MAX(  bfrua(ji,jj) / fse3u(ji,jj,ikbu) , zm1_2dt  ) * ub(ji,jj,ikbu)
83              va(ji,jj,ikbv) = va(ji,jj,ikbv) + MAX(  bfrva(ji,jj) / fse3v(ji,jj,ikbv) , zm1_2dt  ) * vb(ji,jj,ikbv)
84           END DO
85        END DO
86
87        !
88        IF( l_trddyn )   THEN                      ! save the vertical diffusive trends for further diagnostics
89           ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:)
90           ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:)
91           CALL trd_mod( ztrdu(:,:,:), ztrdv(:,:,:), jpdyn_trd_bfr, 'DYN', kt )
92           CALL wrk_dealloc( jpi,jpj,jpk, ztrdu, ztrdv )
93        ENDIF
94        !                                          ! print mean trends (used for debugging)
95        IF(ln_ctl)   CALL prt_ctl( tab3d_1=ua, clinfo1=' bfr  - Ua: ', mask1=umask,               &
96           &                       tab3d_2=va, clinfo2=       ' Va: ', mask2=vmask, clinfo3='dyn' )
97        !
98      ENDIF     ! end explicit bottom friction
99      !
100      IF( nn_timing == 1 )  CALL timing_stop('dyn_bfr')
101      !
102   END SUBROUTINE dyn_bfr
103
104   !!==============================================================================
105END MODULE dynbfr
Note: See TracBrowser for help on using the repository browser.