command: 
 ./cs_solver --mpi

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

                                  (R)
                      Code_Saturne

                      Version 3.2.2


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

  build Wed 16 Apr 2014 16:41:41 BST
  MPI version 2.1 (Open MPI 1.7.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.

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


Local case configuration:

  Date:                Wed 23 Apr 2014 14:11:50 BST
  System:              Linux 3.11.6-4-desktop
  Machine:             mbgnhma2lx
  Processor:           model name	: Intel(R) Core(TM) i7-3820 CPU @ 3.60GHz
  Memory:              32187 MB
  User:                mbgnhma2 (mbgnhma2)
  Directory:           /home/mbgnhma2/Code_saturne/Mixer/TurboModule/2D-Mixer-Keps/RESU/20140423-1411
  MPI ranks:           4 (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

 Reading metadata from file: "mesh_input/mesh_01"
 Reading metadata from file: "mesh_input/mesh_02"
 Reading metadata from file: "mesh_input/mesh_03"
 Reading metadata from file: "mesh_input/mesh_04"
 Reading metadata from file: "mesh_input/mesh_05"

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

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


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




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


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

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

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


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

 --- Physics
       NVAR   =              6 (Nb variables                )
       NSCAL  =              0 (Nb scalars                  )
       NSCAUS =              0 (Nb user scalars             )
       NSCAPP =              0 (Nb specific physics scalars )
       NPROCE =              6 (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         )

       OMEGAX =    0.00000E+00 (Rotation vector x component )
       OMEGAY =    0.00000E+00 (Rotation vector y component )
       OMEGAZ =    0.00000E+00 (Rotation vector z component )
       ICORIO =          0 (Coriolis source terms       )

  -- Continuous phase:

       RO0    =    0.99820E+03 (Reference density           )
       VISCL0 =    0.10030E-02 (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 =              0 (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

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


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

   - Commons
       ITHERM =              0 (0:  no thermal model        )
                               (10: temperature C           )
                               (11: temperature K           )
                               (20: enthalpy                )
                               (30: total energy            )
       ISCALT =             -1 (Thermal scalar number       )


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

 --- Continuous phase:

   - Commons
       ITURB  =             21 (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)     )

   - Linear production k-epsilon (ITURB = 21)
       ALMAX  =   -0.99900E+03 (Characteristic length       )
       UREF   =    0.17233E+00 (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            )


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


 ** 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-03 (Minimum time step           )
       DTMAX  =    0.10000E+01 (Maximum time step           )
       DTREF  =    0.10000E-02 (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
------------------------------------
 Pressure              0  0.1000E+01
 VelocityX             1  0.1000E+01
 VelocityY             1  0.1000E+01
 VelocityZ             1  0.1000E+01
 Turb Kinetic Ene      1  0.1000E+01
 Turb Dissipation      1  0.1000E+01
----------------------------

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

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

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


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

---------------------------------------------------------------------
 Variable          ICONV  IDIFF IDIFFT ISCHCV ISSTPC   BLENCV  THETAV
---------------------------------------------------------------------
 Pressure              0      1      1      1      0 0.00E+00 0.10E+01
 VelocityX             1      1      1      1      0 0.10E+01 0.10E+01
 VelocityY             1      1      1      1      0 0.10E+01 0.10E+01
 VelocityZ             1      1      1      1      0 0.10E+01 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
-------------------------------------------------------------

       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
    ------
       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
-------------------------------------------------------------------
 Pressure            100      2  0.1000E-04  0.1000E-06  0.0000E+00
 VelocityX           100      1  0.1000E-04  0.1000E-06  0.0000E+00
 VelocityY           100      1  0.1000E-04  0.1000E-06  0.0000E+00
 VelocityZ           100      1  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
-----------------------------------------------------------

-------------------------------------------
 Variable         IRCFLU IMLIGR      CLIMGR
-------------------------------------------
 Pressure              1     -1  0.1500E+01
 VelocityX             1     -1  0.1500E+01
 VelocityY             1     -1  0.1500E+01
 VelocityZ             1     -1  0.1500E+01
 Turb Kinetic Ene      1     -1  0.1500E+01
 Turb Dissipation      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                )


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


 ** CALCULATION OF TEMPORAL MEANS (MOMENTS)
    ---------------------------------------

       NBMOMT =              0 (Number of moments           )

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


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

--------------------------------------------------
 Variable         IRESOL NITMAX      EPSILO IDIRCL
--------------------------------------------------
 Pressure             -1  10000  0.1000E-07      1
 VelocityX            -1  10000  0.1000E-07      1
 VelocityY            -1  10000  0.1000E-07      1
 VelocityZ            -1  10000  0.1000E-07      1
 Turb Kinetic Ene     -1  10000  0.1000E-07      1
 Turb Dissipation     -1  10000  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    )


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


 ** MULTIGRID
    ---------

       NCEGRM =             30 (Max nb cells coarsest grid  )
       NGRMAX =             25 (Max number of levels        )
--------------------------------------
 Variable           IMGR NCYMAX NITMGF
--------------------------------------                       
 Pressure              0    100     10
 VelocityX             0    100     10
 VelocityY             0    100     10
 VelocityZ             0    100     10
 Turb Kinetic Ene      0    100     10
 Turb Dissipation      0    100     10
------------------------------

       IMGR   =  0 ou 1        (1: multigrid activated      )
       NCYMAX =                (Max number  of cycles       )
       NITMGF =                (Max nb iter on coarsest grid)


  Multigrid rank merge parameters:
    mean  coarse cells merge threshold: 300
    total coarse cells merge threshold: 500
    minimum ranks merge threshold:      1
    merge stride:                       1

  Multigrid coarsening type: algebraic, face traveral by criteria

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


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

 --- Restarted calculation
       ISUITE =              0 (1: restarted calculation    )
       ILEAUX =              1 (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 =             50 (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

       Number Name
         2 Pressure        
         3 VelocityX       
         4 VelocityY       
         5 VelocityZ       
         6 Turb Kinetic Ene
         7 Turb Dissipation
         8 Density         
        10 Turb Viscosity  
         --           --

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

       Number Name                  Nb. probes (-1: all)
         2 Pressure                      -1
         3 VelocityX                     -1
         4 VelocityY                     -1
         5 VelocityZ                     -1
         6 Turb Kinetic Ene              -1
         7 Turb Dissipation              -1
        10 Turb Viscosity                -1
         --           --                --

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

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

         2 Pressure                       2
         3 VelocityX                      2
         4 VelocityY                      2
         5 VelocityZ                      2
         6 Turb Kinetic Ene               1
         7 Turb Dissipation               1
        10 Turb Viscosity              -999
        11 CFL                         -999
        12 Fourier Number              -999
         --           --                --

 --- Additional post-processing variables (ipstdv)
       ipstfo =              0 (Force exerted  at boundary)
       ipstyp =              0 (Yplus          at boundary)
       ipstft =              0 (Thermal flow   at boundary)
       ipsttb =              0 (Temperature    at boundary)
       ipstnu =              0 (Nusselt        at boundary)


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


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

       IALE   =              0 (1: activated                )
       NALINF =              0 (Fluid initialization
                                                  iterations)


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


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

  -1: name: results
      directory: postprocessing
      format: EnSight Gold
      options: 
      time dependency: transient connectivity
      output: every 25 time steps and at calculation end


 Reading file:        mesh_input/mesh_01
 Finished reading:    mesh_input/mesh_01

 Reading file:        mesh_input/mesh_02
 Finished reading:    mesh_input/mesh_02

 Reading file:        mesh_input/mesh_03
 Finished reading:    mesh_input/mesh_03

 Reading file:        mesh_input/mesh_04
 Finished reading:    mesh_input/mesh_04

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

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

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

 Writing file:        partition_output/domain_number_4
 Finished writing:    partition_output/domain_number_4

 Partitioning finished (1.3 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 =              44711
    maximum value =              44712

      1 : [      44711 ;      44712 ] =          4

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

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

    minimum value =              45075
    maximum value =              45141

      1 : [      45075 ;      45088 [ =          1
      2 : [      45088 ;      45101 [ =          1
      3 : [      45101 ;      45114 [ =          0
      4 : [      45114 ;      45127 [ =          1
      5 : [      45127 ;      45141 ] =          1

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

 Histogram of the number of ghost cells per rank:

    minimum value =                364
    maximum value =                430

      1 : [        364 ;        377 [ =          1
      2 : [        377 ;        390 [ =          1
      3 : [        390 ;        403 [ =          0
      4 : [        403 ;        416 [ =          1
      5 : [        416 ;        430 ] =          1

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

 Histogram of the number of interior faces per rank:

    minimum value =              87138
    maximum value =              87177

      1 : [      87138 ;      87145 [ =          2
      2 : [      87145 ;      87153 [ =          0
      3 : [      87153 ;      87161 [ =          0
      4 : [      87161 ;      87169 [ =          0
      5 : [      87169 ;      87177 ] =          2

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

 Histogram of the number of boundary faces per rank:

    minimum value =              93366
    maximum value =              93451

      1 : [      93366 ;      93383 [ =          1
      2 : [      93383 ;      93400 [ =          0
      3 : [      93400 ;      93417 [ =          1
      4 : [      93417 ;      93434 [ =          1
      5 : [      93434 ;      93451 ] =          1

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

 Histogram of the number of interior faces per cell:

    minimum value =                  2
    maximum value =                  4

      1 : [          2 ;          3 [ =       2016
      2 : [          3 ;          4 ] =     176829

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

 Histogram of the number of neighboring domains per rank:

    minimum value =                  2
    maximum value =                  3

      1 : [          2 ;          3 ] =          4

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

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

 Writing file:        mesh_output
 Finished writing:    mesh_output

 Renumbering mesh:

Checking mesh renumbering for threads:
-------------------------------------


 Histogram of volume mesh matrix bandwidth per rank:

    minimum value =              44821
    maximum value =              44915

      1 : [      44821 ;      44839 [ =          1
      2 : [      44839 ;      44858 [ =          0
      3 : [      44858 ;      44877 [ =          0
      4 : [      44877 ;      44896 [ =          0
      5 : [      44896 ;      44915 ] =          3

 Histogram of volume mesh matrix profile/lines per rank:

    minimum value =               1338
    maximum value =               1916

      1 : [       1338 ;       1453 [ =          1
      2 : [       1453 ;       1569 [ =          1
      3 : [       1569 ;       1684 [ =          0
      4 : [       1684 ;       1800 [ =          0
      5 : [       1800 ;       1916 ] =          2

 Mesh coordinates:               minimum    and maximum
                       X : -4.9222948e-02  4.9222948e-02
                       Y : -4.9215902e-02  4.9217488e-02
                       Z :  0.0000000e+00  1.0000000e-02
 Mesh
     Number of cells:          178845
     Number of interior faces: 347724
     Number of boundary faces: 373638
     Number of vertices:       369362

 Groups:
    "fluid_gap1"
       cells:                 14400
    "fluid_gap2"
       cells:                 31824
    "fluid_holes"
       cells:                  6432
    "fluid_inlet"
       cells:                  3800
    "fluid_outside1"
       cells:                 14174
    "fluid_outside2"
       cells:                 26498
    "fluid_rotor1"
       cells:                 34489
    "fluid_rotor2"
       cells:                 47228
    "inlet"
       boundary faces:          380
    "interface1"
       boundary faces:          760
    "interface2"
       boundary faces:         2400
    "interface3"
       boundary faces:         1552
    "interface4"
       boundary faces:         4032
    "outlet"
       boundary faces:           11
    "outlet_wall"
       boundary faces:          145
    "rotor1_wall"
       boundary faces:         1212
    "rotor2_wall"
       boundary faces:          968
    "stator_wall"
       interior faces:          800
       boundary faces:         4488
    "symmetry1"
       boundary faces:       178845
    "symmetry2"
       boundary faces:       178845

 --- Information on the volumes
       Minimum control volume      =  1.1740955e-11
       Maximum control volume      =  3.1779964e-08
       Total volume for the domain =  6.5027603e-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: 31 -->   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 %

 Warning:
 --------
    Mesh quality issue has been detected

    The mesh should be re-considered using the listed criteria.

    The calculation will run but the solution quality may be degraded...

 Computing geometric quantities (0.0869 s)
                                                             
       ALMAX  =    0.40342E-01 (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
 -----------------------------------------                   
  Pressure          0.0000E+00  0.0000E+00
  VelocityX         0.0000E+00  0.0000E+00
  VelocityY         0.0000E+00  0.0000E+00
  VelocityZ         0.0000E+00  0.0000E+00
  Turb Kinetic Ene  0.1782E-04  0.1782E-04
  Turb Dissipation  0.1678E-06  0.1678E-06
 ---------------------------------

  Local Time Step   0.1000E-02  0.1000E-02
 ---------------------------------



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

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



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



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


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




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


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



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


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

 REINITIALISATION DE LA PRESSION A L'ITERATION          1

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


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

 -----------------------------------------
 Property           Min. value  Max. value
 -----------------------------------------
  density           0.9982E+03  0.9982E+03
  Laminar Viscosit  0.1003E-02  0.1003E-02
  Turb Viscosity    0.1700E+00  0.1700E+00
 -----------------------------------------


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


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


 COURANT MAX=  0.0000E+00 IN -0.3205E-01 -0.3474E-03  0.5000E-02
 COURANT MIN=  0.0000E+00 IN -0.3205E-01 -0.3474E-03  0.5000E-02

 FOURIER MAX=  0.1626E+04 IN  0.2246E-01 -0.2252E-01  0.5000E-02
 FOURIER MIN=  0.1767E+00 IN  0.4840E-01  0.2515E-02  0.5000E-02

 COU/FOU MAX=  0.1626E+04 IN  0.2246E-01 -0.2252E-01  0.5000E-02
 COU/FOU MIN=  0.1767E+00 IN  0.4840E-01  0.2515E-02  0.5000E-02

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


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



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

-------------------------------------------------------------------------
Boundary type          Code    Nb faces
-------------------------------------------------------------------------
Inlet                     2         380
Smooth wall               5        6813
Rough wall                6           0
Symmetry                  4      357690
Free outlet               3          11
Free inlet               13           0
Undefined                 1           0
-------------------------------------------------------------------------


Boundary faces with free inlet/outlet detected
Update of reference point for total pressure
 XYZP0 =   -0.18302E-02   0.49183E-01   0.50000E-02



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

---------------------------------------------------------------
Boundary type          Code    Nb faces           Mass flow
---------------------------------------------------------------
Inlet                     2         380         0.000000000E+00
Smooth wall               5        6813         0.000000000E+00
Rough wall                6           0         0.000000000E+00
Symmetry                  4      357690         0.000000000E+00
Free outlet               3          11         0.000000000E+00
Free inlet               13           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.87245E+00
   Friction velocity        uet   :  0.00000E+00 0.13875E+00
   Friction velocity        uk    :  0.14256E-03 0.23498E+00
   Dimensionless distance   yplus :  0.14150E-02 0.11250E+02
   ------------------------------------------------------
   Nb of reversal of the velocity at the wall   :          0
   Nb of faces within the viscous sub-layer     :       5373
   Total number of wall faces                   :       6813
------------------------------------------------------------


@                                                            
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@                                                            
@ @@ 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.        
@                                                            
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@                                                            

Incoming flow detained for          0 outlet faces on         11

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


  SOLVING NAVIER-STOKES EQUATIONS                            
  ===============================                            


   ** SOLVING VELOCITY
      ----------------

 velocity : CONVECTION IN   CENTERED  BLENDING WITH   0. % OF UPWIND
 _iterative_vector_gradient: isweep = 6, residue norm: 8.853339e-07, norm: 4.772334e+05, var: Var.  2
 velocity :    1038496 FACES WITH UPWIND ON     352096 INTERIOR FACES 
Jacobi [velocity]:
  n_iter :  4822, res_abs :  2.3331e-09, res_nor :  9.9939e-09
 velocity : CONVECTION IN   CENTERED  BLENDING WITH   0. % OF UPWIND
 _iterative_vector_gradient: isweep = 6, residue norm: 1.022081e-06, norm: 3.571213e+05, var: Var.  2
 velocity :      82200 FACES WITH UPWIND ON     352096 INTERIOR FACES 
@
@ @@ WARNING: velocity CONVECTION-DIFFUSION-SOURCE TERMS
@    ========
@  Maximum number of iterations          1 reached
@

 Maximum velocity after prediction   0.2411E+01
 --- Information on the volumes
       Minimum control volume      = -2.1882409e-09
       Maximum control volume      =  3.1779964e-08
       Total volume for the domain =  6.5028220e-05

Abort due to the detection of a negative control volume.

   ** SOLVING CONTINUITY PRESSURE
      ---------------------------

 pressure         : NORMED RESIDUALS =   0.593417E-01
 _iterative_vector_gradient: isweep = 100, residu norm: 2.185700e+84, norm: 3.494363e+09, var: = Var.  2
@ @@ warning: non convergence of grdvec
 pressure         : SWEEP =     1 RIGHT HAND SIDE NORM =   0.801285E+85, RELAXP =   0.100000E+01
Conjugate gradient [pressure]:


Conjugate gradient [pressure]: divergence after 2033 iterations:
  initial residual:  8.0128e+84; current residual:  1.2539e+89



cs_sles.c:2993: Fatal error.

Conjugate gradient: error (divergence) solving for pressure


Call stack:
   1: 0x7f516c9d1766 <reslin_+0x5e6>                  (libsaturne.so.0)
   2: 0x7f516c9eadf6 <invers_+0x1b6>                  (libsaturne.so.0)
   3: 0x7f516c8d0415 <resopv_+0x6225>                 (libsaturne.so.0)
   4: 0x7f516c8b4549 <navstv_+0x4b69>                 (libsaturne.so.0)
   5: 0x7f516c8df851 <tridim_+0x4801>                 (libsaturne.so.0)
   6: 0x7f516c7c7c8c <caltri_+0x3c5c>                 (libsaturne.so.0)
   7: 0x7f516c7b00f5 <cs_run+0x405>                   (libsaturne.so.0)
   8: 0x7f516c7afbda <main+0x14a>                     (libsaturne.so.0)
   9: 0x7f516ac47be5 <__libc_start_main+0xf5>         (libc.so.6)
  10: 0x403a19     <>                               (cs_solver)
End of stack

