1 | !----------------------------------------------------------------------- |
---|
2 | &namdyn_ldf ! lateral diffusion on momentum (default: NO selection) |
---|
3 | !----------------------------------------------------------------------- |
---|
4 | ! ! Type of the operator : |
---|
5 | ln_dynldf_OFF = .false. ! No operator (i.e. no explicit diffusion) |
---|
6 | ln_dynldf_lap = .false. ! laplacian operator |
---|
7 | ln_dynldf_blp = .false. ! bilaplacian operator |
---|
8 | ! ! Direction of action : |
---|
9 | ln_dynldf_lev = .false. ! iso-level |
---|
10 | ln_dynldf_hor = .false. ! horizontal (geopotential) |
---|
11 | ln_dynldf_iso = .false. ! iso-neutral (lap only) |
---|
12 | ! ! Coefficient |
---|
13 | nn_ahm_ijk_t = 0 ! space/time variation of eddy coefficient : |
---|
14 | ! ! =-30 read in eddy_viscosity_3D.nc file |
---|
15 | ! ! =-20 read in eddy_viscosity_2D.nc file |
---|
16 | ! ! = 0 constant |
---|
17 | ! ! = 10 F(k)=c1d |
---|
18 | ! ! = 20 F(i,j)=F(grid spacing)=c2d |
---|
19 | ! ! = 30 F(i,j,k)=c2d*c1d |
---|
20 | ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) |
---|
21 | ! ! = 32 F(i,j,k)=F(local gridscale and deformation rate) |
---|
22 | ! ! time invariant coefficients : ahm = 1/2 Uv*Lv (lap case) |
---|
23 | ! ! or = 1/12 Uv*Lv^3 (blp case) |
---|
24 | rn_Uv = 0.1 ! lateral viscous velocity [m/s] (nn_ahm_ijk_t= 0, 10, 20, 30) |
---|
25 | rn_Lv = 10.e+3 ! lateral viscous length [m] (nn_ahm_ijk_t= 0, 10) |
---|
26 | ! ! Smagorinsky settings (nn_ahm_ijk_t= 32) : |
---|
27 | rn_csmc = 3.5 ! Smagorinsky constant of proportionality |
---|
28 | rn_minfac = 1.0 ! multiplier of theorectical lower limit |
---|
29 | rn_maxfac = 1.0 ! multiplier of theorectical upper limit |
---|
30 | ! ! iso-neutral laplacian operator (ln_dynldf_iso=T) : |
---|
31 | rn_ahm_b = 0.0 ! background eddy viscosity [m2/s] |
---|
32 | / |
---|