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.
ldftra_oce.F90 in trunk/NEMO/OPA_SRC/LDF – NEMO

source: trunk/NEMO/OPA_SRC/LDF/ldftra_oce.F90 @ 605

Last change on this file since 605 was 461, checked in by opalod, 18 years ago

nemo_v1_update_052:RB: update lateral diffusion computation following the reorganization of both dynamics and tracers

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 3.7 KB
RevLine 
[3]1MODULE ldftra_oce
2   !!=====================================================================
3   !!                      ***  MODULE  ldftra_oce  ***
4   !! Ocean physics :  lateral tracer mixing coefficient defined in memory
5   !!=====================================================================
6   !!
7   !! ** Purpose : - Define in memory lateral tracer mixing coefficients
8   !!
9   !! History :
10   !!   9.0  !  02-11  (G. Madec)  Original code (from common.h)
11   !!----------------------------------------------------------------------
[247]12   !!  OPA 9.0 , LOCEAN-IPSL (2005)
13   !! $Header$
14   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
[3]15   !!----------------------------------------------------------------------
16   !! * Modules used
17   USE par_oce         ! ocean parameters
18
19   IMPLICIT NONE
20   PRIVATE
21
22   !!----------------------------------------------------------------------
23   !! Lateral eddy diffusivity coefficients (tracers)
24   !!----------------------------------------------------------------------
25
[32]26   LOGICAL , PUBLIC ::              & !!: ** lateral mixing namelist (nam_traldf) **
27      ln_traldf_lap   = .TRUE.  ,   &  !: laplacian operator
28      ln_traldf_bilap = .FALSE. ,   &  !: bilaplacian operator
29      ln_traldf_level = .FALSE. ,   &  !: iso-level direction
30      ln_traldf_hor   = .FALSE. ,   &  !: horizontal (geopotential) direction
31      ln_traldf_iso   = .TRUE.         !: iso-neutral direction
[3]32
[32]33   REAL(wp), PUBLIC ::              & !!: ** lateral mixing namelist (namldf) **
34      aht0  = 2000._wp     ,        &  !: lateral eddy diffusivity (m2/s)
35      ahtb0 =    0._wp     ,        &  !: lateral background eddy diffusivity (m2/s)
36      aeiv0 = 2000._wp                 !: eddy induced velocity coefficient (m2/s)
[3]37
38#if defined key_traldf_c3d
[32]39   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   &  !: ** 3D coefficients **
[3]40#elif defined key_traldf_c2d
[32]41   REAL(wp), PUBLIC, DIMENSION(jpi,jpj)     ::   &  !: ** 2D coefficients **
[3]42#elif defined key_traldf_c1d
[32]43   REAL(wp), PUBLIC, DIMENSION(jpk)         ::   &  !: ** 1D coefficients **
[3]44#else
[32]45   REAL(wp), PUBLIC                         ::   &  !: ** 0D coefficients **
[3]46#endif
[32]47      ahtt, ahtu, ahtv, ahtw                !: T-, U-, V-, W-points coefficients
[3]48
49
50#if defined key_traldf_eiv
51   !!----------------------------------------------------------------------
52   !!   'key_traldf_eiv'                              eddy induced velocity
53   !!----------------------------------------------------------------------
[32]54   LOGICAL, PUBLIC, PARAMETER ::   lk_traldf_eiv   = .TRUE.   !: eddy induced velocity flag
[3]55     
56# if defined key_traldf_c3d
[32]57   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   &  !: ** 3D coefficients **
[3]58# elif defined key_traldf_c2d
[32]59   REAL(wp), PUBLIC, DIMENSION(jpi,jpj)     ::   &  !: ** 2D coefficients **
[3]60# elif defined key_traldf_c1d
[32]61   REAL(wp), PUBLIC, DIMENSION(jpk)         ::   &  !: ** 1D coefficients **
[3]62# else
[32]63   REAL(wp), PUBLIC                         ::   &  !: ** 0D coefficients **
[3]64# endif
[32]65      aeiu, aeiv, aeiw                              !: U-, V-, W-points  induced velocity coef. (m2/s)
[3]66
67# if defined key_diaeiv
[32]68   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::    &  !:
69      u_eiv, v_eiv, w_eiv     !: The three component of the eddy induced velocity (m/s)
[3]70# endif
71
72#else
73   !!----------------------------------------------------------------------
74   !!   Default option :                           NO eddy induced velocity
75   !!----------------------------------------------------------------------
[32]76   LOGICAL , PUBLIC, PARAMETER ::   lk_traldf_eiv   = .FALSE.   !: eddy induced velocity flag
[3]77   REAL(wp), PUBLIC ::   aeiu, aeiv, aeiw
78#endif
79
80   !!----------------------------------------------------------------------
81END MODULE ldftra_oce
Note: See TracBrowser for help on using the repository browser.