[11305] | 1 | MODULE par_abl |
---|
| 2 | !!====================================================================== |
---|
| 3 | !! *** par_abl *** |
---|
| 4 | !! ABL : set the Atmospheric Boundary Layer parameters |
---|
| 5 | !!====================================================================== |
---|
| 6 | !! History : 4.0 ! 2019-03 (F. Lemarié & G. Samson) Original code |
---|
| 7 | !!---------------------------------------------------------------------- |
---|
| 8 | USE par_kind ! kind parameters |
---|
| 9 | USE par_oce ! to access horizontal domain size & mpp decomposition |
---|
| 10 | |
---|
| 11 | IMPLICIT NONE |
---|
| 12 | PUBLIC |
---|
| 13 | |
---|
| 14 | !!--------------------------------------------------------------------- |
---|
| 15 | !! Active tracer parameters |
---|
| 16 | !!--------------------------------------------------------------------- |
---|
| 17 | INTEGER , PUBLIC, PARAMETER :: jptq = 2 !: Number of active tracers (=2, i.e. T & q ) |
---|
| 18 | INTEGER , PUBLIC, PARAMETER :: jp_ta = 1 !: indice for temperature |
---|
| 19 | INTEGER , PUBLIC, PARAMETER :: jp_qa = 2 !: indice for humidity |
---|
| 20 | INTEGER , PUBLIC, PARAMETER :: jptime = 2 !: number of time indices stored in memory |
---|
[11363] | 21 | |
---|
[11305] | 22 | !!--------------------------------------------------------------------- |
---|
| 23 | !! namABL Namelist options |
---|
| 24 | !!--------------------------------------------------------------------- |
---|
| 25 | INTEGER , PUBLIC :: nn_amxl !: mixing length option |
---|
| 26 | INTEGER , PUBLIC :: nn_dyn_restore !: restoring option for dynamical ABL variables |
---|
| 27 | LOGICAL , PUBLIC :: ln_geos_winds !: large-scale restoring of ABL winds toward geostrophic winds |
---|
| 28 | LOGICAL , PUBLIC :: ln_hpgls_frc !: forcing of ABL winds by large-scale pressure gradient |
---|
[11322] | 29 | LOGICAL , PUBLIC :: ln_smth_pblh !: smoothing of atmospheric PBL height |
---|
[11305] | 30 | |
---|
[11858] | 31 | CHARACTER(len=256), PUBLIC :: cn_ablrst_in !: suffix of abl restart name (input) |
---|
| 32 | CHARACTER(len=256), PUBLIC :: cn_ablrst_out !: suffix of abl restart name (output) |
---|
| 33 | CHARACTER(len=256), PUBLIC :: cn_ablrst_indir !: abl restart input directory |
---|
| 34 | CHARACTER(len=256), PUBLIC :: cn_ablrst_outdir !: abl restart output directory |
---|
| 35 | |
---|
[11305] | 36 | !!--------------------------------------------------------------------- |
---|
| 37 | !! ABL parameters for TKE turbulent closure |
---|
| 38 | !!--------------------------------------------------------------------- |
---|
| 39 | REAL(wp), PUBLIC, PARAMETER :: tke_min = 1.e-6_wp !: minimum TKE |
---|
| 40 | REAL(wp), PUBLIC, PARAMETER :: avm_bak = 1.e-4_wp !: background viscosity |
---|
| 41 | REAL(wp), PUBLIC, PARAMETER :: avt_bak = 1.e-5_wp !: background diffusion |
---|
| 42 | REAL(wp), PUBLIC, PARAMETER :: itvref = 1.0_wp / 288.0_wp !: inverse of reference virtual temperature |
---|
| 43 | !++ TKE closure parameters |
---|
| 44 | REAL(wp), PUBLIC, PARAMETER :: rn_phimax = (1._wp - 2.2_wp) / 2.2_wp !: maximum value for Ri * mxl^2 * N^2 / tke in phiz computation |
---|
| 45 | REAL(wp), PUBLIC, PARAMETER :: rn_Cek = 258._wp !: Ekman constant for Richardson number |
---|
| 46 | REAL(wp), PUBLIC, PARAMETER :: rn_epssfc = 1._wp / ( 1._wp + 2.8_wp * 2.8_wp ) |
---|
| 47 | REAL(wp), PUBLIC :: rn_ceps !: namelist parameter |
---|
| 48 | REAL(wp), PUBLIC :: rn_cm !: namelist parameter |
---|
| 49 | REAL(wp), PUBLIC :: rn_ct !: namelist parameter |
---|
| 50 | REAL(wp), PUBLIC :: rn_ce !: namelist parameter |
---|
| 51 | REAL(wp), PUBLIC :: rn_Rod !: namelist parameter |
---|
| 52 | REAL(wp), PUBLIC :: rn_Sch |
---|
| 53 | REAL(wp), PUBLIC :: mxl_min |
---|
| 54 | REAL(wp), PUBLIC :: rn_ldyn_min !: namelist parameter |
---|
| 55 | REAL(wp), PUBLIC :: rn_ldyn_max !: namelist parameter |
---|
| 56 | REAL(wp), PUBLIC :: rn_ltra_min !: namelist parameter |
---|
| 57 | REAL(wp), PUBLIC :: rn_ltra_max !: namelist parameter |
---|
| 58 | REAL(wp), PUBLIC :: rn_Ric !: critical Richardson number |
---|
| 59 | |
---|
| 60 | !!--------------------------------------------------------------------- |
---|
| 61 | !! ABL parameters for the vertical profile of the restoring term |
---|
| 62 | !!--------------------------------------------------------------------- |
---|
| 63 | REAL(wp), PUBLIC, PARAMETER :: jp_bmin = 0.5_wp |
---|
| 64 | REAL(wp), PUBLIC, PARAMETER :: jp_bmax = 1.5_wp |
---|
| 65 | REAL(wp), PUBLIC :: jp_alp0_tra |
---|
| 66 | REAL(wp), PUBLIC :: jp_alp1_tra |
---|
| 67 | REAL(wp), PUBLIC :: jp_alp2_tra |
---|
| 68 | REAL(wp), PUBLIC :: jp_alp3_tra |
---|
| 69 | REAL(wp), PUBLIC :: jp_alp0_dyn |
---|
| 70 | REAL(wp), PUBLIC :: jp_alp1_dyn |
---|
| 71 | REAL(wp), PUBLIC :: jp_alp2_dyn |
---|
| 72 | REAL(wp), PUBLIC :: jp_alp3_dyn |
---|
| 73 | REAL(wp), PUBLIC :: jp_pblh_min |
---|
| 74 | REAL(wp), PUBLIC :: jp_pblh_max |
---|
| 75 | ! parameter for the semi-implicit treatment of Coriolis term |
---|
| 76 | REAL(wp), PUBLIC, PARAMETER :: gamma_Cor = 0.55_wp |
---|
[11363] | 77 | ! ABL timestep |
---|
| 78 | REAL(wp), PUBLIC :: rdt_abl |
---|
[11305] | 79 | |
---|
| 80 | !!--------------------------------------------------------------------- |
---|
| 81 | !! ABL parameters for the diagnostic mixing length option nn_amxl = 1 |
---|
| 82 | !!--------------------------------------------------------------------- |
---|
| 83 | REAL(wp), PUBLIC, PARAMETER :: amx1 = 4.3995604393911742_wp |
---|
| 84 | REAL(wp), PUBLIC, PARAMETER :: amx2 = -18.159100102732943_wp |
---|
| 85 | REAL(wp), PUBLIC, PARAMETER :: amx3 = 40.241226956967239_wp |
---|
| 86 | REAL(wp), PUBLIC, PARAMETER :: amx4 = -43.603409583363678_wp |
---|
| 87 | REAL(wp), PUBLIC, PARAMETER :: amx5 = 17.121715347554314_wp |
---|
| 88 | REAL(wp), PUBLIC, PARAMETER :: bmx1 = -16.262675447730114_wp |
---|
| 89 | REAL(wp), PUBLIC, PARAMETER :: bmx2 = 85.088728134110781_wp |
---|
| 90 | REAL(wp), PUBLIC, PARAMETER :: bmx3 = -193.46548261141191_wp |
---|
| 91 | REAL(wp), PUBLIC, PARAMETER :: bmx4 = 196.71548261141191_wp |
---|
| 92 | REAL(wp), PUBLIC, PARAMETER :: bmx5 = -72.076052686380677_wp |
---|
| 93 | |
---|
| 94 | !!---------------------------------------------------------------------- |
---|
| 95 | !! NEMO/ABL 4.0 , NEMO Consortium (2018) |
---|
| 96 | !! $Id: sbc_oce.F90 10882 2019-04-17 15:40:17Z clem $ |
---|
| 97 | !! Software governed by the CeCILL license (see ./LICENSE) |
---|
| 98 | !!====================================================================== |
---|
| 99 | END MODULE par_abl |
---|