command: 
 ./cs_solver --mpi

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

                                  (R)
                      Code_Saturne

                      Version 4.0.0


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

  revision 4.0.0
  build Mon 18 May 2015 06:02:44 PM CST
  MPI version 2.1 (Open MPI 1.5.0)


  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.

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


Local case configuration:

  Date:                Tue 19 May 2015 06:50:11 PM CST
  System:              Linux 2.6.18-164.el5
  Machine:             cn032
  Processor:           model name	: Intel(R) Xeon(R) CPU           X5670  @ 2.93GHz
  Memory:              20058 MB
  User:                ysh ()
  Directory:           /store/ysh/hy4.0/CASE1/RESU/20150519-1850
  MPI ranks:           12 (appnum attribute: 0)
  I/O read method:     collective MPI-IO (explicit offsets)
  I/O write method:    collective MPI-IO (explicit offsets)
  I/O rank step:        1

  External libraries for partitioning:
    METIS

 Reading metadata from file: "mesh_input"

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

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


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



 Fractions massiques elementaires / FOL  
  C =   0.846700000000000     
  H =   0.137900000000000     
  O =   5.000000000000000E-003
  S =   5.000000000000000E-004
  In=   9.899999999999996E-003
  Fraction massique de coke / FOL  0.200000000000000     
 Fractions massiques elementaires / coke 
  C =   0.950500000000000     
  H =   0.000000000000000E+000
  O =   0.000000000000000E+000
  S =   0.000000000000000E+000
  In=   4.949999999999997E-002
  nHCFOV =    2.53415792295321        172.335937500000     
   68.0052083333333     
  mtoFOV =    1.00000000000000     
 Fractions massiques sp�cifiques / FOV 
  H2S =   6.640625000000002E-004
  CO  =   1.093750000000000E-002
  CHn =   0.988398437500000     
  ..n =    2.53415792295321     
 Le FOV est equivalent a un melange 
  H2S =   6.640625000000002E-004
  CO  =   1.093750000000000E-002
  CH4 =   0.290604166666667     
 C2H4 =   0.697794270833333     
  nHCFOV 2 =    2.53415792295321     
  Wmole IFOV 1ere formule=   2.294120313704157E-002  0.290604166666667     
  0.697794270833333     
  Wmole IFOV 2eme formule=   1.453415792295321E-002  0.290604166666667     
  0.697794270833333     
 Donnees thermo pour le fuel
 PCIFOL    40467310.0000000     
 H02FOL   -3972045.73656029     
 CP2FOL    1800.00000000000     
 HRFVAP    3929787.35278864     
 H02FOV   -52822.9797145557     
  Verification des enthalpies de formation
  CH4    -4680184.90849523     
  C2H4    1873415.00225193     
  FOV    -52822.9797145557     
  FOL    -3972315.73656029     
  KF     -269.999999999959     
  Classe D =            1  1.000000000000000E-005  5.848035476425733E-006
  2.147229169018941E-006
  Classe D =            2  2.500000000000000E-005  1.462008869106433E-005
  5.368072922547352E-006
  Classe D =            3  5.000000000000000E-005  2.924017738212866E-005
  1.073614584509470E-005
                                                             
 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:
    ----------------

 --- Fuel:            Local Mean Combustible Model
       OPTION =              0


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

 --- Physics
       NVAR   =             22 (Nb variables                )
       NSCAL  =             16 (Nb scalars                  )
       NSCAUS =              0 (Nb user scalars             )
       NSCAPP =             16 (Nb specific physics scalars )
       NPROCE =             45 (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.43941E+00 (Reference density           )
       VISCL0 =    0.37000E-04 (Ref. molecular dyn. visc.   )
       CP0    =    0.10990E+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.80000E+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.42500E-04 (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 =              0 (0: none                     )
                               (1: temperature in Kelvin    )
                               (2: temperature in Celsius   )
       ISCALT =              1 (Thermal scalar number       )


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

 --- Continuous phase:

   - Commons
       ITURB  =             20 (Turbulence model            )
       IDEUCH =              1 (0: one-scale model          )
                               (1: two-scale model          )
                               (2: invariant wall function  )
       IWALLT =              0 (Exch. coeff. correlation    )
                               (0: not activated            )
                               (1: activated                )
       YPLULI =    0.23810E+01 (Limit Y+                    )
       ILOGPO =              1 (0: power law (forbidden for
                                              k-epsilon)    )
                               (1: one-scale log law        )
       IGRHOK =              0 (1: computed Grad(rho k)     )

   - k-epsilon           (ITURB = 20)
       ALMAX  =    0.15000E+00 (Characteristic length       )
       UREF   =    0.10000E+03 (Characteristic velocity     )
       ICLKEP =              0 (k-epsilon clipping model    )
       IKECOU =              1 (k-epsilon coupling mode     )
       IGRAKE =              1 (Account for gravity         )


   - 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))

   - k-epsilon           (ITURB = 20)
       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            )


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


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

    UNSTEADY ALGORITHM

 --- Time step parameters

       IDTVAR =              0 (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-05 (Minimum time step           )
       DTMAX  =    0.10000E-01 (Maximum time step           )
       DTREF  =    0.10000E-04 (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
 Turb Kinetic Ene      1  0.1000E+01
 Turb Dissipation      1  0.1000E+01
 Enthalpy              1  0.1000E+01
 NG_FOL01              1  0.1000E+01
 NG_FOL02              1  0.1000E+01
 NG_FOL03              1  0.1000E+01
 YFOL_FOL01            1  0.1000E+01
 YFOL_FOL02            1  0.1000E+01
 YFOL_FOL03            1  0.1000E+01
 HLF_FOL01             1  0.1000E+01
 HLF_FOL02             1  0.1000E+01
 HLF_FOL03             1  0.1000E+01
 Fr_VAP                1  0.1000E+01
 Fr_HET                1  0.1000E+01
 Var_CB                1  0.1000E+01
 FR_HCN                1  0.1000E+01
 FR_NO                 1  0.1000E+01
 Enth_Ox               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.80000E+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.10E+01 0.10E+01
 Pressure              0      1      1      1      0 0.00E+00 0.10E+01
 Turb Kinetic Ene      1      1      1      1      0 0.00E+00 0.10E+01
 Turb Dissipation      1      1      1      1      0 0.00E+00 0.10E+01
 Enthalpy              1      1      1      1      0 0.00E+00 0.10E+01
 NG_FOL01              1      1      1      1      0 0.00E+00 0.10E+01
 NG_FOL02              1      1      1      1      0 0.00E+00 0.10E+01
 NG_FOL03              1      1      1      1      0 0.00E+00 0.10E+01
 YFOL_FOL01            1      1      1      1      0 0.00E+00 0.10E+01
 YFOL_FOL02            1      1      1      1      0 0.00E+00 0.10E+01
 YFOL_FOL03            1      1      1      1      0 0.00E+00 0.10E+01
 HLF_FOL01             1      1      1      1      0 0.00E+00 0.10E+01
 HLF_FOL02             1      1      1      1      0 0.00E+00 0.10E+01
 HLF_FOL03             1      1      1      1      0 0.00E+00 0.10E+01
 Fr_VAP                1      1      1      1      0 0.00E+00 0.10E+01
 Fr_HET                1      1      1      1      0 0.00E+00 0.10E+01
 Var_CB                1      1      1      1      0 0.00E+00 0.10E+01
 FR_HCN                1      1      1      1      0 0.00E+00 0.10E+01
 FR_NO                 1      1      1      1      0 0.00E+00 0.10E+01
 Enth_Ox               1      1      1      1      0 0.00E+00 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-06  0.0000E+00
 Pressure            100      2  0.1000E-04  0.1000E-06  0.0000E+00
 Turb Kinetic Ene    100      1  0.1000E-04  0.1000E-06  0.0000E+00
 Turb Dissipation    100      1  0.1000E-04  0.1000E-06  0.0000E+00
 Enthalpy            100      1  0.1000E-04  0.1000E-06  0.0000E+00
 NG_FOL01            100      1  0.1000E-04  0.1000E-06  0.0000E+00
 NG_FOL02            100      1  0.1000E-04  0.1000E-06  0.0000E+00
 NG_FOL03            100      1  0.1000E-04  0.1000E-06  0.0000E+00
 YFOL_FOL01          100      1  0.1000E-04  0.1000E-06  0.0000E+00
 YFOL_FOL02          100      1  0.1000E-04  0.1000E-06  0.0000E+00
 YFOL_FOL03          100      1  0.1000E-04  0.1000E-06  0.0000E+00
 HLF_FOL01           100      1  0.1000E-04  0.1000E-06  0.0000E+00
 HLF_FOL02           100      1  0.1000E-04  0.1000E-06  0.0000E+00
 HLF_FOL03           100      1  0.1000E-04  0.1000E-06  0.0000E+00
 Fr_VAP              100      1  0.1000E-04  0.1000E-06  0.0000E+00
 Fr_HET              100      1  0.1000E-04  0.1000E-06  0.0000E+00
 Var_CB              100      1  0.1000E-04  0.1000E-06  0.0000E+00
 FR_HCN              100      1  0.1000E-04  0.1000E-06  0.0000E+00
 FR_NO               100      1  0.1000E-04  0.1000E-06  0.0000E+00
 Enth_Ox             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
 Turb Kinetic Ene      1     -1  0.1500E+01
 Turb Dissipation      1     -1  0.1500E+01
 Enthalpy              0     -1  0.1500E+01
 NG_FOL01              0     -1  0.1500E+01
 NG_FOL02              0     -1  0.1500E+01
 NG_FOL03              0     -1  0.1500E+01
 YFOL_FOL01            0     -1  0.1500E+01
 YFOL_FOL02            0     -1  0.1500E+01
 YFOL_FOL03            0     -1  0.1500E+01
 HLF_FOL01             0     -1  0.1500E+01
 HLF_FOL02             0     -1  0.1500E+01
 HLF_FOL03             0     -1  0.1500E+01
 Fr_VAP                0     -1  0.1500E+01
 Fr_HET                0     -1  0.1500E+01
 Var_CB                0     -1  0.1500E+01
 FR_HCN                0     -1  0.1500E+01
 FR_NO                 0     -1  0.1500E+01
 Enth_Ox               0     -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-07      1
 Pressure          0.1000E-07      1
 Turb Kinetic Ene  0.1000E-07      1
 Turb Dissipation  0.1000E-07      1
 Enthalpy          0.1000E-07      1
 NG_FOL01          0.1000E-07      1
 NG_FOL02          0.1000E-07      1
 NG_FOL03          0.1000E-07      1
 YFOL_FOL01        0.1000E-07      1
 YFOL_FOL02        0.1000E-07      1
 YFOL_FOL03        0.1000E-07      1
 HLF_FOL01         0.1000E-07      1
 HLF_FOL02         0.1000E-07      1
 HLF_FOL03         0.1000E-07      1
 Fr_VAP            0.1000E-07      1
 Fr_HET            0.1000E-07      1
 Var_CB            0.1000E-07      1
 FR_HCN            0.1000E-07      1
 FR_NO             0.1000E-07      1
 Enth_Ox           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
-------------------------------------------------------------
 Enthalpy              1      0      0  0.4250E-04  0.7000E+00
 NG_FOL01              2      0      0  0.3700E-04  0.7000E+00
 NG_FOL02              3      0      0  0.3700E-04  0.7000E+00
 NG_FOL03              4      0      0  0.3700E-04  0.7000E+00
 YFOL_FOL01            5      0      0  0.3700E-04  0.7000E+00
 YFOL_FOL02            6      0      0  0.3700E-04  0.7000E+00
 YFOL_FOL03            7      0      0  0.3700E-04  0.7000E+00
 HLF_FOL01             8      0      0  0.3700E-04  0.7000E+00
 HLF_FOL02             9      0      0  0.3700E-04  0.7000E+00
 HLF_FOL03            10      0      0  0.3700E-04  0.7000E+00
 Fr_VAP               11      0      0  0.3700E-04  0.7000E+00
 Fr_HET               12      0      0  0.3700E-04  0.7000E+00
 Var_CB               13      0      0  0.3700E-04  0.7000E+00
 FR_HCN               14      0      0  0.3700E-04  0.7000E+00
 FR_NO                15      0      0  0.3700E-04  0.7000E+00
 Enth_Ox              16      0      0  0.3700E-04  0.7000E+00
---------------------------------------------------------------------

------------------------------------
 Variable         Number      RVARFL
------------------------------------
 Enthalpy              1  0.8000E+00
 NG_FOL01              2  0.8000E+00
 NG_FOL02              3  0.8000E+00
 NG_FOL03              4  0.8000E+00
 YFOL_FOL01            5  0.8000E+00
 YFOL_FOL02            6  0.8000E+00
 YFOL_FOL03            7  0.8000E+00
 HLF_FOL01             8  0.8000E+00
 HLF_FOL02             9  0.8000E+00
 HLF_FOL03            10  0.8000E+00
 Fr_VAP               11  0.8000E+00
 Fr_HET               12  0.8000E+00
 Var_CB               13  0.8000E+00
 FR_HCN               14  0.8000E+00
 FR_NO                15  0.8000E+00
 Enth_Ox              16  0.8000E+00
-------------------------------------------

-------------------------------------------------------
 Variable         Number ICLVFL      SCAMIN      SCAMAX
-------------------------------------------------------
 Enthalpy              1     -1 -0.1000E+13  0.1000E+13
 NG_FOL01              2     -1 -0.3972E+07  0.1000E+13
 NG_FOL02              3     -1 -0.3972E+07  0.1000E+13
 NG_FOL03              4     -1 -0.3972E+07  0.1000E+13
 YFOL_FOL01            5     -1  0.0000E+00  0.4000E+00
 YFOL_FOL02            6     -1  0.0000E+00  0.4000E+00
 YFOL_FOL03            7     -1  0.0000E+00  0.4000E+00
 HLF_FOL01             8     -1 -0.1000E+13  0.1000E+13
 HLF_FOL02             9     -1 -0.1000E+13  0.1000E+13
 HLF_FOL03            10     -1 -0.1000E+13  0.1000E+13
 Fr_VAP               11     -1  0.0000E+00  0.1000E+01
 Fr_HET               12     -1  0.0000E+00  0.1000E+01
 Var_CB               13     -1  0.0000E+00  0.2500E+00
 FR_HCN               14     -1  0.0000E+00  0.1000E+01
 FR_NO                15     -1  0.0000E+00  0.1000E+01
 Enth_Ox              16     -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
          4  0.0000E+00         0  0.0000E+00
          5  0.0000E+00         0  0.0000E+00
          6  0.0000E+00         0  0.0000E+00
          7  0.0000E+00         0  0.0000E+00
          8  0.0000E+00         0  0.0000E+00
          9  0.0000E+00         0  0.0000E+00
         10  0.0000E+00         0  0.0000E+00
         11  0.0000E+00         0  0.0000E+00
         12  0.0000E+00         0  0.0000E+00
         13  0.0000E+00         0  0.0000E+00
         14  0.0000E+00         0  0.0000E+00
         15  0.0000E+00         0  0.0000E+00
         16  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 =          10000 (Final time step required    )

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


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


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

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

 --- Post-processed variables


Velocity
Pressure
Turb Kinetic Energy
Turb Dissipation
Enthalpy
NG_FOL01
NG_FOL02
NG_FOL03
YFOL_FOL01
YFOL_FOL02
YFOL_FOL03
HLF_FOL01
HLF_FOL02
HLF_FOL03
Fr_VAP
Fr_HET
Var_CB
FR_HCN
FR_NO
Enth_Ox
Density
Boundary Density
Laminar Viscosity
Turb Viscosity
CFL
Fourier Number
Total Pressure
T_Gas
Rho_Gas
Ym_FO0
Ym_FOV
Ym_CO
Ym_H2S
Ym_H2
Ym_HCN
Ym_NH3
Ym_O2
Ym_CO2
Ym_H2O
Ym_SO2
Ym_N2
EXP1
EXP1
EXP3
T_Fuel_01
T_Fuel_02
T_Fuel_03
Rho_Fuel_01
Rho_Fuel_02
Rho_Fuel_03
Diam_Drop_01
Diam_Drop_02
Diam_Drop_03
H1-Hlf_01
H1-Hlf_02
H1-Hlf_03
EVA_Fuel_01
EVA_Fuel_02
EVA_Fuel_03
Het_TS_Fuel_01
Het_TS_Fuel_02
Het_TS_Fuel_03
Balance_C
Balance_O
Balance_H
         --           --

 --- 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 Turb Kinetic Ene              -1
         7 Turb Dissipation              -1
         8 Enthalpy                      -1
         9 NG_FOL01                      -1
        10 NG_FOL02                      -1
        11 NG_FOL03                      -1
        12 YFOL_FOL01                    -1
        13 YFOL_FOL02                    -1
        14 YFOL_FOL03                    -1
        15 HLF_FOL01                     -1
        16 HLF_FOL02                     -1
        17 HLF_FOL03                     -1
        18 Fr_VAP                        -1
        19 Fr_HET                        -1
        20 Var_CB                        -1
        21 FR_HCN                        -1
        22 FR_NO                         -1
        23 Enth_Ox                       -1
        26 Turb Viscosity                -1
         --           --                --

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

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

         2 Velocity                       2
         5 Pressure                       2
         6 Turb Kinetic Ene               1
         7 Turb Dissipation               1
         8 Enthalpy                       1
         9 NG_FOL01                       1
        10 NG_FOL02                       1
        11 NG_FOL03                       1
        12 YFOL_FOL01                     1
        13 YFOL_FOL02                     1
        14 YFOL_FOL03                     1
        15 HLF_FOL01                      1
        16 HLF_FOL02                      1
        17 HLF_FOL03                      1
        18 Fr_VAP                         1
        19 Fr_HET                         1
        20 Var_CB                         1
        21 FR_HCN                         1
        22 FR_NO                          1
        23 Enth_Ox                        1
        24 Density                     -999
        25 Laminar Viscosit            -999
        26 Turb Viscosity              -999
        27 CFL                         -999
        28 Fourier Number              -999
        29 Total Pressure              -999
        30 T_Gas                       -999
        31 Rho_Gas                     -999
        32 Ym_FO0                      -999
        33 Ym_FOV                      -999
        34 Ym_CO                       -999
        35 Ym_H2S                      -999
        36 Ym_H2                       -999
        37 Ym_HCN                      -999
        38 Ym_NH3                      -999
        39 Ym_O2                       -999
        40 Ym_CO2                      -999
        41 Ym_H2O                      -999
        42 Ym_SO2                      -999
        43 Ym_N2                       -999
        45 EXP1                        -999
        46 EXP1                        -999
        47 EXP3                        -999
        48 T_Fuel_01                   -999
        49 T_Fuel_02                   -999
        50 T_Fuel_03                   -999
        51 Rho_Fuel_01                 -999
        52 Rho_Fuel_02                 -999
        53 Rho_Fuel_03                 -999
        54 Diam_Drop_01                -999
        55 Diam_Drop_02                -999
        56 Diam_Drop_03                -999
        57 H1-Hlf_01                   -999
        58 H1-Hlf_02                   -999
        59 H1-Hlf_03                   -999
        60 EVA_Fuel_01                 -999
        61 EVA_Fuel_02                 -999
        62 EVA_Fuel_03                 -999
        63 Het_TS_Fuel_01              -999
        64 Het_TS_Fuel_02              -999
        65 Het_TS_Fuel_03              -999
        66 Balance_C                   -999
        67 Balance_O                   -999
        68 Balance_H                   -999
         --           --                --

 --- Additional post-processing variables (ipstdv)
       ipstfo =              0 (Force exerted by the
                                     fluid on the boundary)
       ipstft =              0 (Thermal flux   at boundary)
       ipsttb =              0 (Temperature    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)


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


Postprocessing output writers:
------------------------------

  -1: name: results
      directory: postprocessing
      format: EnSight Gold
      options: 
      time dependency: fixed mesh
      output: at calculation end


 Reading file:        mesh_input
 Finished reading:    mesh_input
 No "partition_input/domain_number_12" file available;

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

 Partitioning by space-filling curve: Morton (in bounding box).
  Number of cells per domain (histogramm):
    [      11758 ;      11759 ] =         12

 Partitioning finished (0.764 s)


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

 Halo construction with standard neighborhood
 ============================================

 Face interfaces creation
 Vertex interfaces creation
 Halo creation
 Halo definition
    Local halo definition
    Distant halo creation
    Updating the faces -> cells connectivity

 Histogram of the number of cells per rank:

    minimum value =              11758
    maximum value =              11759

      1 : [      11758 ;      11759 ] =         12

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

 Histogram of the number of standard + halo cells per rank:

    minimum value =              13283
    maximum value =              14474

      1 : [      13283 ;      13521 [ =          1
      2 : [      13521 ;      13759 [ =          2
      3 : [      13759 ;      13997 [ =          3
      4 : [      13997 ;      14235 [ =          1
      5 : [      14235 ;      14474 ] =          5

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

 Histogram of the number of ghost cells per rank:

    minimum value =               1525
    maximum value =               2716

      1 : [       1525 ;       1763 [ =          1
      2 : [       1763 ;       2001 [ =          2
      3 : [       2001 ;       2239 [ =          3
      4 : [       2239 ;       2477 [ =          1
      5 : [       2477 ;       2716 ] =          5

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

 Histogram of the number of interior faces per rank:

    minimum value =              29673
    maximum value =              30785

      1 : [      29673 ;      29895 [ =          3
      2 : [      29895 ;      30117 [ =          2
      3 : [      30117 ;      30340 [ =          2
      4 : [      30340 ;      30562 [ =          2
      5 : [      30562 ;      30785 ] =          3

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

 Histogram of the number of boundary faces per rank:

    minimum value =                395
    maximum value =               1575

      1 : [        395 ;        631 [ =          3
      2 : [        631 ;        867 [ =          4
      3 : [        867 ;       1103 [ =          1
      4 : [       1103 ;       1339 [ =          3
      5 : [       1339 ;       1575 ] =          1

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

 Histogram of the number of interior faces per cell:

    minimum value =                  3
    maximum value =                  5

      1 : [          3 ;          4 [ =        200
      2 : [          4 ;          5 ] =     140900

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

 Histogram of the number of neighboring domains per rank:

    minimum value =                  5
    maximum value =                 11

      1 : [          5 ;          6 [ =          2
      2 : [          6 ;          7 [ =          1
      3 : [          7 ;          8 [ =          1
      4 : [          8 ;          9 [ =          3
      5 : [          9 ;         11 ] =          5

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

 Global definition of the number of elements (cells, vertices, faces...)
Synchronizing cell families

 Renumbering mesh:

   renumbering for cells:
     pre-numbering:                       none
     cells adjacent to ghost cells last:  no
     numbering:                           none

   renumbering for interior faces:
     cell adjacency pre-ordering:         lowest id first
     faces adjacent to ghost cells last:  no
     numbering:                           none

   renumbering for boundary faces:
     numbering:                           none

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

 Histogram of volume mesh matrix bandwidth per rank:

    minimum value =              12393
    maximum value =              14073

      1 : [      12393 ;      12729 [ =          2
      2 : [      12729 ;      13065 [ =          1
      3 : [      13065 ;      13401 [ =          2
      4 : [      13401 ;      13737 [ =          5
      5 : [      13737 ;      14073 ] =          2

 Histogram of volume mesh matrix profile/lines per rank:

    minimum value =                994
    maximum value =               2027

      1 : [        994 ;       1200 [ =          2
      2 : [       1200 ;       1407 [ =          0
      3 : [       1407 ;       1613 [ =          3
      4 : [       1613 ;       1820 [ =          5
      5 : [       1820 ;       2027 ] =          2

 Mesh coordinates:               minimum    and maximum
                       X : -1.0000001e-02  1.0000001e-02
                       Y : -1.0000001e-02  1.0000001e-02
                       Z :  0.0000000e+00  1.0000000e-01
 Mesh
     Number of cells:          141100
     Number of interior faces: 347428
     Number of boundary faces: 10644
     Number of vertices:       74562

 Groups:
    "airin"
       boundary faces:         2806
    "fluid"
       cells:                141100
    "fuelin"
       boundary faces:           16
    "out"
       boundary faces:         2822
    "wall"
       boundary faces:         5000

 --- Information on the volumes
       Minimum control volume      =  7.9293499e-12
       Maximum control volume      =  1.7309725e-09
       Total volume for the domain =  3.1395263e-05

  Criterion 1: Orthogonality:
    Number of bad cells detected: 0 -->   0 %

  Criterion 2: Offset:
    Number of bad cells detected: 0 -->   0 %

  Criterion 3: Least-Squares Gradient Quality:
    Number of bad cells detected: 0 -->   0 %

  Criterion 4: Cells Volume Ratio:
    Number of bad cells detected: 0 -->   0 %

  Criterion 5: Guilt by Association:
    Number of bad cells detected: 0 -->   0 %

 Computing geometric quantities (0.025 s)
                                                             

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


 ** 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
  Turb Kinetic Ene  0.1000E-09  0.1000E-09
  Turb Dissipation  0.1000E-09  0.1000E-09
  Enthalpy          0.5275E+06  0.5275E+06
  NG_FOL01          0.0000E+00  0.0000E+00
  NG_FOL02          0.0000E+00  0.0000E+00
  NG_FOL03          0.0000E+00  0.0000E+00
  YFOL_FOL01        0.0000E+00  0.0000E+00
  YFOL_FOL02        0.0000E+00  0.0000E+00
  YFOL_FOL03        0.0000E+00  0.0000E+00
  HLF_FOL01         0.0000E+00  0.0000E+00
  HLF_FOL02         0.0000E+00  0.0000E+00
  HLF_FOL03         0.0000E+00  0.0000E+00
  Fr_VAP            0.0000E+00  0.0000E+00
  Fr_HET            0.0000E+00  0.0000E+00
  Var_CB            0.0000E+00  0.0000E+00
  FR_HCN            0.0000E+00  0.0000E+00
  FR_NO             0.0000E+00  0.0000E+00
  Enth_Ox           0.5275E+06  0.5275E+06
 ---------------------------------

                dt  0.1000E-04  0.1000E-04
 ---------------------------------



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

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



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



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


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




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


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



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


  INITIALISATIONS                                            
  ===============                                            

 REINITIALISATION DE LA PRESSION A L'ITERATION          1

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


  COMPUTATION OF PHYSICAL QUANTITIES                         
  ==================================                         

  Values of f3 min and max:    1.00000000000000        1.00000000000000     
  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                         = 141100
 Nb points avec PDF rectangle sans Dirac                =      0
 - - - - - - - - - -- - - -  et Dirac en FMINI          =      0
 - - - - - - - - - -- - - - - - - - - -  FMAXI          =      0
 - - - - - - - - - - - - - - - Diracs en FMINI et FMAXI =      0


MODELISATION DE LA COMBUSTION AVEC LE MODELE DE DIFFUSION TURBULENTE (CPCYM2)
CHIMIE RAPIDE A 3 CORPS - EXTENSION A 3 COMBUSTIBLES (Application au FUEL)
============================================================================
 Nb de points de calculs                                     =    141100
 Nb de points turbulents (passage par les PDF)               =         0

CONTROLE DES VALEURS DES FRACTIONS MASSIQUES
 Nb de points de calculs qui respectent somme des Yi = 1    =    141100
 Nb de points YCHX1 , YCHX2   < 0 ou > 1 =         0        0
 Nb de points YC0   , YH2S    < 0 ou > 1 =         0        0
 Nb de points YH2   , YHCN    < 0 ou > 1 =         0        0
 Nb de points YNH3            < 0 ou > 1 =         0
 Nb de points YO2   , YN2     < 0 ou > 1 =         0        0
 Nb de points YCO2  , YH2O    < 0 ou > 1 =         0        0
 Nb de points YSO2            < 0 ou > 1 =         0
  Somme Min MAX =    1.00000000000000        1.00000000000000     
 -----------------------------------------
 Property           Min. value  Max. value
 -----------------------------------------
  density           0.4394E+00  0.4394E+00
  molecular_viscos  0.3700E-04  0.3700E-04
  turbulent_viscos  0.3955E-11  0.3955E-11
 -----------------------------------------

 --- Diffusivity:
 -----------------------------------------------
 Scalar           Number  Min. value  Max. value
 -----------------------------------------------
 Enthalpy              1  0.4250E-04  0.4250E-04
 NG_FOL01              2  0.3700E-04  0.3700E-04
 NG_FOL02              3  0.3700E-04  0.3700E-04
 NG_FOL03              4  0.3700E-04  0.3700E-04
 YFOL_FOL01            5  0.3700E-04  0.3700E-04
 YFOL_FOL02            6  0.3700E-04  0.3700E-04
 YFOL_FOL03            7  0.3700E-04  0.3700E-04
 HLF_FOL01             8  0.3700E-04  0.3700E-04
 HLF_FOL02             9  0.3700E-04  0.3700E-04
 HLF_FOL03            10  0.3700E-04  0.3700E-04
 Fr_VAP               11  0.3700E-04  0.3700E-04
 Fr_HET               12  0.3700E-04  0.3700E-04
 Var_CB               13  0.3700E-04  0.3700E-04
 FR_HCN               14  0.3700E-04  0.3700E-04
 FR_NO                15  0.3700E-04  0.3700E-04
 Enth_Ox              16  0.3700E-04  0.3700E-04
 -----------------------------------------------


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


  COMPUTATION OF CFL, FOURIER AND VARIABLE DT                
  ===========================================                


 COURANT MAX=  0.0000E+00 IN -0.6552E-02 -0.5295E-03  0.4701E-01
 COURANT MIN=  0.0000E+00 IN -0.6552E-02 -0.5295E-03  0.4701E-01

 FOURIER MAX=  0.1549E+00 IN -0.9630E-04  0.1199E-03  0.4509E-03
 FOURIER MIN=  0.1791E-01 IN -0.6505E-02 -0.7412E-02  0.9618E-01

 COU/FOU MAX=  0.1549E+00 IN -0.9630E-04  0.1199E-03  0.4509E-03
 COU/FOU MIN=  0.1791E-01 IN -0.6505E-02 -0.7412E-02  0.9618E-01

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


  SETTING UP THE BOUNDARY CONDITIONS                         
  ==================================                         

