9.0
general documentation
Loading...
Searching...
No Matches
atmcls.f90 File Reference

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

Functions/Subroutines

subroutine atmcls (ifac, utau, rough_d, duplus, dtplus, yplus_t, uet, gredu, cfnns, cfnnk, cfnne, dlmo, temp, totwt, liqwt, icodcl, rcodcl)

Detailed Description

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)

Function/Subroutine Documentation

◆ 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]ifactreated boundary face
[in]utautangential mean velocity
[in]rough_droughness z0
[in]duplus1 over dimensionless velocity in neutral conditions
[in]dtplus1 over dimensionless temperature in neutral conditions
[in]yplus_tthermal dimensionless wall distance
[out]uetfriction velocity
[out]gredureduced gravity for non horizontal wall
[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]dlmoinverse Monin Obukhov length (for log only)
[in]temppotential temperature in boundary cell
[in]totwttotal water content in boundary cell
[in]liqwtliquid water content in boundary cell
[in]icodclface boundary condition code:
  • 1 Dirichlet
  • 2 Radiative outlet
  • 3 Neumann
  • 4 sliding and $ \vect{u} \cdot \vect{n} = 0 $
  • 5 smooth wall and $ \vect{u} \cdot \vect{n} = 0 $
  • 6 rough wall and $ \vect{u} \cdot \vect{n} = 0 $
  • 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]rcodclboundary 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
    1. for the velocity $ (\mu+\mu_T)
                                       \gradv \vect{u} \cdot \vect{n}  $
    2. for the pressure $ \Delta t
                                       \grad P \cdot \vect{n}  $
    3. for a scalar $ cp \left( K +
                                        \dfrac{K_T}{\sigma_T} \right)
                                        \grad T \cdot \vect{n} $