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.
dynldf.F90 in NEMO/branches/UKMO/dev_r10448_bdyvol/src/OCE/DYN – NEMO

source: NEMO/branches/UKMO/dev_r10448_bdyvol/src/OCE/DYN/dynldf.F90 @ 10455

Last change on this file since 10455 was 10068, checked in by nicolasmartin, 6 years ago

First part of modifications to have a common default header : fix typos and SVN keywords properties

  • Property svn:keywords set to Id
File size: 5.4 KB
Line 
1MODULE dynldf
2   !!======================================================================
3   !!                       ***  MODULE  dynldf  ***
4   !! Ocean physics:  lateral diffusivity trends
5   !!=====================================================================
6   !! History :  2.0  ! 2005-11  (G. Madec)  Original code (new step architecture)
7   !!            3.7  ! 2014-01  (F. Lemarie, G. Madec)  restructuration/simplification of ahm specification,
8   !!                 !                                  add velocity dependent coefficient and optional read in file
9   !!----------------------------------------------------------------------
10
11   !!----------------------------------------------------------------------
12   !!   dyn_ldf      : update the dynamics trend with the lateral diffusion
13   !!   dyn_ldf_init : initialization, namelist read, and parameters control
14   !!----------------------------------------------------------------------
15   USE oce            ! ocean dynamics and tracers
16   USE dom_oce        ! ocean space and time domain
17   USE phycst         ! physical constants
18   USE ldfdyn         ! lateral diffusion: eddy viscosity coef.
19   USE dynldf_lap_blp ! lateral mixing   (dyn_ldf_lap & dyn_ldf_blp routines)
20   USE dynldf_iso     ! lateral mixing                 (dyn_ldf_iso routine )
21   USE trd_oce        ! trends: ocean variables
22   USE trddyn         ! trend manager: dynamics   (trd_dyn      routine)
23   !
24   USE prtctl         ! Print control
25   USE in_out_manager ! I/O manager
26   USE lib_mpp        ! distribued memory computing library
27   USE lbclnk         ! ocean lateral boundary conditions (or mpp link)
28   USE timing         ! Timing
29
30   IMPLICIT NONE
31   PRIVATE
32
33   PUBLIC   dyn_ldf       ! called by step module
34   PUBLIC   dyn_ldf_init  ! called by opa  module
35
36   !! * Substitutions
37#  include "vectopt_loop_substitute.h90"
38   !!----------------------------------------------------------------------
39   !! NEMO/OCE 4.0 , NEMO Consortium (2018)
40   !! $Id$
41   !! Software governed by the CeCILL license (see ./LICENSE)
42   !!----------------------------------------------------------------------
43CONTAINS
44
45   SUBROUTINE dyn_ldf( kt )
46      !!----------------------------------------------------------------------
47      !!                  ***  ROUTINE dyn_ldf  ***
48      !!
49      !! ** Purpose :   compute the lateral ocean dynamics physics.
50      !!----------------------------------------------------------------------
51      INTEGER, INTENT(in) ::   kt   ! ocean time-step index
52      !
53      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::   ztrdu, ztrdv
54      !!----------------------------------------------------------------------
55      !
56      IF( ln_timing )   CALL timing_start('dyn_ldf')
57      !
58      IF( l_trddyn )   THEN                      ! temporary save of momentum trends
59         ALLOCATE( ztrdu(jpi,jpj,jpk) , ztrdv(jpi,jpj,jpk) )
60         ztrdu(:,:,:) = ua(:,:,:) 
61         ztrdv(:,:,:) = va(:,:,:) 
62      ENDIF
63
64      SELECT CASE ( nldf_dyn )                   ! compute lateral mixing trend and add it to the general trend
65      !
66      CASE ( np_lap   )    ;   CALL dyn_ldf_lap( kt, ub, vb, ua, va, 1 )      ! iso-level    laplacian
67      CASE ( np_lap_i )    ;   CALL dyn_ldf_iso( kt )                         ! rotated      laplacian
68      CASE ( np_blp   )    ;   CALL dyn_ldf_blp( kt, ub, vb, ua, va    )      ! iso-level bi-laplacian
69      !
70      END SELECT
71
72      IF( l_trddyn ) THEN                        ! save the horizontal diffusive trends for further diagnostics
73         ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:)
74         ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:)
75         CALL trd_dyn( ztrdu, ztrdv, jpdyn_ldf, kt )
76         DEALLOCATE ( ztrdu , ztrdv )
77      ENDIF
78      !                                          ! print sum trends (used for debugging)
79      IF(ln_ctl)   CALL prt_ctl( tab3d_1=ua, clinfo1=' ldf  - Ua: ', mask1=umask,   &
80         &                       tab3d_2=va, clinfo2=       ' Va: ', mask2=vmask, clinfo3='dyn' )
81      !
82      IF( ln_timing )   CALL timing_stop('dyn_ldf')
83      !
84   END SUBROUTINE dyn_ldf
85
86
87   SUBROUTINE dyn_ldf_init
88      !!----------------------------------------------------------------------
89      !!                  ***  ROUTINE dyn_ldf_init  ***
90      !!
91      !! ** Purpose :   initializations of the horizontal ocean dynamics physics
92      !!----------------------------------------------------------------------
93      !
94      IF(lwp) THEN                     !==  Namelist print  ==!
95         WRITE(numout,*)
96         WRITE(numout,*) 'dyn_ldf_init : Choice of the lateral diffusive operator on dynamics'
97         WRITE(numout,*) '~~~~~~~~~~~~'
98         WRITE(numout,*) '   Namelist namdyn_ldf: already read in ldfdyn module'
99         WRITE(numout,*) '      see ldf_dyn_init report for lateral mixing parameters'
100         WRITE(numout,*)
101         !
102         SELECT CASE( nldf_dyn )             ! print the choice of operator
103         CASE( np_no_ldf )   ;   WRITE(numout,*) '   ==>>>   NO lateral viscosity'
104         CASE( np_lap    )   ;   WRITE(numout,*) '   ==>>>   iso-level laplacian operator'
105         CASE( np_lap_i  )   ;   WRITE(numout,*) '   ==>>>   rotated laplacian operator with iso-level background'
106         CASE( np_blp    )   ;   WRITE(numout,*) '   ==>>>   iso-level bi-laplacian operator'
107         END SELECT
108      ENDIF
109      !
110   END SUBROUTINE dyn_ldf_init
111
112   !!======================================================================
113END MODULE dynldf
Note: See TracBrowser for help on using the repository browser.