programmer's documentation
Functions/Subroutines
atmcls.f90 File Reference

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...

Functions/Subroutines

subroutine atmcls (ifac, iel, utau, yplus, uet, gredu, q0, e0, rib, lmo, cfnnu, cfnns, cfnnk, cfnne, icodcl, rcodcl)
 

Detailed Description

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)

Function/Subroutine Documentation

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]ifactreated boundary face
[in]ielboundary cell
[in]utautangentiel mean
[in]yplusadim distance to he boundary faces
[out]uetfriction velocity
[out]gredureduced gravity for non horizontal wall
[out]q0latent heat flux
[out]e0sensible heat flux
[out]ribRichardson number
[out]lmoMonin-Obukhov length
[out]cfnnunon neutral correction coefficients for profiles of wind
[out]cfnnsnon neutral correction coefficients for profiles of scalar
[out]cfnnknon neutral correction coefficients for profiles of k
[out]cfnnenon neutral correction coefficients for profiles of eps
[out]icodclcode 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]rcodclvaleur 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