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
|