| [3] | 1 | MODULE ldftra_oce |
|---|
| 2 | !!===================================================================== |
|---|
| 3 | !! *** MODULE ldftra_oce *** |
|---|
| 4 | !! Ocean physics : lateral tracer mixing coefficient defined in memory |
|---|
| 5 | !!===================================================================== |
|---|
| [1601] | 6 | !! History : 9.0 ! 02-11 (G. Madec) Original code |
|---|
| [3] | 7 | !!---------------------------------------------------------------------- |
|---|
| 8 | USE par_oce ! ocean parameters |
|---|
| 9 | |
|---|
| 10 | IMPLICIT NONE |
|---|
| 11 | PRIVATE |
|---|
| 12 | |
|---|
| 13 | !!---------------------------------------------------------------------- |
|---|
| 14 | !! Lateral eddy diffusivity coefficients (tracers) |
|---|
| 15 | !!---------------------------------------------------------------------- |
|---|
| [1601] | 16 | ! !!* Namelist namtra_ldf : lateral mixing * |
|---|
| 17 | LOGICAL , PUBLIC :: ln_traldf_lap = .TRUE. !: laplacian operator |
|---|
| 18 | LOGICAL , PUBLIC :: ln_traldf_bilap = .FALSE. !: bilaplacian operator |
|---|
| 19 | LOGICAL , PUBLIC :: ln_traldf_level = .FALSE. !: iso-level direction |
|---|
| 20 | LOGICAL , PUBLIC :: ln_traldf_hor = .FALSE. !: horizontal (geopotential) direction |
|---|
| 21 | LOGICAL , PUBLIC :: ln_traldf_iso = .TRUE. !: iso-neutral direction |
|---|
| [2236] | 22 | LOGICAL , PUBLIC :: ln_traldf_grif = .FALSE. !: griffies skew flux |
|---|
| [1601] | 23 | REAL(wp), PUBLIC :: rn_aht_0 = 2000._wp !: lateral eddy diffusivity (m2/s) |
|---|
| 24 | REAL(wp), PUBLIC :: rn_ahtb_0 = 0._wp !: lateral background eddy diffusivity (m2/s) |
|---|
| 25 | REAL(wp), PUBLIC :: rn_aeiv_0 = 2000._wp !: eddy induced velocity coefficient (m2/s) |
|---|
| [2236] | 26 | REAL(wp), PUBLIC :: rn_slpmax = 0.01_wp !: slope limit |
|---|
| [3] | 27 | |
|---|
| [1601] | 28 | REAL(wp), PUBLIC :: aht0, ahtb0, aeiv0 !!: OLD namelist names |
|---|
| [2236] | 29 | REAL(wp), PUBLIC :: slpmax !: slope limit |
|---|
| [3] | 30 | |
|---|
| 31 | #if defined key_traldf_c3d |
|---|
| [1601] | 32 | REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: ahtt, ahtu, ahtv, ahtw !: ** 3D coefficients ** at T-, U-, V-, W-points |
|---|
| [3] | 33 | #elif defined key_traldf_c2d |
|---|
| [1601] | 34 | REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: ahtt, ahtu, ahtv, ahtw !: ** 2D coefficients ** at T-, U-, V-, W-points |
|---|
| [3] | 35 | #elif defined key_traldf_c1d |
|---|
| [1601] | 36 | REAL(wp), PUBLIC, DIMENSION(jpk) :: ahtt, ahtu, ahtv, ahtw !: ** 1D coefficients ** at T-, U-, V-, W-points |
|---|
| [3] | 37 | #else |
|---|
| [1601] | 38 | REAL(wp), PUBLIC :: ahtt, ahtu, ahtv, ahtw !: ** 0D coefficients ** at T-, U-, V-, W-points |
|---|
| [3] | 39 | #endif |
|---|
| 40 | |
|---|
| 41 | |
|---|
| 42 | #if defined key_traldf_eiv |
|---|
| 43 | !!---------------------------------------------------------------------- |
|---|
| 44 | !! 'key_traldf_eiv' eddy induced velocity |
|---|
| 45 | !!---------------------------------------------------------------------- |
|---|
| [2236] | 46 | LOGICAL, PUBLIC, PARAMETER :: lk_traldf_eiv = .TRUE. !: eddy induced velocity flag |
|---|
| [3] | 47 | |
|---|
| 48 | # if defined key_traldf_c3d |
|---|
| [1601] | 49 | REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: aeiu, aeiv, aeiw !: ** 3D coefficients ** at U-, V-, W-points [m2/s] |
|---|
| [3] | 50 | # elif defined key_traldf_c2d |
|---|
| [1601] | 51 | REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: aeiu, aeiv, aeiw !: ** 2D coefficients ** at U-, V-, W-points [m2/s] |
|---|
| [3] | 52 | # elif defined key_traldf_c1d |
|---|
| [1601] | 53 | REAL(wp), PUBLIC, DIMENSION(jpk) :: aeiu, aeiv, aeiw !: ** 1D coefficients ** at U-, V-, W-points [m2/s] |
|---|
| [3] | 54 | # else |
|---|
| [1601] | 55 | REAL(wp), PUBLIC :: aeiu, aeiv, aeiw !: ** 0D coefficients ** at U-, V-, W-points [m2/s] |
|---|
| [3] | 56 | # endif |
|---|
| 57 | # if defined key_diaeiv |
|---|
| [1601] | 58 | REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: u_eiv, v_eiv, w_eiv !: eddy induced velocity [m/s] |
|---|
| [3] | 59 | # endif |
|---|
| 60 | |
|---|
| 61 | #else |
|---|
| 62 | !!---------------------------------------------------------------------- |
|---|
| 63 | !! Default option : NO eddy induced velocity |
|---|
| 64 | !!---------------------------------------------------------------------- |
|---|
| [32] | 65 | LOGICAL , PUBLIC, PARAMETER :: lk_traldf_eiv = .FALSE. !: eddy induced velocity flag |
|---|
| [1601] | 66 | REAL(wp), PUBLIC :: aeiu, aeiv, aeiw !: eddy induced coef. (not used) |
|---|
| [3] | 67 | #endif |
|---|
| 68 | |
|---|
| 69 | !!---------------------------------------------------------------------- |
|---|
| [1601] | 70 | !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009) |
|---|
| 71 | !! $Id$ |
|---|
| [2236] | 72 | !! Software governed by the CeCILL licence (NEMOGCM/License_CeCILL.txt) |
|---|
| [1601] | 73 | !!===================================================================== |
|---|
| [3] | 74 | END MODULE ldftra_oce |
|---|
