Compute friction velocity u* and surface sensible heat flux q0 for a non neutral atmospheric surface layer using the explicit formula developed for the ECMWF by Louis (1982)
More...
subroutine atmcls (ifac, utau, rough_d, duplus, dtplus, yplus_t, uet, gredu, cfnns, cfnnk, cfnne, dlmo, temp, totwt, liqwt, icodcl, rcodcl)
Compute friction velocity u* and surface sensible heat flux q0 for a non neutral atmospheric surface layer using the explicit formula developed for the ECMWF by Louis (1982)
◆ atmcls()
subroutine atmcls
(
integer ifac ,
double precision utau ,
double precision rough_d ,
double precision duplus ,
double precision dtplus ,
double precision yplus_t ,
double precision uet ,
double precision gredu ,
double precision cfnns ,
double precision cfnnk ,
double precision cfnne ,
double precision dlmo ,
double precision temp ,
double precision totwt ,
double precision liqwt ,
integer, dimension(nfabor,nvar) icodcl ,
double precision, dimension(nfabor,nvar,3) rcodcl )
Parameters
[in] ifac treated boundary face
[in] utau tangential mean velocity
[in] rough_d roughness z0
[in] duplus 1 over dimensionless velocity in neutral conditions
[in] dtplus 1 over dimensionless temperature in neutral conditions
[in] yplus_t thermal dimensionless wall distance
[out] uet friction velocity
[out] gredu reduced gravity for non horizontal wall
[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] dlmo inverse Monin Obukhov length (for log only)
[in] temp potential temperature in boundary cell
[in] totwt total water content in boundary cell
[in] liqwt liquid water content in boundary cell
[in] icodcl face boundary condition code:
1 Dirichlet
2 Radiative outlet
3 Neumann
4 sliding and
5 smooth wall and
6 rough wall and
9 free inlet/outlet (input mass flux blocked to 0)
11 Boundary value related to the next cell value by an affine function
13 Dirichlet for the advection operator and Neumann for the diffusion operator
[in] rcodcl boundary condition values:
rcodcl(1) value of the dirichlet
rcodcl(2) value of the exterior exchange coefficient (infinite if no exchange)
rcodcl(3) value flux density (negative if gain) in w/m2
for the velocity
for the pressure
for a scalar