Compute u*, q0, e0, (momentum, sensible heat and latent heat fluxes) for a non neutral atmospheric surface layer using the explicit formula developed for the ECMWF by Louis (1982)
More...
|
subroutine | atmcls (ifac, iel, utau, yplus, uet, gredu, q0, e0, rib, lmo, cfnnu, cfnns, cfnnk, cfnne, icodcl, rcodcl) |
|
Compute u*, q0, e0, (momentum, sensible heat and latent heat fluxes) for a non neutral atmospheric surface layer using the explicit formula developed for the ECMWF by Louis (1982)
subroutine atmcls |
( |
integer |
ifac, |
|
|
integer |
iel, |
|
|
double precision |
utau, |
|
|
double precision |
yplus, |
|
|
double precision |
uet, |
|
|
double precision |
gredu, |
|
|
double precision |
q0, |
|
|
double precision |
e0, |
|
|
double precision |
rib, |
|
|
double precision |
lmo, |
|
|
double precision |
cfnnu, |
|
|
double precision |
cfnns, |
|
|
double precision |
cfnnk, |
|
|
double precision |
cfnne, |
|
|
integer, dimension(nfabor,nvarcl) |
icodcl, |
|
|
double precision, dimension(nfabor,nvarcl,3) |
rcodcl |
|
) |
| |
- Parameters
-
[in] | ifac | treated boundary face |
[in] | iel | boundary cell |
[in] | utau | tangentiel mean |
[in] | yplus | adim distance to he boundary faces |
[out] | uet | friction velocity |
[out] | gredu | reduced gravity for non horizontal wall |
[out] | q0 | latent heat flux |
[out] | e0 | sensible heat flux |
[out] | rib | Richardson number |
[out] | lmo | Monin-Obukhov length |
[out] | cfnnu | non neutral correction coefficients for profiles of wind |
[out] | cfnns | non neutral correction coefficients for profiles of scalar |
[out] | cfnnk | non neutral correction coefficients for profiles of k |
[out] | cfnne | non neutral correction coefficients for profiles of eps |
[out] | icodcl | code for boundary conditions at boundary faces (nfabor,nvarcl)
- = 1 -> dirichlet
- = 3 -> densite de flux
- = 4 -> glissemt et u.n=0 (vitesse)
- = 5 -> frottemt et u.n=0 (vitesse)
- = 6 -> rugosite et u.n=0 (vitesse)
- = 9 -> entree/sortie libre (vitesse entrante eventuelle bloquee
|
[out] | rcodcl | valeur des conditions aux limites (nfabor,nvarcl) aux faces de bord
- rcodcl(1) = valeur du dirichlet
- rcodcl(2) = valeur du coef. d'echange ext. (infinie si pas d'echange)
- rcodcl(3) = valeur de la densite de flux (negatif si gain) w/m2 ou hauteur de rugosite (m) si icodcl=6 – pour les vitesses (vistl+visct)*gradu – pour la pression dt*gradp – pour les scalaires cp*(viscls+visct/sigmas)*gradt
|