command: 
 /opt/saturne-3.0/libexec/code_saturne/cs_solver --param LES.xml --mpi

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

                                  (R)
                      Code_Saturne

                      Version 3.0.0


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

  build Thu 02 May 2013 08:06:16 AM PKT
  MPI version 2.1 (Open MPI 1.4.3)


  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 07 Jun 2016 01:14:37 PM PKT
  System:              Linux 3.2.0-60-generic
  Machine:             fahad-virtual-machine
  Processor:           model name	: AMD Opteron(tm) Processor 6348
  Memory:              122453 MB
  User:                fahad (Fahad)
  Directory:           /home/fahad/Desktop/Code_Saturne_CircularCylinder/CIRCULARCYLINDER/LES/RESU/20160607-1311
  MPI ranks:           5 (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 3.1
    PT-SCOTCH 1.12.2147483631

 Reading metadata from file: "mesh_input"

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

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


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



@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@
@ @@   WARNING :      WHEN READING INPUT DATA
@    =========
@   CHOICE OF TIME-SCHEME
@
@     TIME-SCHEME FOR VELOCITY IS SECOND ORDER
@       (THETAV =   0.50E+00)
@     CERTAIN TERMES ARE HOWEVER FIRST ORDER IN TIME  WITH
@       THE FOLLOWING SETTINGS:
@
@ parameters       ISTMPF ISNO2T ISTO2T IROEXT IVIEXT ICPEXT
@ Values  entered      2      1      0      0      0      0
@
@  computation will go on.
@
@ Check the input data given through the User Interface
@   or in cs_user_parameters.f90.
@
@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@


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


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

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

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


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

 --- Geometry
       NPRFML =              1 (Nb max. family properties   )
       NFML   =              7 (Nb families                 )

 --- Physics
       NVAR   =              4 (Nb variables                )
       NSCAL  =              0 (Nb scalars                  )
       NSCAUS =              0 (Nb user scalars             )
       NSCAPP =              0 (Nb specific physics scalars )
       NPROCE =              9 (Nb cell properties          )
       NPROFA =              2 (Nb internal face properties )
       NPROFB =              3 (Nb boundary face 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.11786E+01 (Reference density           )
       VISCL0 =    0.18300E-04 (Ref. molecular dyn. visc.   )
       CP0    =    0.10172E+04 (Ref. specific heat          )
       ICP    =              0 (> 0: variable CP (usphyv)   )
       P0     =    0.10132E+06 (Ref. total pressure         )
       PRED0  =    0.00000E+00 (Ref. reduced pressure       )
       T0     =    0.29315E+03 (Ref. temperature            )

       IROVAR =              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

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


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

 --- Continuous phase:

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

   - LES                 (ITURB = 40, 41, 42)
                               (Sub-grid scale model        )
                               (40 Smagorinsky model        )
                               (41 Dynamic model            )
                               (42 WALE model               )
       CSMAGO =    0.65000E-01 (Smagorinsky constant        )
       CWALE  =    0.25000E+00 (WALE model constant         )
       XLESFL =    0.20000E+01 (Filter with in a cell is    )
       ALES   =    0.10000E+01 (written as                  )
       BLES   =    0.33333E+00 (XLESFL*(ALES*VOLUME)**(BLES))
       IDRIES =              1 (=1 Van Driest damping       )
       CDRIES =    0.26000E+02 (Van Driest constant         )
       XLESFD =    0.15000E+01 (Ratio between the explicit  )
                               (filter and LES filter       )
                               (recommended value: 1.5      )
       SMAGMX =    0.65000E+00 (Max Smagonsky in the        )
                               (dynamic model case          )
       IVRTEX =              0 (Use the vortex method       )
 --- 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))


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


 ** 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-04 (Minimum time step           )
       DTMAX  =    0.10000E+00 (Maximum time step           )
       DTREF  =    0.10000E-03 (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
----------------------------

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

 --- Order of base time stepping scheme
       ISCHTP =              2 (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      1 0.10E+01 0.50E+00
 VelocityY             1      1      1      1      1 0.10E+01 0.50E+00
 VelocityZ             1      1      1      1      1 0.10E+01 0.50E+00
-------------------------------------------------------------

       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
    ------
       IVELCO =              1 (0: segregated solving of the
                                   velocity components during
                                   the prediction step
                                1: coupled solving of the
                                   velocity components       )
       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 =              2 (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 =              2 (time scheme for flow
                               (0: explicit (THETFL = 0     )
                               (1: std scheme (Saturne 1.0  )
                               (2: 2nd-order (THETFL = 0.5  )
       THETFL =    0.50000E+00 (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.50000E+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      5  0.1000E-04  0.1000E-03  0.0000E+00
 VelocityX           100     10  0.1000E-04  0.1000E-03  0.0000E+00
 VelocityY           100     10  0.1000E-04  0.1000E-03  0.0000E+00
 VelocityZ           100     10  0.1000E-04  0.1000E-03  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
-----------------------------------

       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           )

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


 ** WALL DISTANCE COMPUTATION
    -------------------------

       ICDPAR =             -1 ( 1: std, reread if restart
                               (-1: std, recomputed if restrt
                               ( 2: old, reread if restart
                               (-2: old, recomputed if restrt


       NITMAY =          10000 (Nb iter for iter resolution )
       NSWRSY =              1 (Nb iter for rhs reconstr.   )
       NSWRGY =            100 (Nb iter for grad. reconstr. )
       IMLIGY =             -1 (Gradient limitation method  )
       IRCFLY =              1 (Conv. Diff. flow reconstr.  )
       ISCHCY =              1 (Convective scheme           )
       ISSTPY =              0 (Slope tet use               )
       IMGRPY =              0 (Multigrid algorithm         )
       IWARNY =              0 (Verbosity level             )
       NTCMXY =           1000 (Nb iter for steady convect. )

       BLENCY =    0.00000E+00 (2nd order conv. scheme prop.)
       EPSILY =    0.10000E-07 (Iterative solver precision  )
       EPSRSY =    0.10000E-04 (rhs reconstruction precision)
       EPSRGY =    0.10000E-04 (Gradient reconstr. precision)
       CLIMGY =    0.15000E+01 (Coeff. for grad. limitation )
       EXTRAY =    0.00000E+00 (Coeff. for grad. extrapolat.)
       COUMXY =    0.50000E+04 (Max CFL for convection      )
       EPSCVY =    0.10000E-07 (Precision for steady conv.  )
       YPLMXY =    0.20000E+03 (y+ max w. damping influence )


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


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

--------------------------------------------------
 Variable         IRESOL NITMAX      EPSILO IDIRCL
--------------------------------------------------
 Pressure              0  10000  0.1000E-04      1
 VelocityX            -1  10000  0.1000E-04      1
 VelocityY            -1  10000  0.1000E-04      1
 VelocityZ            -1  10000  0.1000E-04      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
------------------------------

       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 =              1 (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       
         9 CourantNb       
        10 FourierNb       
        11 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 Density                       -1
         7 LamVisc                       -1
         9 CourantNb                     -1
        10 FourierNb                     -1
        11 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
         9 CourantNb                   -999
        10 FourierNb                   -999
        11 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)


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


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

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

 Partitioning 11275730 cells to 5 domains on 5 ranks
   (ParMETIS_V3_PartKway).

  Total number of faces on parallel boundaries: 2296
  wall-clock time: 15.631742 s

  Number of cells per domain (histogramm):
    [    2232027 ;    2238608 [ =          3
    [    2238608 ;    2245189 [ =          0
    [    2245189 ;    2251771 [ =          0
    [    2251771 ;    2258352 [ =          0
    [    2258352 ;    2264934 [ =          0
    [    2264934 ;    2271515 [ =          1
    [    2271515 ;    2278096 [ =          0
    [    2278096 ;    2284678 [ =          0
    [    2284678 ;    2291259 [ =          0
    [    2291259 ;    2297841 ] =          1

 Writing file:        partition_output/domain_number_5
 Finished writing:    partition_output/domain_number_5

 Partitioning finished (39.9 s)


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

  -1: name: results
      directory: postprocessing
      format: EnSight Gold
      options: binary
      time dependency: fixed mesh
      output: every 1 time steps


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

 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 =            2232027
    maximum value =            2297841

      1 : [    2232027 ;    2245189 [ =          3
      2 : [    2245189 ;    2258352 [ =          0
      3 : [    2258352 ;    2271515 [ =          1
      4 : [    2271515 ;    2284678 [ =          0
      5 : [    2284678 ;    2297841 ] =          1

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

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

    minimum value =            2232962
    maximum value =            2298755

      1 : [    2232962 ;    2246120 [ =          3
      2 : [    2246120 ;    2259279 [ =          0
      3 : [    2259279 ;    2272437 [ =          1
      4 : [    2272437 ;    2285596 [ =          0
      5 : [    2285596 ;    2298755 ] =          1

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

 Histogram of the number of ghost cells per rank:

    minimum value =                737
    maximum value =               1018

      1 : [        737 ;        793 [ =          1
      2 : [        793 ;        849 [ =          0
      3 : [        849 ;        905 [ =          0
      4 : [        905 ;        961 [ =          2
      5 : [        961 ;       1018 ] =          2

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

 Histogram of the number of interior faces per rank:

    minimum value =            3338371
    maximum value =            3436044

      1 : [    3338371 ;    3357905 [ =          3
      2 : [    3357905 ;    3377440 [ =          0
      3 : [    3377440 ;    3396974 [ =          1
      4 : [    3396974 ;    3416509 [ =          0
      5 : [    3416509 ;    3436044 ] =          1

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

 Histogram of the number of boundary faces per rank:

    minimum value =            4484328
    maximum value =            4618031

      1 : [    4484328 ;    4511068 [ =          3
      2 : [    4511068 ;    4537809 [ =          0
      3 : [    4537809 ;    4564549 [ =          1
      4 : [    4564549 ;    4591290 [ =          0
      5 : [    4591290 ;    4618031 ] =          1

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

 Histogram of the number of interior faces per cell:

    minimum value =                  2
    maximum value =                  3

      1 : [          2 ;          3 ] =   11275730

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

 Histogram of the number of neighboring domains per rank:

    minimum value =                  2
    maximum value =                  3

      1 : [          2 ;          3 ] =          5

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

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

 Renumbering mesh:

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


 Histogram of volume mesh matrix bandwidth per rank:

    minimum value =            2232007
    maximum value =            2297772

      1 : [    2232007 ;    2245160 [ =          3
      2 : [    2245160 ;    2258313 [ =          0
      3 : [    2258313 ;    2271466 [ =          1
      4 : [    2271466 ;    2284619 [ =          0
      5 : [    2284619 ;    2297772 ] =          1

 Histogram of volume mesh matrix profile/lines per rank:

    minimum value =               1353
    maximum value =               1746

      1 : [       1353 ;       1431 [ =          1
      2 : [       1431 ;       1510 [ =          1
      3 : [       1510 ;       1588 [ =          0
      4 : [       1588 ;       1667 [ =          1
      5 : [       1667 ;       1746 ] =          2

 Mesh coordinates:               minimum    and maximum
                       X : -7.0000000e+00  1.5000000e+01
                       Y : -5.0000000e+00  5.0000000e+00
                       Z :  0.0000000e+00  1.0000000e+00
 Mesh
     Number of cells:          11275730
     Number of interior faces: 16863169
     Number of boundary faces: 22652312
     Number of vertices:       11376582

 Groups:
    "Back"
       boundary faces:     11275730
    "Bottom"
       boundary faces:          275
    "Front"
       boundary faces:     11275730
    "Left"
       boundary faces:          125
    "Top"
       boundary faces:          275
    "right"
       boundary faces:          125

 --- Information on the volumes
       Minimum control volume      =  5.5578274e-11
       Maximum control volume      =  5.2654104e-03
       Total volume for the domain =  2.1685841e+02

  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: 2113 -->   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 (8.94 s)
                                                             
       ALMAX  =    0.59972E+01 (Characteristic length       )
       ALMAX is the cubic root of the domain volume.


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


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

  Local Time Step   0.1000E-03  0.1000E-03
 ---------------------------------



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

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



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



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


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




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


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


 -----------------------------------------
 Property           Min. value  Max. value
 -----------------------------------------
  Density           0.1179E+01  0.1179E+01
  LamVisc           0.1830E-04  0.1830E-04
  Turb Viscosity    0.0000E+00  0.0000E+00
 -----------------------------------------



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

-------------------------------------------------------------------------
Boundary type          Code    Nb faces
-------------------------------------------------------------------------
Inlet                     2         125
Smooth wall               5           0
Rough wall                6      100052
Symmetry                  4    22552010
Free outlet               3         125
Undefined                 1           0
-------------------------------------------------------------------------


Boundary faces with free inlet/outlet detected
Update of reference point for total pressure
 XYZP0 =    0.15000E+02   0.61542E-11   0.50000E+00



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

---------------------------------------------------------------
Boundary type          Code    Nb faces           Mass flow
---------------------------------------------------------------
Inlet                     2         125         0.000000000E+00
Smooth wall               5           0         0.000000000E+00
Rough wall                6      100052         0.000000000E+00
Symmetry                  4    22552010         0.000000000E+00
Free outlet               3         125         0.000000000E+00
Undefined                 1           0         0.000000000E+00
---------------------------------------------------------------


   ** BOUNDARY CONDITIONS FOR ROUGH WALLS
   --------------------------------------
------------------------------------------------------------
                                         Minimum     Maximum
------------------------------------------------------------
   Rel velocity at the wall uiptn :  0.00000E+00 0.00000E+00
   Friction velocity        uet   :  0.14672E-11 0.93805E-10
   Friction velocity        uk    :  0.14672E-11 0.93805E-10
   Rough dimensionless dist yplus :  0.44874E-02 0.33143E+00
   ------------------------------------------------------   
   Nb of reversal of the velocity at the wall   :          0
------------------------------------------------------------


                                                             
 ** WALL DISTANCE                                            
    -------------                                            
                                                             
  Min distance =    0.44965E-06 Max distance =    0.15965E+02
                                                             
     (Distance calculation done in       7146 iterations)



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

---------------------------------------------------------------
Boundary type          Code    Nb faces           Mass flow
---------------------------------------------------------------
Inlet                     2         125        -0.785746667E+02
Smooth wall               5           0         0.000000000E+00
Rough wall                6      100052         0.000000000E+00
Symmetry                  4    22552010         0.000000000E+00
Free outlet               3         125         0.785408695E+02
Undefined                 1           0         0.000000000E+00
---------------------------------------------------------------


   ** BOUNDARY CONDITIONS FOR ROUGH WALLS
   --------------------------------------
------------------------------------------------------------
                                         Minimum     Maximum
------------------------------------------------------------
   Rel velocity at the wall uiptn :  0.00000E+00 0.00000E+00
   Friction velocity        uet   :  0.14672E-11 0.93805E-10
   Friction velocity        uk    :  0.14672E-11 0.93805E-10
   Rough dimensionless dist yplus :  0.44874E-02 0.33143E+00
   ------------------------------------------------------   
   Nb of reversal of the velocity at the wall   :          0
------------------------------------------------------------




 Write final restart files
   checkpoint at iteration          1, Physical time    0.10000E-03


   ** Writing the main restart file
      -----------------------------

 Start writing
   End writing the dimensions                           
   End writing the options                              
   End writing the main variables                       
 End writing

   ** Writing the auxilliary restart file
      -----------------------------------

 Start writing
   End writing the dimensions and options               
   End writing the physical properties                  
   End writing the time step                            
   End writing the mass fluxes                          
   End writing the boundary conditions                  
   End writing the source terms                         
   End writing the wall distance                        
 End writing


 Time for restart files:    0.15249E+02



   ** INFORMATION ON CONVERGENCE
      --------------------------
---------------------------------------------------------------
   Variable    Rhs norm      N_iter  Norm. residual      derive
---------------------------------------------------------------
c  Pressure     0.24368E+01     286   0.14808E-04   0.23125E+00 
c  VelocityX    0.75714E+05       2   0.45212E-04   0.10354E+07 
c  VelocityY    0.75714E+05       2   0.45212E-04   0.72642E+04 
c  VelocityZ    0.75714E+05       2   0.45212E-04   0.15014E-23 
---------------------------------------------------------------
  
  

   ** INFORMATION ON VARIABLES
      ------------------------
---------------------------------------------------------------
   Variable      Min. value    Max. value   Min clip   Max clip
---------------------------------------------------------------
v  Pressure     0.11442E+04   0.26694E+07         --         -- 
v  VelocityX   -0.10197E+02   0.25971E+02         --         -- 
v  VelocityY   -0.15013E+02   0.14818E+02         --         -- 
v  VelocityZ   -0.11562E-12   0.13425E-12         --         -- 
v  CourantNb    0.00000E+00   0.00000E+00          0          0 
v  FourierNb    0.17084E-05   0.17161E+03          0          0 
v  total_pressu 0.10247E+06   0.27707E+07          0          0 
---------------------------------------------------------------
  
  

   ** INFORMATION ON CLIPPINGS
      ------------------------
---------------------------------------------------------------
   Variable    Min wo clips  Max wo clips   Min clip   Max clip
---------------------------------------------------------------
---------------------------------------------------------------
  
  

   ** INFORMATION ON THE TIME STEP
      -----------------------------
---------------------------------------------------------------
Criterion    Value      at         xc           yc           zc
---------------------------------------------------------------
Courant min  0.00000E+00 -0.99751E+00 -0.70840E-01  0.50000E+00  
Courant max  0.00000E+00 -0.99751E+00 -0.70840E-01  0.50000E+00  
Fourier min  0.17084E-05 -0.48000E+01 -0.49794E+01  0.50000E+00  
Fourier max  0.17161E+03 -0.90219E+00 -0.43133E+00  0.50000E+00  
Cou_Fou min  0.17084E-05 -0.48000E+01 -0.49794E+01  0.50000E+00  
Cou_Fou max  0.17161E+03 -0.90219E+00 -0.43133E+00  0.50000E+00  
---------------------------------------------------------------
  
  

 TIME FOR THE TIME STEP                1:        0.44479E+04

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



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



                 FINAL STAGE OF THE CALCULATION              
                 ==============================              


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




   ** TIME FOR FINAL WRITING:    0.18985E+00
      -----------------------                                



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




                      END OF CALCULATION                     
                      ==================                     


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

 Destroying structures and ending computation


Calculation time summary:

  User CPU time:           3955.995 s
  System CPU time:          363.091 s
  Total CPU time:         22742.257 s

  Elapsed time:            4627.867 s
  CPU / elapsed time          0.933

Memory use summary:

  Total memory used:                              32.116 GiB
                             local minimum:        6.360 GiB  (rank 2)
                             local maximum:        6.542 GiB  (rank 3)
  Theoretical instrumented dynamic memory:        10.815 GiB
                             local minimum:        2.141 GiB  (rank 2)
                             local maximum:        2.204 GiB  (rank 3)
