Local variables to be added
integer          ifac, ii
integer          izone
integer          ilelt, nlelt
double precision d2s3
double precision zref, xuref
double precision ustar, rugd, rugt
double precision zent, xuent, xvent
double precision xkent, xeent
integer, allocatable, dimension(:) :: lstelt
Initialization and finalization
Initialization and finalization is similar to that of the base examples
Example 1
For boundary faces of color 11, assign an inlet boundary condition prescribed from the meteo profile with automatic choice between inlet/ outlet according to the meteo profile.
izone = 1
do ilelt = 1, nlelt
  ifac = lstelt(ilelt)
  
  
  
  
  iprofm(izone) = 1
enddo
  
Example 2
For boundary faces of color 21, assign an inlet boundary condition prescribed from the meteo profile.
izone = 2
do ilelt = 1, nlelt
  ifac = lstelt(ilelt)
  
  
  iprofm(izone) = 1
  
enddo
  
Example 3
For boundary faces of color 31, assign an inlet boundary condition prescribed from the meteo profile except for dynamical variables which are prescribed with a rough log law.
izone = 3
do ilelt = 1, nlelt
  ifac = lstelt(ilelt)
  
  
  iprofm(izone) = 1
  
  ustar=
xkappa*xuref/log((zref+rugd)/rugd)
  xuent=ustar/
xkappa*log((zent+rugd)/rugd)
  xvent = 0.d0
  xeent=ustar**3/
xkappa/(zent+rugd)
  rcodcl(ifac,
iu,1) = xuent
  rcodcl(ifac,
iv,1) = xvent
  
    rcodcl(ifac,
ik,1)  = xkent
    rcodcl(ifac,
iep,1) = xeent
    rcodcl(ifac,
ir11,1) = d2s3*xkent
    rcodcl(ifac,
ir22,1) = d2s3*xkent
    rcodcl(ifac,
ir33,1) = d2s3*xkent
    rcodcl(ifac,
ir12,1) = 0.d0
    rcodcl(ifac,
ir13,1) = 0.d0
    rcodcl(ifac,
ir23,1) = 0.d0
    rcodcl(ifac,
iep,1)  = xeent
    rcodcl(ifac,
ik,1)   = xkent
    rcodcl(ifac,
iep,1)  = xeent
    rcodcl(ifac,
iphi,1) = d2s3
    rcodcl(ifac,
ifb,1)  = 0.d0
    rcodcl(ifac,
ik,1)   = xkent
    rcodcl(ifac,
iomg,1) = xeent/
cmu/xkent
    rcodcl(ifac,
inusa,1) = 
cmu*xkent**2/xeent
  endif
enddo
  
Example 4
Prescribe at boundary faces of color '12' an outlet.
izone = 5
do ilelt = 1, nlelt
  ifac = lstelt(ilelt)
  
  
  
  
  
  
  
    
    
    
    
    
    
    
    
    
  endif
enddo
  
Example 5
Prescribe at boundary faces of color 4 a symmetry.
izone = 6
do ilelt = 1, nlelt
  ifac = lstelt(ilelt)
  
enddo