Local variables to be added
integer          ifac, izone, ii
integer          ilelt, nlelt
double precision uref2, d2s3
double precision xkent, xeent
integer, allocatable, dimension(:) :: lstelt
Initialization and finalization
Initialization and finalization is similar to that of the base examples
Example 1
Definition of a fuel flow inlet for each face of colour 11
do ilelt = 1, nlelt
  ifac = lstelt(ilelt)
  izone = 1
  ientfu(izone) = 1
  qimp(izone)   = 2.62609d-4 / 72.d0
  rcodcl(ifac,
iw,1) = 21.47d0
    uref2 = rcodcl(ifac,
iu,1)**2                           &
    uref2 = max(uref2,1.d-12)
    
        xkent, xeent )
      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
    elseif (
iturb.eq.50) 
then 
      rcodcl(ifac,
ik,1)   = xkent
      rcodcl(ifac,
iep,1)  = xeent
      rcodcl(ifac,
iphi,1) = d2s3
      rcodcl(ifac,
ifb,1)  = 0.d0
    elseif (
iturb.eq.60) 
then 
      rcodcl(ifac,
ik,1)   = xkent
      rcodcl(ifac,
iomg,1) = xeent/
cmu/xkent
    elseif (
iturb.eq.70) 
then 
      rcodcl(ifac,
inusa,1) = 
cmu*xkent**2/xeent
    endif
  endif
enddo
  
Example 2
Definition of an air flow inlet for each face of colour 21
do ilelt = 1, nlelt
  ifac = lstelt(ilelt)
  izone = 2
  ientox(izone) = 1
  qimp(izone)   = 4.282d-3 / 72.d0
  rcodcl(ifac,
iw,1) = 0.097d0
    uref2 = rcodcl(ifac,
iu,1)**2                           &
    uref2 = max(uref2,1.d-12)
    
        xkent, xeent )
      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
    elseif (
iturb.eq.50) 
then 
      rcodcl(ifac,
ik,1)   = xkent
      rcodcl(ifac,
iep,1)  = xeent
      rcodcl(ifac,
iphi,1) = d2s3
      rcodcl(ifac,
ifb,1)  = 0.d0
    elseif (
iturb.eq.60) 
then 
      rcodcl(ifac,
ik,1)   = xkent
      rcodcl(ifac,
iomg,1) = xeent/
cmu/xkent
    elseif (
iturb.eq.70) 
then 
      rcodcl(ifac,
inusa,1) = 
cmu*xkent**2/xeent
    endif
  endif
enddo
  
Example 3
Definition of a wall for each face of colour 51 up to 59
do ilelt = 1, nlelt
  ifac = lstelt(ilelt)
  izone = 3
enddo
  
Example 4
Definition of an exit for each face of colour 91
do ilelt = 1, nlelt
  ifac = lstelt(ilelt)
  izone = 4
enddo
  
Example 5
Definition of symmetric boundary conditions for each face of colour 41 and 4.
do ilelt = 1, nlelt
  ifac = lstelt(ilelt)
  izone = 5
enddo