My Project
programmer's documentation
Functions/Subroutines
attycl.f90 File Reference

Automatic boundary conditions for atmospheric module (based on meteo file) More...

Functions/Subroutines

subroutine attycl (itypfb, izfppp, icodcl, rcodcl)
 Automatically compute the boundary conditions from the meteo file or from the imbrication profiles. More...
 

Detailed Description

Automatic boundary conditions for atmospheric module (based on meteo file)

Function/Subroutine Documentation

◆ attycl()

subroutine attycl ( integer, dimension(nfabor)  itypfb,
integer, dimension(nfabor)  izfppp,
integer, dimension(nfabor,nvar icodcl,
double precision, dimension(nfabor,nvar,3)  rcodcl 
)

Automatically compute the boundary conditions from the meteo file or from the imbrication profiles.

Parameters
[in]itypfbboundary face types
[in]izfpppboundary face zone number for atmospheric module
[out]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)
  • 13 Dirichlet for the advection operator and Neumann for the diffusion operator
[out]rcodclBoundary conditions value
  • 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 or roughness in m if icodcl=6
    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} $