command: 
 ./cs_solver --param Calc --mpi

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

                                  (R)
                      Code_Saturne

                      Version 3.2.1


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

  build Mi 15 Jan 2014 03:54:49 CET
  MPI version 2.2 (MPICH2 1.4.1)


  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:                Mi 15 Jan 2014 04:39:18 CET
  System:              Linux 3.11.0-15-generic
  Machine:             sebastian-GA-78LMT-USB3
  Processor:           model name	: AMD FX(tm)-4300 Quad-Core Processor
  Memory:              3431 MB
  User:                sebastian (sebastian)
  Directory:           /home/sebastian/Code_Saturne/3.2.1/arch/Linux_x86_64/bin/RPM2500/Case1/RESU/20140115-0439
  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

  External libraries for partitioning:
    ParMETIS 4.0.3
    SCOTCH 5.1.10

 Reading metadata from file: "mesh_input/mesh_01"
 Reading metadata from file: "mesh_input/mesh_02"

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

                   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.10350E+04 (Reference density           )
       VISCL0 =    0.35000E-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  =             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.99900E+03 (Characteristic length       )
       UREF   =    0.10000E+01 (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-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
 TurbEner              1  0.1000E+01
 Dissip                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
 TurbEner              1      1      1      1      0 0.00E+00 0.10E+01
 Dissip                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
 TurbEner            100      1  0.1000E-04  0.1000E-06  0.0000E+00
 Dissip              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
 TurbEner              1     -1  0.1500E+01
 Dissip                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              0  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
 TurbEner             -1  10000  0.1000E-07      1
 Dissip               -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              1    100     10
 VelocityX             0    100     10
 VelocityY             0    100     10
 VelocityZ             0    100     10
 TurbEner              0    100     10
 Dissip                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 =             10 (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 TurbEner        
         7 Dissip          
        10 TurbVisc        
        11 CourantNb       
        12 FourierNb       
        13 total_pressure  
         --           --

 --- 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 Pressure                      -1
         3 VelocityX                     -1
         4 VelocityY                     -1
         5 VelocityZ                     -1
         6 TurbEner                      -1
         7 Dissip                        -1
         8 Density                       -1
         9 LamVisc                       -1
        10 TurbVisc                      -1
        11 CourantNb                     -1
        12 FourierNb                     -1
        13 total_pressure                -1
         --           --                --

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

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

         2 Pressure                       0
         3 VelocityX                      0
         4 VelocityY                      0
         5 VelocityZ                      0
         6 TurbEner                       0
         7 Dissip                         0
        10 TurbVisc                    -999
        11 CourantNb                   -999
        12 FourierNb                   -999
        13 total_pressure              -999
         --           --                --

 --- Additional post-processing variables (ipstdv)
       ipstfo =              1 (Force exerted  at boundary)
       ipstyp =              1 (Yplus          at boundary)
       ipstft =              1 (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: binary
      time dependency: fixed mesh
      output: 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
 No "partition_input/domain_number_4" file available;

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

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

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

 Partitioning finished (3.35 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 =             118176
    maximum value =             118177

      1 : [     118176 ;     118177 ] =          4

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

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

    minimum value =             121237
    maximum value =             123327

      1 : [     121237 ;     121655 [ =          2
      2 : [     121655 ;     122073 [ =          0
      3 : [     122073 ;     122491 [ =          1
      4 : [     122491 ;     122909 [ =          0
      5 : [     122909 ;     123327 ] =          1

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

 Histogram of the number of ghost cells per rank:

    minimum value =               3060
    maximum value =               5150

      1 : [       3060 ;       3478 [ =          2
      2 : [       3478 ;       3896 [ =          0
      3 : [       3896 ;       4314 [ =          1
      4 : [       4314 ;       4732 [ =          0
      5 : [       4732 ;       5150 ] =          1

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

 Histogram of the number of interior faces per rank:

    minimum value =             223171
    maximum value =             230224

      1 : [     223171 ;     224581 [ =          1
      2 : [     224581 ;     225992 [ =          1
      3 : [     225992 ;     227402 [ =          1
      4 : [     227402 ;     228813 [ =          0
      5 : [     228813 ;     230224 ] =          1

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

 Histogram of the number of boundary faces per rank:

    minimum value =              18046
    maximum value =              30419

      1 : [      18046 ;      20520 [ =          1
      2 : [      20520 ;      22995 [ =          1
      3 : [      22995 ;      25469 [ =          0
      4 : [      25469 ;      27944 [ =          0
      5 : [      27944 ;      30419 ] =          2

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

 Histogram of the number of interior faces per cell:

    minimum value =                  2
    maximum value =                  4

      1 : [          2 ;          3 [ =        322
      2 : [          3 ;          4 ] =     472385

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

 Histogram of the number of neighboring domains per rank:

    minimum value =                  3
    maximum value =                  3

      1 : [          3 ;          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 =             120729
    maximum value =             123027

      1 : [     120729 ;     121188 [ =          2
      2 : [     121188 ;     121648 [ =          1
      3 : [     121648 ;     122107 [ =          0
      4 : [     122107 ;     122567 [ =          0
      5 : [     122567 ;     123027 ] =          1

 Histogram of volume mesh matrix profile/lines per rank:

    minimum value =               4607
    maximum value =              35263

      1 : [       4607 ;      10738 [ =          1
      2 : [      10738 ;      16869 [ =          1
      3 : [      16869 ;      23000 [ =          1
      4 : [      23000 ;      29131 [ =          0
      5 : [      29131 ;      35263 ] =          1

 Mesh coordinates:               minimum    and maximum
                       X : -2.9999674e+01  6.1540000e+01
                       Y : -3.3796631e+01  2.9998701e+01
                       Z :  0.0000000e+00  2.9000000e+01
 Mesh
     Number of cells:          472707
     Number of interior faces: 895724
     Number of boundary faces: 99380
     Number of vertices:       105536

 Groups:
    "Inlet"
       boundary faces:           27
    "Outlet"
       boundary faces:           27
    "boundary"
       boundary faces:        12872
    "cells"
       cells:                255133
    "wall"
       boundary faces:        86454

 --- Information on the volumes
       Minimum control volume      =  5.8030726e-08
       Maximum control volume      =  2.2106987e+01
       Total volume for the domain =  2.1645694e+04

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

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

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

  Criterion 4: Cells Volume Ratio:
    Number of bad cells detected: 10 -->   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.782 s)
  Global number of boundary faces selected for joining:      12872

  Global number of single vertices found:     21
  Global number of single edges found:         8
  Global min/max. tolerance:

 Glob. Num. |  Tolerance  |              Coordinates

       6980 |    0.060034 |  2.5340295569e+01   1.0615101752e-01   8.5000000000e+00 | ORI
       6807 |    0.596570 | -1.0292240615e+01   1.3522274114e+01   8.5000000000e+00 | ORI
  Determination of possible face intersections:

    bounding-box tree layout: 3D

  Global number of intersections detected:        56171
    Vertex-Vertex intersections:            5986
    Other intersections:                   50185

  Global number of new vertices to create:      38534

  Equivalences broken for 29 edges.

  Edge intersections and vertex creation done.

  Pre-merge for 51990 global element couples.

  Max size of a merge set of vertices: 9

  Excessive transitivity for 0 set(s) of vertices.

  Synchronization of the "single" elements after the merge step.

  Edge removed for 273 faces (global).
  Join mesh cleaning done.

  Edge removed for 143 faces (global).
  Join mesh cleaning done.

  Merge vertices and mesh update done.

  *** WARNING ***
  Globally, 2 problem(s) found during the face splitting
                1  open cycles,
                1  edges traversed twice,
                0  faces split into more than max_subfaces (= 200)

    => Eventually modify joining parameters


  Mesh cleaning done for degenerate faces.
    Global number of cleaned interior faces:        0
    Global number of cleaned border faces:          0
                                                             
       ALMAX  =    0.27775E+02 (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
  TurbEner          0.6000E-03  0.6000E-03
  Dissip            0.4762E-07  0.4762E-07
 ---------------------------------

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



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

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



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



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


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




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


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


 -----------------------------------------
 Property           Min. value  Max. value
 -----------------------------------------
  density           0.1035E+04  0.1035E+04
  LamVisc           0.3500E-02  0.3500E-02
  TurbVisc          0.7041E+03  0.7041E+03
 -----------------------------------------



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

-------------------------------------------------------------------------
Boundary type          Code    Nb faces
-------------------------------------------------------------------------
Inlet                     2          27
Smooth wall               5       86454
Rough wall                6           0
Symmetry                  4          32
Free outlet               3          27
Free inlet               13           0
Undefined                 1           0
-------------------------------------------------------------------------


Boundary faces with free inlet/outlet detected
Update of reference point for total pressure
 XYZP0 =    0.61540E+02  -0.23178E+02   0.45732E+01



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

---------------------------------------------------------------
Boundary type          Code    Nb faces           Mass flow
---------------------------------------------------------------
Inlet                     2          27         0.000000000E+00
Smooth wall               5       86454         0.000000000E+00
Rough wall                6           0         0.000000000E+00
Symmetry                  4          32         0.000000000E+00
Free outlet               3          27         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.56978E+04
   Friction velocity        uet   :  0.00000E+00 0.44454E+03
   Friction velocity        uk    :  0.11879E-01 0.11935E+01
   Dimensionless distance   yplus :  0.81721E+01 0.78397E+04
   ------------------------------------------------------
   Nb of reversal of the velocity at the wall   :          0
   Nb of faces within the viscous sub-layer     :          0
   Total number of wall faces                   :      86454
------------------------------------------------------------


  Global number of boundary faces selected for joining:      12872

  Global number of single vertices found:     21
  Global number of single edges found:         8
  Global min/max. tolerance:

 Glob. Num. |  Tolerance  |              Coordinates

       6980 |    0.060034 |  2.4504322345e+01  -6.4560075898e+00   8.5000000000e+00 | ORI
       6807 |    0.596570 | -6.4417250450e+00   1.5725339187e+01   8.5000000000e+00 | ORI
  Determination of possible face intersections:

    bounding-box tree layout: 3D

  Global number of intersections detected:        58367
    Vertex-Vertex intersections:            6121
    Other intersections:                   52246

  Global number of new vertices to create:      38840

  Equivalences broken for 28 edges.

  Edge intersections and vertex creation done.

  Pre-merge for 54342 global element couples.

  Max size of a merge set of vertices: 10

  Excessive transitivity for 0 set(s) of vertices.

  Synchronization of the "single" elements after the merge step.

  Edge removed for 303 faces (global).
  Join mesh cleaning done.

  Edge removed for 149 faces (global).
  Join mesh cleaning done.

  Merge vertices and mesh update done.

  *** WARNING ***
  Globally, 1 problem(s) found during the face splitting
                0  open cycles,
                1  edges traversed twice,
                0  faces split into more than max_subfaces (= 200)

    => Eventually modify joining parameters


  Mesh cleaning done for degenerate faces.
    Global number of cleaned interior faces:        0
    Global number of cleaned border faces:          0

   ** INFORMATION ON UNSTEADY ROTOR/STATOR TREATMENT
      ----------------------------------------------
 Time dedicated to mesh update (s):      0.7931
 Global time                   (s):      0.8711

 Warning:
 --------
   _iterative_scalar_gradient; variable: Var.  0; sweeps: 100
                               normed residual:  7.3785e-01; norm:  5.1664e-12

   ** INFORMATION ON CONVERGENCE
      --------------------------
---------------------------------------------------------------
   Variable    Rhs norm      N_iter  Norm. residual      derive
---------------------------------------------------------------
c  Pressure     0.21465E+08    2845   0.47868E+01   0.10208E+01 
c  VelocityX    0.19063E+09     293   0.38787E-01   0.41896E+10 
c  VelocityY    0.19063E+09     293   0.38787E-01   0.79500E+09 
c  VelocityZ    0.19063E+09     293   0.38787E-01   0.30913E+09 
c  TurbEner     0.47517E+11     646   0.51996E-02   0.37231E+09 
c  Dissip       0.67272E+15     781   0.88294E-02   0.16078E+18 
---------------------------------------------------------------
  
  

  ** Clippings for computed fields
     -----------------------------

   field        initial min     initial max  clips to min  clips to max
-  --------  --------------  --------------  ------------  ------------
a  TurbEner
a  Dissip  


  ** Computed fields on cells
     ------------------------

   field                  minimum         maximum        set mean    spatial mean
-  --------------  --------------  --------------  --------------  --------------
v  Pressure            1.0728e+09      3.9163e+11      3.4908e+11      3.1596e+11
v  Velocity[x]             -71210      1.1622e+05          3214.4          723.71
v  Velocity[y]        -1.3891e+05           11379         -4529.2         -280.63
v  Velocity[z]             -55774           63766          570.19          34.898
v  TurbEner             0.0005985           20805          493.12          208.37
v  Dissip              4.1146e-08      3.1529e+09       1.747e+07      3.0106e+06
v  TurbVisc                704.15          704.15          704.15          704.15
v  CourantNb                    0               0               0               0
v  FourierNb            0.0010572          375.54          7.4543        0.026763
v  total_pressure      1.0729e+09      3.9163e+11      3.4908e+11      3.1596e+11


  ** Computed values on cells
     ------------------------
