[3] | 1 | MODULE zdf_oce |
---|
| 2 | !!====================================================================== |
---|
| 3 | !! *** MODULE zdf_oce *** |
---|
| 4 | !! Ocean physics : define vertical mixing variables |
---|
| 5 | !!===================================================================== |
---|
[1492] | 6 | !! history : 1.0 ! 2002-06 (G. Madec) Original code |
---|
| 7 | !! 3.2 ! 2009-07 (G.Madec) addition of avm |
---|
[3] | 8 | !!---------------------------------------------------------------------- |
---|
[2715] | 9 | USE par_oce ! ocean parameters |
---|
| 10 | USE in_out_manager ! I/O manager |
---|
| 11 | USE lib_mpp ! MPP library |
---|
[3] | 12 | |
---|
| 13 | IMPLICIT NONE |
---|
| 14 | PRIVATE |
---|
| 15 | |
---|
[2715] | 16 | PUBLIC zdf_oce_alloc ! Called in nemogcm.F90 |
---|
| 17 | |
---|
[3] | 18 | #if defined key_zdfcst || defined key_esopa |
---|
[1492] | 19 | LOGICAL, PARAMETER, PUBLIC :: lk_zdfcst = .TRUE. !: constant vertical mixing flag |
---|
[3] | 20 | #else |
---|
[1492] | 21 | LOGICAL, PARAMETER, PUBLIC :: lk_zdfcst = .FALSE. !: constant vertical mixing flag |
---|
[3] | 22 | #endif |
---|
| 23 | |
---|
[4147] | 24 | ! !!* namelist namzdf: vertical diffusion * |
---|
| 25 | REAL(wp), PUBLIC :: rn_avm0 !: vertical eddy viscosity (m2/s) |
---|
| 26 | REAL(wp), PUBLIC :: rn_avt0 !: vertical eddy diffusivity (m2/s) |
---|
| 27 | INTEGER , PUBLIC :: nn_avb !: constant or profile background on avt (=0/1) |
---|
| 28 | INTEGER , PUBLIC :: nn_havtb !: horizontal shape or not for avtb (=0/1) |
---|
| 29 | LOGICAL , PUBLIC :: ln_zdfexp !: explicit vertical diffusion scheme flag |
---|
| 30 | INTEGER , PUBLIC :: nn_zdfexp !: number of sub-time step (explicit time stepping) |
---|
| 31 | LOGICAL , PUBLIC :: ln_zdfevd !: convection: enhanced vertical diffusion flag |
---|
| 32 | INTEGER , PUBLIC :: nn_evdm !: =0/1 flag to apply enhanced avm or not |
---|
| 33 | REAL(wp), PUBLIC :: rn_avevd !: vertical eddy coeff. for enhanced vert. diff. (m2/s) |
---|
| 34 | LOGICAL , PUBLIC :: ln_zdfnpc !: convection: non-penetrative convection flag |
---|
| 35 | INTEGER , PUBLIC :: nn_npc !: non penetrative convective scheme call frequency |
---|
| 36 | INTEGER , PUBLIC :: nn_npcp !: non penetrative convective scheme print frequency |
---|
[3] | 37 | |
---|
[1537] | 38 | |
---|
[2715] | 39 | REAL(wp), PUBLIC, SAVE, ALLOCATABLE, DIMENSION(:) :: avmb , avtb !: background profile of avm and avt |
---|
| 40 | REAL(wp), PUBLIC, SAVE, ALLOCATABLE, DIMENSION(:,:) :: avtb_2d !: horizontal shape of background Kz profile |
---|
| 41 | REAL(wp), PUBLIC, SAVE, ALLOCATABLE, DIMENSION(:,:) :: bfrua, bfrva !: Bottom friction coefficients set in zdfbfr |
---|
[4990] | 42 | REAL(wp), PUBLIC, SAVE, ALLOCATABLE, DIMENSION(:,:) :: tfrua, tfrva !: top friction coefficients set in zdfbfr |
---|
[2715] | 43 | REAL(wp), PUBLIC, SAVE, ALLOCATABLE, DIMENSION(:,:,:) :: avmu , avmv !: vertical viscosity coef at uw- & vw-pts [m2/s] |
---|
| 44 | REAL(wp), PUBLIC, SAVE, ALLOCATABLE, DIMENSION(:,:,:) :: avm , avt !: vertical viscosity & diffusivity coef at w-pt [m2/s] |
---|
[7256] | 45 | REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: avt_k , avm_k ! not enhanced Kz |
---|
| 46 | REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: avmu_k, avmv_k ! not enhanced Kz |
---|
| 47 | REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: en !: now turbulent kinetic energy [m2/s2] |
---|
[3] | 48 | |
---|
[1492] | 49 | !!---------------------------------------------------------------------- |
---|
[2715] | 50 | !! NEMO/OPA 4.0 , NEMO Consortium (2011) |
---|
[1492] | 51 | !! $Id$ |
---|
[2715] | 52 | !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) |
---|
| 53 | !!---------------------------------------------------------------------- |
---|
| 54 | CONTAINS |
---|
| 55 | |
---|
| 56 | INTEGER FUNCTION zdf_oce_alloc() |
---|
| 57 | !!---------------------------------------------------------------------- |
---|
| 58 | !! *** FUNCTION zdf_oce_alloc *** |
---|
| 59 | !!---------------------------------------------------------------------- |
---|
| 60 | ! |
---|
| 61 | ALLOCATE(avmb(jpk) , bfrua(jpi,jpj) , & |
---|
| 62 | & avtb(jpk) , bfrva(jpi,jpj) , avtb_2d(jpi,jpj) , & |
---|
[4990] | 63 | & tfrua(jpi, jpj), tfrva(jpi, jpj) , & |
---|
[2715] | 64 | & avmu(jpi,jpj,jpk), avm(jpi,jpj,jpk) , & |
---|
[7256] | 65 | & avmv (jpi,jpj,jpk), avt (jpi,jpj,jpk) , & |
---|
| 66 | & avt_k (jpi,jpj,jpk), avm_k (jpi,jpj,jpk) , & |
---|
| 67 | & avmu_k(jpi,jpj,jpk), avmv_k(jpi,jpj,jpk) , & |
---|
| 68 | & en (jpi,jpj,jpk), STAT = zdf_oce_alloc ) |
---|
[2715] | 69 | ! |
---|
| 70 | IF( zdf_oce_alloc /= 0 ) CALL ctl_warn('zdf_oce_alloc: failed to allocate arrays') |
---|
| 71 | ! |
---|
| 72 | END FUNCTION zdf_oce_alloc |
---|
| 73 | |
---|
[3] | 74 | !!====================================================================== |
---|
| 75 | END MODULE zdf_oce |
---|