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 Dirichlet2 Radiative outlet3 Neumann4 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 function13 Dirichlet for the advection operator and Neumann for the diffusion operator  |  | [in] | rcodcl | boundary condition values: 
rcodcl(1) value of the dirichletrcodcl(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   |