commande : 
 ./cs_solver --param Case-1-2 --mpi

***************************************************************

                                  (R)
                      Code_Saturne

                      Version 4.2.0


  Copyright (C) 1998-2015 EDF S.A., France

  build dim. 07 févr. 2016 13:11:58 CET
  Version MPI 3.0 (Open MPI 1.10.2)


  The Code_Saturne CFD tool  is free software;
  you can redistribute it and/or modify it under the terms
  of the GNU General Public License as published by the
  Free Software Foundation; either version 2 of the License,
  or (at your option) any later version.

  The Code_Saturne CFD tool is distributed in the hope that
  it will be useful, but WITHOUT ANY WARRANTY; without even
  the implied warranty of MERCHANTABILITY or FITNESS FOR A
  PARTICULAR PURPOSE.  See the GNU General Public License
  for more details.

***************************************************************


Configuration locale du cas :

  Date :               mer. 28 déc. 2016 18:41:24 CET
  Système :            Linux 4.4.0-28-generic
  Machine :            Atanor-Station-Calcul
  Processeur :         model name	: Intel(R) Xeon(R) CPU           E5506  @ 2.13GHz
  Mémoire :            12009 Mo
  Utilisateur :        atanor (Atanor)
  Répertoire :         /home/atanor/Documents/Tests/TestBruleurGaz/CASE1/RESU/20161228-1841
  Rangs MPI :          8 (attribut appnum : 0)
  Threads OpenMP :     1
  Processeurs/noeud :  4
  méthode de lecture : MPI-IO collectif (positions explicites)
  méthode d'écriture : MPI-IO collectif (positions explicites)
   pas des rangs E/S :  1

 Lecture des méta-données du fichier: "mesh_input"

===============================================================

                   CALCULATION PREPARATION
                   =======================


 ===========================================================



                                                             
 Pas d erreur detectee lors de la verification des donnees   
                                        (cs_user_combustion).


 No error detected during the data verification
                          cs_user_parameters.f90 and others).


 ===========================================================

               CALCULATION PARAMETERS SUMMARY
               ==============================

 -----------------------------------------------------------


 ** SPECIFIC PHYSICS:
    ----------------

 --- Diffusion Flame: 3 Point Chemistry
       OPTION =              1

 --- Janaf or not (user tabulation required in this case)
       INDJON =              1 (1: Janaf, 0: user)

 --- Combustible characteristics
       Combustible :     CH4         
       PCI =        0.50141E+08 J/kg


 ** DIMENSIONS
    ----------

 --- Physics
       NVAR   =              9 (Nb variables                )
       NSCAL  =              3 (Nb scalars                  )
       NSCAUS =              0 (Nb user scalars             )
       NSCAPP =              3 (Nb specific physics scalars )
       NPROCE =             22 (Nb cell properties          )


 -----------------------------------------------------------


 ** PHYSICAL PROPERTIES
    -------------------

       GX     =    0.00000E+00 (Gravity x component         )
       GY     =    0.00000E+00 (Gravity y component         )
       GZ     =    0.00000E+00 (Gravity z component         )

       ICORIO =          0 (Coriolis source terms       )

  -- Continuous phase:

       RO0    =    0.11989E+01 (Reference density           )
       VISCL0 =    0.18300E-04 (Ref. molecular dyn. visc.   )
       CP0    =    0.10172E+04 (Ref. specific heat          )
       ICP    =              0 (> 0: variable CP (usphyv)   )
       P0     =    0.10132E+06 (Ref. total pressure         )
       PRED0  =    0.00000E+00 (Ref. reduced pressure       )
       T0     =    0.29315E+03 (Ref. temperature            )

       IROVAR =              1 (Density constant(0) or not(1)
       IVIVAR =              0 (Molec. visc cst.(0) or not(1)

       Initial reference point for pressure
       XYZP0  =    0.00000E+00   0.00000E+00   0.00000E+00
       DIFTL0 =    0.10000E-01 (Ref. dynamic diffusivity    )


 -----------------------------------------------------------


 ** HOMOGENEOUS MIXTURE MODEL FOR CAVITATION
    ----------------------------------------

       ICAVIT =             -1 (-1: single phase flow       )
                               ( 0: no vap./cond. model     )
                               ( 1: Merkle's model        )


 -----------------------------------------------------------


 ** THERMAL MODEL
    -------------
 --- Continuous phase:

   - Commons
       ITHERM =              2 (0: no thermal model         )
                               (1: temperature              )
                               (2: enthalpy                 )
                               (3: total energy             )
       ITPSCL =              1 (0: none                     )
                               (1: temperature in Kelvin    )
                               (2: temperature in Celsius   )
       ISCALT =              3 (Thermal scalar number       )


 ** TURBULENCE
    ----------

 --- Continuous phase:

   - Commons
       ITURB  =             21 (Turbulence model            )
       IWALLF =              3 (wall function               )
                               (0: disabled                 )
                               (1: one scale power law
                                   (forbidden for k-epsilon))
                               (2: one scale log law        )
                               (3: two scales log law       )
                               (4: scalable wall function   )
                               (5: two scales V. Driest     )
                               (6: two scales smooth/rough  )
       IWALLT =              0 (Exch. coeff. correlation    )
                               (0: not activated            )
                               (1: activated                )
       YPLULI =    0.23810E+01 (Limit Y+                    )
       IGRHOK =              0 (1: computed Grad(rho k)     )

   - Linear production k-epsilon (ITURB = 21)
       ALMAX  =   -0.99900E+03 (Characteristic length       )
       UREF   =    0.10000E+01 (Characteristic velocity     )
       ICLKEP =              0 (k-epsilon clipping model    )
       IKECOU =              0 (k-epsilon coupling mode     )
       IGRAKE =              1 (Account for gravity         )
       RELAXV =    0.10000E+01 for k        (Relaxation)
       RELAXV =    0.10000E+01 for epsilon  (Relaxation)

   - Rotation/curvature correction
       IRCCOR =              0 (0: desactivated             )
                               (1: activated                )

 --- Constants

   - Commons
       XKAPPA =    0.42000E+00 (Von Karman constant         )
       CSTLOG =    0.52000E+01 (U+=Log(y+)/kappa +CSTLOG    )
       APOW   =    0.83000E+01 (U+=APOW (y+)**BPOW (W&W law))
       BPOW   =    0.14286E+00 (U+=APOW (y+)**BPOW (W&W law))

   - Linear production k-epsilon (ITURB = 21)
       Ce1    =    0.14400E+01 (Cepsilon 1: production coef.)
       CE2    =    0.19200E+01 (Cepsilon 2: dissipat.  coef.)
       SIGMAK =    0.10000E+01 (Prandtl relative to k       )
       SIGMAE =    0.13000E+01 (Prandtl relative to epsilon )
       CMU    =    0.90000E-01 (Cmu constant                )


 -----------------------------------------------------------


 ** SECONDARY VISCOSITY
    -------------------

 --- Continuous phase:
       IVISSE =              1 (1: accounted for            )


 -----------------------------------------------------------


 ** RADIATIVE THERMAL TRANSFER
    --------------------------

 --- Continuous phase:
       IIRAYO =              1 (0: no; 1: DOM; 2: P-1       )

 --- Options:
       ISUIRD =              0 (0: no restart; 1: restart   )
       NFREQR =              1 (Radiation pass frequency    )
       i_quadrature =        1 (quadrature number (if DOM)
                                 1: S4 (24 directions)
                                 2: S6 (48 directions)
                                 3: S8 (80 directions)
                                 4: T2 (32 directions)
                                 5: T4 (128 directions)
                                 6: Tn (8*ndirec^2 directions))
       IDIVER =              2 (0 1 or 2: compute radiat. ST)
       IMODAK =              0 (1: modak absor coef; 0 else )
       IIMPAR =              1 (0 1 or 2: print wall temp.  )
       IIMLUM =              0 (0 1 or 2: print solver info )
       IMOADF =              0 (0 1 or 2: none,ADF08,ADF50  )
       IMFSCK =              0 (0 1 : no FSCK, FSCK         )


 -----------------------------------------------------------


 ** TIME STEPPING
    -------------

    UNSTEADY ALGORITHM

 --- Time step parameters

       IDTVAR =              2 (0 cst; 1,2 var (t, t-space  )
       IPTLRO =              0 (1: rho-related DT clipping  )
       COUMAX =    0.10000E+01 (Maximum target CFL          )
       FOUMAX =    0.10000E+02 (Maximum target Fourier      )
       VARRDT =    0.10000E+00 (For var. DT, max. increase  )
       DTMIN  =    0.10000E-02 (Minimum time step           )
       DTMAX  =    0.10000E+02 (Maximum time step           )
       DTREF  =    0.10000E-01 (Reference time step         )

       With a non-constant time step (IDTVAR = 1 or 2),
         when the value of COUMAX or FOUMAX is negative
         or zero, the associated time step limitation (for
         CFL and Fourier respectively) is ignored.

 --- Frozen velocity field

       ICCVFG =              0 (1: frozen velocity field    )

 --- Per-variable properties

------------------------------------
 Variable          ISTAT      CDTVAR
------------------------------------
 Velocity              1  0.1000E+01
 Pressure              0  0.1000E+01
 k                     1  0.1000E+01
 epsilon               1  0.1000E+01
 mixture_fraction      1  0.1000E+01
 mixture_fraction      1  0.1000E+01
 Enthalpy              1  0.1000E+01
----------------------------

       ISTAT  =  0 ou  1       (1 for unsteady              )
       CDTVAR >  0             (time step multiplier        )

--- Relaxation coefficient
    RHO(n+1)=SRROM*RHO(n)+(1-SRROM)*RHO(n+1)
       SRROM  =    0.95000E+00

 --- Order of base time stepping scheme
       ISCHTP =              1 (1: order 1; 2: order 2      )
                                                             

 -----------------------------------------------------------


 ** CONVECTION - DIFFUSION
    ----------------------

---------------------------------------------------------------------
 Variable          ICONV  IDIFF IDIFFT ISCHCV ISSTPC   BLENCV  THETAV
---------------------------------------------------------------------
 Velocity              1      1      1      1      0 0.00E+00 0.10E+01
 Pressure              0      1      1      1      0 0.00E+00 0.10E+01
 k                     1      1      1      1      0 0.00E+00 0.10E+01
 epsilon               1      1      1      1      0 0.00E+00 0.10E+01
 mixture_fraction      1      1      1      1      0 0.10E+01 0.10E+01
 mixture_fraction      1      1      1      1      0 0.10E+01 0.10E+01
 Enthalpy              1      1      1      1      0 0.10E+01 0.10E+01
-------------------------------------------------------------

       ICONV  =  0 ou  1       (1 for convection active     )
       IDIFF  =  0 ou  1       (1 for total diffusion active)
       IDIFFT =  0 ou  1       (1 for turbulent diff. active)
       ISCHCV =  0 ou  1       (SOLU or CD                  )
       ISSTPC =  0 ou  1       (1: no slope test            )
       BLENCV =  [0.;1.]       (1-upwind proportion         )
       THETAV =  [0.;1.]       (0.5 Crank-Nicolson/AB       )
                               (theta for convection-       )
                               (diffusion terms uses        )
                               ((1-theta).old+theta.new     )


 -----------------------------------------------------------


 ** STOKES
    ------
       IDILAT =              1 (1 : without unsteady term
                                    in the continuity equation
                                2 : with unsteady term in 
                                    the continuity equation)
       IPOROS =              0 (0 : without porous media
                                1 : with porous media)
       IPHYDR =              0 (1: account for explicit
                                   balance between pressure
                                   gradient, gravity source
                                   terms, and head losses     )
       ICALHY =              0 (1: compute hydrastatic
                                   pressure for Dirichlet
                                   conditions for pressure
                                   on outlet                  )
       IPRCO  =              1 (1: pressure-continuity        )
       IPUCOU =              0 (1: reinforced U-P coupling    )
       NTERUP =              1 (n: n sweeps on navsto for
                                   velocity/pressure coupling )

  -- Continuous phase:

       IREVMC =              0 (Velocity reconstruction mode)

       RELAXV =    0.10000E+01 for pressure (relaxation)
       ARAK   =    0.10000E+01 (Arakawa factor              )

       ISTMPF =              1 (time scheme for flow
                               (0: explicit (THETFL = 0     )
                               (1: std scheme (Saturne 1.0  )
                               (2: 2nd-order (THETFL = 0.5  )
       THETFL =   -0.99900E+03 (theta for mass flow         )
       IROEXT =              0 (density extrapolation
                               (0: explicit
                               (1: n+thetro with thetro=1/2
                               (2: n+thetro with thetro=1
       THETRO =    0.00000E+00 (theta for density
                               ((1+theta).new-theta.old
       IVIEXT =              0 (total viscosity extrapolation
                               (0: explicit
                               (1: n+thetvi with thetro=1/2
                               (2: n+thetvi with thetro=1
       THETVI =    0.00000E+00 (theta for total viscosity
                               ((1+theta).new-theta.old
       ICPEXT =              0 (specific heat extrapolation
                               (0: explicit
                               (1: n+thetcp with thetro=1/2
                               (2: n+thetcp with thetro=1
       THETCP =    0.00000E+00 (specific heat theta-scheme
                               ((1+theta).new-theta.old
       THETSN =    0.00000E+00 (Nav-Stokes S.T. theta scheme)
                               ((1+theta).new-theta.old
       THETST =    0.00000E+00 (Turbulence S.T. theta-scheme)
                               ((1+theta).new-theta.old
       EPSUP  =    0.10000E-04 (Velocity/pressure coupling
                                stop test                   )


 -----------------------------------------------------------


 ** GRADIENTS CALCULATION
    ---------------------

       IMRGRA =              0 (Reconstruction mode         )
       ANOMAX =    0.78540E+00 (Non-ortho angle: limit for  )
                               (least squares ext. neighbors)

-------------------------------------------------------------------
 Variable         NSWRGR NSWRSM      EPSRGR      EPSRSM      EXTRAG
-------------------------------------------------------------------
 Velocity            100      1  0.1000E-04  0.1000E-05  0.0000E+00
 Pressure            100      2  0.1000E-04  0.1000E-05  0.0000E+00
 k                   100      1  0.1000E-04  0.1000E-05  0.0000E+00
 epsilon             100      1  0.1000E-04  0.1000E-05  0.0000E+00
 mixture_fraction    100      1  0.1000E-04  0.1000E-06  0.0000E+00
 mixture_fraction    100      1  0.1000E-04  0.1000E-06  0.0000E+00
 Enthalpy            100      1  0.1000E-04  0.1000E-06  0.0000E+00
-----------------------------------------------------------

-------------------------------------------
 Variable         IRCFLU IMLIGR      CLIMGR
-------------------------------------------
 Velocity              1     -1  0.1500E+01
 Pressure              1     -1  0.1500E+01
 k                     1     -1  0.1500E+01
 epsilon               1     -1  0.1500E+01
 mixture_fraction      1     -1  0.1500E+01
 mixture_fraction      1     -1  0.1500E+01
 Enthalpy              1     -1  0.1500E+01
-----------------------------------

       NSWRGR =                (nb sweep gradient reconstr. )
       NSWRSM =                (nb sweep rhs reconstrcution )
       EPSRGR =                (grad. reconstruction prec.  )
       EPSRSM =                (rhs   reconstruction prec.  )
       EXTRAG =  [0.;1.]       (gradients extrapolation     )
       IRCFLU =  0 ou  1       (flow reconstruction         )
       IMLIGR =  < 0, 0 ou 1   (gradient limitation method  )
       CLIMGR =  > 1 ou 1      (gradient limitation coeff.  )


 -----------------------------------------------------------


 ** FACE INTERPOLATION
    ------------------
       IMVISF =              0 (0 arithmetic                )


 -----------------------------------------------------------


 ** BASE ITERATIVE SOLVERS
    ----------------------

------------------------------------
 Variable              EPSILO IDIRCL
------------------------------------
 Velocity          0.1000E-06      1
 Pressure          0.1000E-06      1
 k                 0.1000E-06      1
 epsilon           0.1000E-06      1
 mixture_fraction  0.1000E-07      1
 mixture_fraction  0.1000E-07      1
 Enthalpy          0.1000E-07      1
------------------------------------

       IRESOL =            -1  (automatic solver choice     )
                IPOL*1000 + 0  (p conjuguate gradient       )
                            1  (Jacobi                      )
                IPOL*1000 + 2  (bicgstab                    )
                  avec IPOL    (preconditioning degree      )
       NITMAX =                (max number of iterations    )
       EPSILO =                (resolution precision        )
       IDIRCL = 0 ou 1         (shift diagonal if   
                                ISTAT=0 and no Dirichlet    )


 -----------------------------------------------------------


 ** SCALARS
    -------

       ITBRRB =              0 (T or H reconstruction at bdy)

-------------------------------------------------------------
 Variable         Number ISCACP  ITURT     VISLS0      SIGMAS
-------------------------------------------------------------
 mixture_fraction      1      0      0  0.1830E-04  0.7000E+00
 mixture_fraction      2      0      0  0.1830E-04  0.7000E+00
 Enthalpy              3      0      0  0.1830E-04  0.7000E+00
---------------------------------------------------------------------

------------------------------------
 Variable         Number      RVARFL
------------------------------------
 mixture_fraction      1  0.8000E+00
 mixture_fraction      2  0.8000E+00
 Enthalpy              3  0.8000E+00
-------------------------------------------

-------------------------------------------------------
 Variable         Number ICLVFL      SCAMIN      SCAMAX
-------------------------------------------------------
 mixture_fraction      1     -1  0.0000E+00  0.1000E+01
 mixture_fraction      2      1 -0.1000E+13  0.1000E+13
 Enthalpy              3     -1 -0.1000E+13  0.1000E+13
-------------------------------------------------------

-------------------------------------------------------------

       For each scalar, the number indicates it's rank
         in the list of all scalars. User scalars are placed
         first, from 1 to NSCAUS. Specific physics scalars
         are placed at the end, from
         NSCAUS+1 to NSCAPP+NSCAUS=NSCAL.

       ISCACP = 0 or 1     2   (use Cp or not               )
       VISLS0 = >0             (Reference viscosity         )
       SIGMAS = >0             (Schmidt                     )
       RVARFL = >0             (Rf, cf variance dissipation )
       ICLVFL = 0, 1 or 2      (Variance clipping mode      )
       SCAMIN =                (Min authorized value        )
       SCAMAX =                (Max authorized value        )
        For variances, SCAMIN is ignored and SCAMAX is used
          only if ICLVFL = 2

------------------------------------------------------
   Scalar        THETSS    IVSEXT      THETVS
------------------------------------------------------
          1  0.0000E+00         0  0.0000E+00
          2  0.0000E+00         0  0.0000E+00
          3  0.0000E+00         0  0.0000E+00
------------------------------------------------------

       THETSS =                (theta for source terms      )
                               ((1+theta).new-theta.old     )
       IVSEXT =                (extrap. total viscosity     )
                               (0: explicit                 )
                               (1: n+thetvs with thetvs=1/2 )
                               (2: n+thetvs with thetvs=1   )
       THETVS =                (theta for scalar diffusivity
                               ((1+theta).new-theta.old     )


 -----------------------------------------------------------


 ** CALCULATION MANAGEMENT
    ----------------------

 --- Restarted calculation
       ISUITE =              0 (1: restarted calculation    )
       ILEAUX =              0 (1: read  restart/auxiliary  )
       IECAUX =              1 (1: write checkpoint/auxiliary)


 --- Calculation time
     The numbering of time steps and the measure of simulated
       physical time are absolute values, and not values
       relative to the current calculation.

       INPDT0 =              0 (1: 0 time step calcuation   )
       NTMABS =           1000 (Final time step required    )

 --- CPU time margin
       TMARUS =   -0.10000E+01 (CPU time margin before stop )


 -----------------------------------------------------------


 ** INPUT-OUTPUT
    ------------

 --- Restart file
       NTSUIT =             10 (Checkpoint frequency )

 --- Post-processed variables


Velocity
Pressure
k
epsilon
mixture_fraction
mixture_fraction_variance
Enthalpy
temperature
Density
Boundary Density
TurbVisc
CourantNb
FourierNb
total_pressure
temperature
ym_fuel
ym_oxyd
ym_prod
KABS
Temp4
Temp3
Local Time Step
Flux_incident
Th_conductivity
Thickness
Emissivity
Flux_net
Flux_convectif
Coeff_ech_convectif
Yplus
TimeAverage1
TimeAverage2
         --           --

 --- Probe history files
       NTHIST =             -1 (Output frequency     )
       FRHIST =     -.10000E+01 (Output frequency (s) )
       NCAPT  =              0 (Number of probes     )
       NTHSAV =             -1 (Checkpoint frequency )

       Number Name                  Nb. probes (-1: all)
         2 Velocity[X]                   -1
         3 Velocity[Y]                   -1
         4 Velocity[Z]                   -1
         5 Pressure                      -1
         6 k                             -1
         7 epsilon                       -1
         8 mixture_fraction              -1
         9 mixture_fraction              -1
        10 Enthalpy                      -1
        11 Density                       -1
        12 LamVisc                       -1
        13 TurbVisc                      -1
        14 CourantNb                     -1
        15 FourierNb                     -1
        16 total_pressure                -1
        17 temperature                   -1
        18 ym_fuel                       -1
        19 ym_oxyd                       -1
        20 ym_prod                       -1
        21 KABS                          -1
        22 Temp4                         -1
        23 Temp3                         -1
        24 Local Time Step               -1
        25 Luminance                     -1
        26 Qrad[X]                       -1
        27 Qrad[Y]                       -1
        28 Qrad[Z]                       -1
        29 Srad                          -1
        30 ITSRI                         -1
        31 Absorp                        -1
        32 Emiss                         -1
        33 CoefAb                        -1
        34 TimeAverage1                  -1
        35 TimeAverage2                  -1
         --           --                --

 --- Log files
       NTLIST =             10 (Output frequency     )

       Number Name                IWARNI verbosity level
                                      (-999: not applicable)

         2 Velocity                       0
         5 Pressure                       0
         6 k                              0
         7 epsilon                        0
         8 mixture_fraction               1
         9 mixture_fraction               1
        10 Enthalpy                       1
        11 Density                     -999
        13 TurbVisc                    -999
        14 CourantNb                   -999
        15 FourierNb                   -999
        16 total_pressure              -999
        17 temperature                 -999
        18 ym_fuel                     -999
        19 ym_oxyd                     -999
        20 ym_prod                     -999
        21 KABS                        -999
        22 Temp4                       -999
        23 Temp3                       -999
        24 Local Time Step             -999
        34 TimeAverage1                -999
        35 TimeAverage2                -999
         --           --                --

 --- Additional post-processing variables (ipstdv)
       ipstfo =              1 (Force exerted by the
                                     fluid on the boundary)
       ipstft =              1 (Thermal flux   at boundary)
       ipstnu =              0 (Dimensionless thermal
                                          flux at boundary)


 -----------------------------------------------------------


 ** ALE METHOD (MOVING MESH)
    -----------

       IALE   =              0 (1: activated                )
       NALINF =              0 (Fluid initialization
                                                  iterations)
       IFLXMW =              0 (ALE mass flux computation
                                0: thanks to vertices
                                1: thanks to mesh velocity)


 -----------------------------------------------------------


Gestionnaires de sorties pour post traitement :
-----------------------------------------------

  -1: nom: results
      répertoire: postprocessing
      format: EnSight Gold
      options: binary
      dépendence en temps: maillage fixe
      sorties: tous les 10 pas de temps


 Lecture du fichier :  mesh_input
 Fin de la lecture :  mesh_input
 Pas de fichier "partition_input/domain_number_8" disponible;

 ----------------------------------------------------------

 Partitionnement par courbe fractale : Morton (dans boîte englobante).
  Nombre de cellules par domaine (histogramme) :
    [      45213 ;      45214 ] =          8

 Partitionnement terminé (2.95 s)


 ----------------------------------------------------------

 Construction du halo avec voisinage standard
 ============================================

 Création de l'interface entre faces
 Création de l'interface entre sommets
 Création des halos
 Définition des halos
    Définition du halo local
    Création du halo distant
    Mise à jour de la connectivité faces -> cellules

 Histogramme de la distribution du nombre de cellules sur les rangs :

    valeur minimale =            45213
    valeur maximale =            45214

      1 : [      45213 ;      45214 ] =          8

 ----------------------------------------------------------

 Histogramme de la distribution du nombre de cellules
 standard + halo sur les rangs :

    valeur minimale =            47073
    valeur maximale =            50927

      1 : [      47073 ;      47843 [ =          1
      2 : [      47843 ;      48614 [ =          1
      3 : [      48614 ;      49385 [ =          2
      4 : [      49385 ;      50156 [ =          3
      5 : [      50156 ;      50927 ] =          1

 ----------------------------------------------------------

 Histogramme de la distribution des cellules fantômes sur les rangs :

    valeur minimale =             1859
    valeur maximale =             5714

      1 : [       1859 ;       2630 [ =          1
      2 : [       2630 ;       3401 [ =          1
      3 : [       3401 ;       4172 [ =          2
      4 : [       4172 ;       4943 [ =          3
      5 : [       4943 ;       5714 ] =          1

 ----------------------------------------------------------

 Histogramme de la distribution du nombre de faces interne sur les rangs :

    valeur minimale =            89651
    valeur maximale =            92772

      1 : [      89651 ;      90275 [ =          1
      2 : [      90275 ;      90899 [ =          0
      3 : [      90899 ;      91523 [ =          1
      4 : [      91523 ;      92147 [ =          1
      5 : [      92147 ;      92772 ] =          5

 ----------------------------------------------------------

 Histogramme de la distribution du nombre de faces de bord sur les rangs :

    valeur minimale =             2491
    valeur maximale =             3941

      1 : [       2491 ;       2781 [ =          1
      2 : [       2781 ;       3071 [ =          2
      3 : [       3071 ;       3361 [ =          0
      4 : [       3361 ;       3651 [ =          3
      5 : [       3651 ;       3941 ] =          2

 ----------------------------------------------------------

 Histogramme du nombre de faces internes par cellule :

    valeur minimale =                2
    valeur maximale =                6

      1 : [          2 ;          3 [ =        418
      2 : [          3 ;          4 [ =      25468
      3 : [          4 ;          5 [ =     331902
      4 : [          5 ;          6 ] =       3920

 ----------------------------------------------------------

 Histogramme de la distribution du nombre de domaines voisins sur les rangs :

    valeur minimale =                5
    valeur maximale =                7

      1 : [          5 ;          6 [ =          1
      2 : [          6 ;          7 ] =          7

 ----------------------------------------------------------

 Définition globale du nombre d'éléments (cellules, sommets, faces...)
 Synchronisation des familles des cellules

 Renumérotation du maillage:

   renumérotation des cellules :
     pré-numérotation                              none
     cellules adjacentes aux halos en dernier :    non
     numérotation :                                none

   renumérotation des faces internes :
     pré-ordonnancement par cellules adjacentes :  id le plus bas en premier
     faces adjacentes aux halos en dernier :       non
     numérotation :                                none

   renumérotation des faces de bord :
     numérotation :                                none

 ----------------------------------------------------------

 Histogramme de la largeur de bande de matrice maillage volumique sur les rangs :

    valeur minimale =            46529
    valeur maximale =            50548

      1 : [      46529 ;      47332 [ =          1
      2 : [      47332 ;      48136 [ =          1
      3 : [      48136 ;      48940 [ =          1
      4 : [      48940 ;      49744 [ =          4
      5 : [      49744 ;      50548 ] =          1

 Histogramme du profil/lignes de matrice maillage volumique sur les rangs :

    valeur minimale =             7831
    valeur maximale =            12785

      1 : [       7831 ;       8821 [ =          2
      2 : [       8821 ;       9812 [ =          0
      3 : [       9812 ;      10803 [ =          2
      4 : [      10803 ;      11794 [ =          1
      5 : [      11794 ;      12785 ] =          3

 Coordonnées du maillage         minimale    et maximale
                       X :  3.4694470e-18  6.2740123e-01
                       Y : -2.3316136e-01  1.6683863e-01
                       Z : -1.8723129e-01  1.6800664e-01
 Maillage
     Nombre de cellules :       361708
     Nombre de faces internes : 713442
     Nombre de faces de bord :  26498
     Nombre de sommets :        69379

 Groupes :
    "ENTREEAIR"
       faces de bord :              107
    "ENTREEGAZ"
       faces de bord :               44
    "PART_1_1"
       faces intérieures :         1400
       faces de bord :             8838
    "PART_1_2"
       faces intérieures :          655
       faces de bord :              176
    "PART_1_3"
       faces intérieures :          655
       faces de bord :              184
    "PART_1_4"
       faces intérieures :          400
       faces de bord :            15982
    "PART_1_5"
       faces de bord :              927
    "SORTIEAIR"
       faces de bord :              240

 --- Information sur les volumes
       Volume de controle minimal =  2.7596738e-11
       Volume de controle maximal =  1.8763277e-06
       Volume total du domaine    =  6.2931571e-02

  Critère 1 : orthogonalité :
    Nombre de mauvaises cellules détecté : 0 -->   0 %

  Critère 2 : décentrement :
    Nombre de mauvaises cellules détecté : 8169 -->   2 %

  Critère 3 : qualité du gradient moindres-carrés :
    Nombre de mauvaises cellules détecté : 91 -->   0 %

  Critère 4 : ratio des volumes de cellules :
    Nombre de mauvaises cellules détecté : 0 -->   0 %

  Critère 5 : culpabilité par association :
    Nombre de mauvaises cellules détecté : 0 -->   0 %

 Attention :
 ---------
    Un défaut de qualité de maillage a été détecté

    Le maillage devrait être revu en fonction des critères indiqués.

    Le calcul sera effectué mais la qualité de la solution peut être dégradée...

 Computing geometric quantities (0.136 s)
                                                             
       ALMAX  =    0.39776E+00 (Characteristic length       )
       ALMAX is the cubic root of the domain volume.

       ALMAX is the length used to initialize the turbulence.

 -----------------------------------------------------------


 ** VARIABLES INITIALIZATION
    ------------------------

 -----------------------------------------
  Variable          Min. value  Max. value
 -----------------------------------------                   
  Velocity          0.0000E+00  0.0000E+00
  Velocity          0.0000E+00  0.0000E+00
  Velocity          0.0000E+00  0.0000E+00
  Pressure          0.0000E+00  0.0000E+00
  k                 0.6000E-03  0.6000E-03
  epsilon           0.3325E-05  0.3325E-05
  mixture_fraction  0.0000E+00  0.0000E+00
  mixture_fraction  0.0000E+00  0.0000E+00
  Enthalpy          0.1903E+04  0.1903E+04
 ---------------------------------

                dt  0.1000E-01  0.1000E-01
 ---------------------------------



-------------------------------------------------------------

 ------------------------------------------------------------- 



===============================================================



                       MAIN CALCULATION                      
                       ================                      


===============================================================




===============================================================


 INSTANT    0.100000000E-01   TIME STEP NUMBER               1
 ============================================================= 


  PPPDFR : Points de clipping de la variance :            0
  PPPDFR : Basculement sans les PDF            0

PDF RECTANGLE - PICS DE DIRAC COPDFR
MOYENNE, VARIANCE DU TRACEUR TRANPORTES
NOMBRE DE POINTS TURBULENTS (PASSAGE PAR LES PDF)   =      0
NOMBRE DE POINTS DE CALCULS                         = 361708
 Nb points avec PDF rectangle sans Dirac                =      0
 - - - - - - - - - -- - - -  et Dirac en FMINI          =      0
 - - - - - - - - - -- - - - - - - - - -  FMAXI          =      0
 - - - - - - - - - - - - - - - Diracs en FMINI et FMAXI =      0

 -----------------------------------------
 Property           Min. value  Max. value
 -----------------------------------------
  density           0.1189E+01  0.1189E+01
  molecular_viscos  0.1830E-04  0.1830E-04
  turbulent_viscos  0.1159E-01  0.1159E-01
 -----------------------------------------

 --- Diffusivity:
 -----------------------------------------------
 Scalar           Number  Min. value  Max. value
 -----------------------------------------------
 mixture_fraction      1  0.1830E-04  0.1830E-04
 mixture_fraction      2  0.1830E-04  0.1830E-04
 Enthalpy              3  0.1830E-04  0.1830E-04
 -----------------------------------------------


   ** INFORMATIONS SUR LE MODULE DE RAYONNEMENT
      ------------------------------------------
    Initialisation de la temperature de paroi
    (TPAROI) avec le profil utilisateur (TINTP)
    et du flux incident aux parois (QINCID).



   ** INFORMATION ON BOUNDARY FACES TYPE
      ----------------------------------

-------------------------------------------------------------------------
Boundary type          Code    Nb faces
-------------------------------------------------------------------------
Inlet                     2         151
Smooth wall               5       26107
Rough wall                6           0
Symmetry                  4           0
Free outlet               3         240
Free inlet               13           0
Convective inlet         15           0
Free surface             14           0
Undefined                 1           0
-------------------------------------------------------------------------


Boundary faces with free inlet/outlet detected
Update of reference point for total pressure
 XYZP0 =    0.52705E+00  -0.23316E+00  -0.13540E-01



   ** BOUNDARY MASS FLOW INFORMATION
      ------------------------------

---------------------------------------------------------------
Boundary type          Code    Nb faces           Mass flow
---------------------------------------------------------------
Inlet                     2         151         0.000000000E+00
Smooth wall               5       26107         0.000000000E+00
Rough wall                6           0         0.000000000E+00
Symmetry                  4           0         0.000000000E+00
Free outlet               3         240         0.000000000E+00
Free inlet               13           0         0.000000000E+00
Convective inlet         15           0         0.000000000E+00
Free surface             14           0         0.000000000E+00
Undefined                 1           0         0.000000000E+00
---------------------------------------------------------------


   ** BOUNDARY CONDITIONS FOR SMOOTH WALLS
   ---------------------------------------
------------------------------------------------------------
                                         Minimum     Maximum
------------------------------------------------------------
   Rel velocity at the wall uiptn :  0.00000E+00 0.00000E+00
   Friction velocity        uet   :  0.00000E+00 0.00000E+00
   Friction velocity        uk    :  0.13925E-02 0.10035E-01
   Dimensionless distance   yplus :  0.67737E-02 0.36936E+01
   Friction thermal sca.    tstar :  0.00000E+00 0.00000E+00
   Rough dim-less th. sca.  tplus :  0.67737E-02 0.36936E+01
   ------------------------------------------------------
   Nb of reversal of the velocity at the wall   :          0
   Nb of faces within the viscous sub-layer     :      25849
   Total number of wall faces                   :      26107
------------------------------------------------------------


@                                                            
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@                                                            
@ @@ WARNING: MESH TOO REFINED AT THE WALL                   
@    ========                                                
@    The mesh seems to be too refined at the wall to use     
@      a high-Reynolds turbulence model.                     
@                                                            
@    The last time step at which too small values for the    
@      dimensionless distance to the wall (yplus) have been  
@      observed is the time step          1
@                                                            
@    The minimum value for yplus must be greater than the    
@      limit value YPLULI =    0.23810E+01
@                                                            
@    Have a look at the distribution of yplus at the wall    
@      (with EnSight or ParaView for example) to conclude on 
@      the way the results quality might be affected.        
@                                                            
@    This warning is only printed at the first two           
@      occurences of the problem and at the last time step   
@      of the calculation. The vanishing of the message does 
@      not necessarily mean the vanishing of the problem.    
@                                                            
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@                                                            


   ** INFORMATIONS SUR LE TERME SOURCE RADIATIF
      -----------------------------------------
-------------------------------------------------------------------
Zone         Flux net radiatif (Watt) (normale unitaire sortante)
     1              0.1806E+02
     2              0.9468E-02
     3              0.1145E+02
-------------------------------------------------------------------
Flux net radiatif sur toutes les frontieres  Fnet =  0.2952E+02 Watt
Integrale volumique du terme source radiatif Srad = -0.2952E+02 Watt
(Si IDIVER = 1 ou 2 alors on doit avoir Srad = -Fnet)
-------------------------------------------------------------------

   ** SOLVING VARIABLE mixture_fraction
      ----------------

@
@ @@ WARNING: mixture_ CONVECTION-DIFFUSION-SOURCE TERMS
@    ========
@  Maximum number of iterations          1 reached
@

   ** SOLVING VARIABLE mixture_fraction
      ----------------

@
@ @@ WARNING: mixture_ CONVECTION-DIFFUSION-SOURCE TERMS
@    ========
@  Maximum number of iterations          1 reached
@

   ** SOLVING VARIABLE enthalpy        
      ----------------

@
@ @@ WARNING: enthalpy CONVECTION-DIFFUSION-SOURCE TERMS
@    ========
@  Maximum number of iterations          1 reached
@

 INSTANT    0.200000000E-01   TIME STEP NUMBER               2
 ============================================================= 


  PPPDFR : Points de clipping de la variance :            0
  PPPDFR : Basculement sans les PDF          222

PDF RECTANGLE - PICS DE DIRAC COPDFR
MOYENNE, VARIANCE DU TRACEUR TRANPORTES
NOMBRE DE POINTS TURBULENTS (PASSAGE PAR LES PDF)   = 173017
NOMBRE DE POINTS DE CALCULS                         = 361708
 Nb points avec PDF rectangle sans Dirac                =      0
 - - - - - - - - - -- - - -  et Dirac en FMINI          = 171083
 - - - - - - - - - -- - - - - - - - - -  FMAXI          =      0
 - - - - - - - - - - - - - - - Diracs en FMINI et FMAXI =   1934

  Clipping de HSTOE EN MIN EN   144570 POINTS
      Valeur Max :  -0.3599421E+11
      Valeur De Clipping :   -2959254.    

 --- Diffusivity:
 -----------------------------------------------
 Scalar           Number  Min. value  Max. value
 -----------------------------------------------
 mixture_fraction      1  0.1830E-04  0.1830E-04
 mixture_fraction      2  0.1830E-04  0.1830E-04
 Enthalpy              3  0.1830E-04  0.1830E-04
 -----------------------------------------------


   ** Information on wall temperature
      -------------------------------
-----------------------------------------------------------------------
Fixed profiles   Temp max (C)   Temp min (C)   Temp mean (C)  Net flux (W)
         1         0.2685E+02     0.2685E+02     0.2484E+02    -0.3278E+13
-----------------------------------------------------------------------


   ** BOUNDARY MASS FLOW INFORMATION
      ------------------------------

---------------------------------------------------------------
Boundary type          Code    Nb faces           Mass flow
---------------------------------------------------------------
Inlet                     2         151        -0.130900000E+00
Smooth wall               5       26107         0.000000000E+00
Rough wall                6           0         0.000000000E+00
Symmetry                  4           0         0.000000000E+00
Free outlet               3         240         0.130899986E+00
Free inlet               13           0         0.000000000E+00
Convective inlet         15           0         0.000000000E+00
Free surface             14           0         0.000000000E+00
Undefined                 1           0         0.000000000E+00
---------------------------------------------------------------


   ** BOUNDARY CONDITIONS FOR SMOOTH WALLS
   ---------------------------------------
------------------------------------------------------------
                                         Minimum     Maximum
------------------------------------------------------------
   Rel velocity at the wall uiptn : -0.35683E+03 0.80291E+02
   Friction velocity        uet   :  0.88031E-03 0.89040E+01
   Friction velocity        uk    :  0.12507E-02 0.75229E+01
   Dimensionless distance   yplus :  0.29537E+00 0.18514E+03
   Friction thermal sca.    tstar : -0.85393E+01 0.46700E+05
   Rough dim-less th. sca.  tplus :  0.29537E+00 0.17538E+02
   ------------------------------------------------------
   Nb of reversal of the velocity at the wall   :      22929
   Nb of faces within the viscous sub-layer     :        352
   Total number of wall faces                   :      26107
------------------------------------------------------------


@                                                            
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@                                                            
@ @@ WARNING: MESH TOO REFINED AT THE WALL                   
@    ========                                                
@    The mesh seems to be too refined at the wall to use     
@      a high-Reynolds turbulence model.                     
@                                                            
@    The last time step at which too small values for the    
@      dimensionless distance to the wall (yplus) have been  
@      observed is the time step          2
@                                                            
@    The minimum value for yplus must be greater than the    
@      limit value YPLULI =    0.23810E+01
@                                                            
@    Have a look at the distribution of yplus at the wall    
@      (with EnSight or ParaView for example) to conclude on 
@      the way the results quality might be affected.        
@                                                            
@    This warning is only printed at the first two           
@      occurences of the problem and at the last time step   
@      of the calculation. The vanishing of the message does 
@      not necessarily mean the vanishing of the problem.    
@                                                            
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@                                                            


   ** INFORMATIONS SUR LE TERME SOURCE RADIATIF
      -----------------------------------------
-------------------------------------------------------------------
Zone         Flux net radiatif (Watt) (normale unitaire sortante)
     1              0.2153E+02
     2              0.1189E-01
     3              0.1114E+02
-------------------------------------------------------------------
Flux net radiatif sur toutes les frontieres  Fnet =  0.3268E+02 Watt
Integrale volumique du terme source radiatif Srad = -0.3268E+02 Watt
(Si IDIVER = 1 ou 2 alors on doit avoir Srad = -Fnet)
-------------------------------------------------------------------

   ** SOLVING VARIABLE mixture_fraction
      ----------------

@
@ @@ WARNING: mixture_ CONVECTION-DIFFUSION-SOURCE TERMS
@    ========
@  Maximum number of iterations          1 reached
@

   ** SOLVING VARIABLE mixture_fraction
      ----------------

@
@ @@ WARNING: mixture_ CONVECTION-DIFFUSION-SOURCE TERMS
@    ========
@  Maximum number of iterations          1 reached
@

   ** SOLVING VARIABLE enthalpy        
      ----------------

@
@ @@ WARNING: enthalpy CONVECTION-DIFFUSION-SOURCE TERMS
@    ========
@  Maximum number of iterations          1 reached
@

 INSTANT    0.300000000E-01   TIME STEP NUMBER               3
 ============================================================= 


  PPPDFR : Points de clipping de la variance :            0
  PPPDFR : Basculement sans les PDF          166

PDF RECTANGLE - PICS DE DIRAC COPDFR
MOYENNE, VARIANCE DU TRACEUR TRANPORTES
NOMBRE DE POINTS TURBULENTS (PASSAGE PAR LES PDF)   = 110609
NOMBRE DE POINTS DE CALCULS                         = 361708
 Nb points avec PDF rectangle sans Dirac                =  28921
 - - - - - - - - - -- - - -  et Dirac en FMINI          =  81346
 - - - - - - - - - -- - - - - - - - - -  FMAXI          =      7
 - - - - - - - - - - - - - - - Diracs en FMINI et FMAXI =    335

  Clipping de HSTOE EN MAX EN     3031 POINTS
      Valeur Max :    950528.6    
      Valeur De Clipping :   -187205.9    

  Clipping de HSTOE EN MIN EN    81842 POINTS
      Valeur Max :  -0.1011810E+12
      Valeur De Clipping :   -2959254.    

 --- Diffusivity:
 -----------------------------------------------
 Scalar           Number  Min. value  Max. value
 -----------------------------------------------
 mixture_fraction      1  0.1830E-04  0.1830E-04
 mixture_fraction      2  0.1830E-04  0.1830E-04
 Enthalpy              3  0.1830E-04  0.1830E-04
 -----------------------------------------------


   ** Information on wall temperature
      -------------------------------
-----------------------------------------------------------------------
Fixed profiles   Temp max (C)   Temp min (C)   Temp mean (C)  Net flux (W)
         1         0.2685E+02     0.2685E+02     0.2484E+02    -0.3186E+13
-----------------------------------------------------------------------
