command: 
 ./cs_solver --param restart_3.xml --mpi

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

                                  (R)
                      Code_Saturne

                      Version 4.0-beta


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

  revision 4.0-beta-r7278
  build Sun 15 Feb 2015 21:11:07 GMT
  MPI version 3.0 (Open MPI 1.8.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:                Mon 23 Mar 2015 16:44:15 GMT
  System:              Linux 3.0.101-0.29-default
  Machine:             node5
  Processor:           model name	: Intel(R) Xeon(R) CPU           X5650  @ 2.67GHz
  Memory:              24148 MB
  User:                ccytbf (Andrew Heffron)
  Directory:           /gpfs/home/external/ccytbf/run/cs-4.0-r7278/rae2822/CASE1/RESU/20150323-1644
  MPI ranks:           48 (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

 Reading metadata from file: "mesh_input"

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

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


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




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


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

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

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


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

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

       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.10133E+06 (Ref. total pressure         )
       PRED0  =    0.00000E+00 (Ref. reduced pressure       )
       T0     =    0.30000E+03 (Ref. temperature            )

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

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


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


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

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


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


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

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


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

 --- Continuous phase:

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

   - Spalart-Allmaras    (ITURB = 70)
       ALMAX  =   -0.99900E+03 (Characteristic length       )
       UREF   =    0.10000E+01 (Characteristic velocity     )
       RELAXV =    0.10000E+01 for nu (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))

   - Spalart-Allmaras    (ITURB = 70)
       CSAB1  =    0.13550E+00 (b1 constant                 )
       CSAB2  =    0.62200E+00 (b2 constant                 )
       CSASIG =    0.66667E+00 (sigma constant              )
       CSAV1  =    0.71000E+01 (v1 constant                 )
       CSAW1  =    0.32011E+01 (w1 constant                 )
       CSAW2  =    0.30000E+00 (w2 constant                 )
       CSAW3  =    0.20000E+01 (w3 constant                 )


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


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

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


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


 ** COMPRESSIBLE: additional data
    ------------

 --- Continuous phase :
       ICV    =              0 (0: Cv cst; 1: variable      )
       IVISCV =              0 (0: kappa cst; 1: variable
                                kappa: volume viscosity
                                in kg/(m.s)                 )
       VISCV0 =    0.00000E+00 (kappa value if constant     )
       ICFGRP =              1 (1: pressure BC with dominant
                                hydrostatic effect          )


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


 ** 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.50000E-07 (Minimum time step           )
       DTMAX  =    0.50000E-03 (Maximum time step           )
       DTREF  =    0.50000E-06 (Reference time step         )

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

 --- Frozen velocity field

       ICCVFG =              0 (1: frozen velocity field    )

 --- Per-variable properties

------------------------------------
 Variable          ISTAT      CDTVAR
------------------------------------
 Velocity              1  0.1000E+01
 Pressure              1  0.1000E+01
 nu_tilda              1  0.1000E+01
 TotEner               1  0.1000E+01
 temperature           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
---------------------------------------------------------------------
 Velocity              1      1      1      0      0 0.00E+00 0.10E+01
 Pressure              0      1      1      1      0 0.00E+00 0.10E+01
 nu_tilda              1      1      1      1      0 0.00E+00 0.10E+01
 TotEner               1      1      1      1      0 0.00E+00 0.10E+01
 temperature           1      1      1      1      0 0.00E+00 0.10E+01
-------------------------------------------------------------

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


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


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

  -- Continuous phase:

       IREVMC =              0 (Velocity reconstruction mode)

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

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


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


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

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

-------------------------------------------------------------------
 Variable         NSWRGR NSWRSM      EPSRGR      EPSRSM      EXTRAG
-------------------------------------------------------------------
 Velocity            100      1  0.1000E-04  0.1000E-07  0.0000E+00
 Pressure            100      2  0.1000E-04  0.1000E-06  0.0000E+00
 nu_tilda            100      1  0.1000E-04  0.1000E-06  0.0000E+00
 TotEner             100      1  0.1000E-04  0.1000E-07  0.0000E+00
 temperature         100      1  0.1000E-04  0.1000E-06  0.0000E+00
-----------------------------------------------------------

-------------------------------------------
 Variable         IRCFLU IMLIGR      CLIMGR
-------------------------------------------
 Velocity              1      1  0.1500E+01
 Pressure              1      1  0.1500E+01
 nu_tilda              1      1  0.1500E+01
 TotEner               1      1  0.1500E+01
 temperature           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                )


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


 ** 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               )
       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              EPSILO IDIRCL
------------------------------------
 Velocity          0.1000E-08      1
 Pressure          0.1000E-07      1
 nu_tilda          0.1000E-07      1
 TotEner           0.1000E-08      1
 temperature       0.1000E-07      1
------------------------------------

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


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


 ** SCALARS
    -------

       ITBRRB =              0 (T or H reconstruction at bdy)

-------------------------------------------------------------
 Variable         Number ISCACP ITURT      VISLS0      SIGMAS
-------------------------------------------------------------
 TotEner               1      0      0  0.2453E-04  0.1000E+01
 temperature           2      0      0  0.2495E-01  0.1000E+01
---------------------------------------------------------------------

------------------------------------
 Variable         Number      RVARFL
------------------------------------
 TotEner               1  0.8000E+00
 temperature           2  0.8000E+00
-------------------------------------------

-------------------------------------------------------
 Variable         Number ICLVFL      SCAMIN      SCAMAX
-------------------------------------------------------
 TotEner               1     -1 -0.1000E+13  0.1000E+13
 temperature           2     -1 -0.1000E+13  0.1000E+13
-------------------------------------------------------

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

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

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

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

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


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


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

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


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

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

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


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


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

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

 --- Post-processed variables


Velocity
Pressure
nu_tilda
TotEner
temperature
Density
Boundary Density
TurbVisc
CourantNb
FourierNb
         --           --

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

       Number Name                  Nb. probes (-1: all)
         2 Velocity[X]                   -1
         3 Velocity[Y]                   -1
         4 Velocity[Z]                   -1
         5 Pressure                      -1
         6 nu_tilda                      -1
         7 TotEner                       -1
         8 temperature                   -1
         9 Density                       -1
        10 LamVisc                       -1
        11 TurbVisc                      -1
        12 CourantNb                     -1
        13 FourierNb                     -1
         --           --                --

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

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

         2 Velocity                       0
         5 Pressure                       0
         6 nu_tilda                       0
         7 TotEner                        0
         8 temperature                    0
         9 Density                     -999
        11 TurbVisc                    -999
        12 CourantNb                   -999
        13 FourierNb                   -999
         --           --                --

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


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


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

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


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


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

  -1: name: results
      directory: postprocessing
      format: CGNS 3.2.1

      options: 
      time dependency: fixed mesh
      output: at calculation end


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

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

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

 Partitioning finished (10.2 s)


 -------------------------------------------------------
  Joining number 1:

  Periodicity type: translation
  Transformation matrix:             1            0            0            0
                                     0            1            0            0
                                     0            0            1        0.061

  Selection criteria: "SIDE_1 or SIDE_2"

  Parameters for the joining operation:
    Shortest incident edge fraction:           0.10000
    Maximum angle between joined face planes:  5.00000

  Advanced joining parameters:
    Verbosity level:                                 1
    Visualization level:                             1
    Deepest level reachable in tree building:       30
    Max boxes by leaf:                              25
    Max ratio of linked boxes / init. boxes:   5.00000
    Max ratio of boxes for distribution:       2.00000
    Merge step tolerance multiplier:           1.00000
    Pre-merge factor:                          0.05000
    Tolerance computation mode:                      1
    Intersection computation mode:                   1
    Max. number of equiv. breaks:                  500
    Max. number of subfaces by face:               200

  Before joining
     Number of cells:          7509150
     Number of interior faces: 22192659
     Number of boundary faces: 669582
     Number of vertices:       7845318

  Global number of boundary faces selected for joining:     600732

  Element selection successfully done.
  Global min/max. tolerance:

 Glob. Num. |  Tolerance  |              Coordinates

     278103 |    0.000016 | -1.6827923363e-04   2.6547694277e-05  -9.0897085046e-10 | ORI
     939953 |    0.000244 | -1.3863347637e+00  -1.4029420317e+00   0.0000000000e+00 | ORI
  Determination of possible face intersections:

    bounding-box tree layout: 3D

  Global number of intersections detected:     10153947
    Vertex-Vertex intersections:        10153947
    Other intersections:                       0

  Joining operation is conforming.

  Edge intersections and vertex creation done.

  Pre-merge for 656078 global element couples.

  No need to merge vertices.

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

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

  Merge vertices and mesh update done.

  Global number of faces after splitting:     901098

  Global configuration after the joining operation:
     Global number of border faces to add:          0

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

  After joining
     Number of cells:          7509150
     Number of interior faces: 22793391
     Number of boundary faces: 68850
     Number of vertices:       7845318


  Joining  1 completed (3.74 s).

  All joining operations successfully finished:

    Wall clock time:                   3.8


 ----------------------------------------------------------
 Composing periodicities

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

 Face interfaces creation
 Definition of periodic vertices
 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 =             156440
    maximum value =             156441

      1 : [     156440 ;     156441 ] =         48

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

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

    minimum value =             191132
    maximum value =             258329

      1 : [     191132 ;     204571 [ =         12
      2 : [     204571 ;     218010 [ =         15
      3 : [     218010 ;     231450 [ =         13
      4 : [     231450 ;     244889 [ =          4
      5 : [     244889 ;     258329 ] =          4

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

 Histogram of the number of ghost cells per rank:

    minimum value =              34692
    maximum value =             101889

      1 : [      34692 ;      48131 [ =         12
      2 : [      48131 ;      61570 [ =         15
      3 : [      61570 ;      75010 [ =         13
      4 : [      75010 ;      88449 [ =          4
      5 : [      88449 ;     101889 ] =          4

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

 Histogram of the number of interior faces per rank:

    minimum value =             485873
    maximum value =             521159

      1 : [     485873 ;     492930 [ =         12
      2 : [     492930 ;     499987 [ =         14
      3 : [     499987 ;     507044 [ =         13
      4 : [     507044 ;     514101 [ =          5
      5 : [     514101 ;     521159 ] =          4

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

 Histogram of the number of boundary faces per rank:

    minimum value =                  0
    maximum value =               4252

      1 : [          0 ;        850 [ =         12
      2 : [        850 ;       1700 [ =         20
      3 : [       1700 ;       2551 [ =          8
      4 : [       2551 ;       3401 [ =          7
      5 : [       3401 ;       4252 ] =          1

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

 Histogram of the number of interior faces per cell:

    minimum value =                  4
    maximum value =                  6

      1 : [          4 ;          5 [ =         50
      2 : [          5 ;          6 ] =    7509100

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

 Histogram of the number of neighboring domains per rank:

    minimum value =                 11
    maximum value =                 26

      1 : [         11 ;         14 [ =          4
      2 : [         14 ;         17 [ =         11
      3 : [         17 ;         20 [ =          6
      4 : [         20 ;         23 [ =         14
      5 : [         23 ;         26 ] =         13

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

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

 Writing file:        mesh_output
 Finished writing:    mesh_output

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

 Ignoring periodicity for graph-based partitioning.

 Partitioning 7509150 cells to 48 domains on 48 ranks
   (ParMETIS_V3_PartKway).

  Total number of faces on parallel boundaries: 209512
  wall-clock time: 0.327948 s

  Number of cells per domain (histogramm):
    [     106823 ;     115943 [ =          1
    [     115943 ;     125064 [ =          1
    [     125064 ;     134184 [ =          1
    [     134184 ;     143305 [ =          3
    [     143305 ;     152425 [ =         12
    [     152425 ;     161546 [ =         15
    [     161546 ;     170666 [ =          7
    [     170666 ;     179787 [ =          5
    [     179787 ;     188907 [ =          2
    [     188907 ;     198028 ] =          1

 Writing file:        partition_output/domain_number_48
 Finished writing:    partition_output/domain_number_48

 Partitioning finished (1.54 s)

 ----------------------------------------------------------
 Composing periodicities

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

 Face interfaces creation
 Definition of periodic vertices
 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 =             106823
    maximum value =             198028

      1 : [     106823 ;     125064 [ =          2
      2 : [     125064 ;     143305 [ =          4
      3 : [     143305 ;     161546 [ =         27
      4 : [     161546 ;     179787 [ =         12
      5 : [     179787 ;     198028 ] =          3

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

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

    minimum value =             121971
    maximum value =             222989

      1 : [     121971 ;     142174 [ =          2
      2 : [     142174 ;     162378 [ =          4
      3 : [     162378 ;     182581 [ =         30
      4 : [     182581 ;     202785 [ =          9
      5 : [     202785 ;     222989 ] =          3

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

 Histogram of the number of ghost cells per rank:

    minimum value =              15148
    maximum value =              24961

      1 : [      15148 ;      17110 [ =          3
      2 : [      17110 ;      19073 [ =         24
      3 : [      19073 ;      21035 [ =         13
      4 : [      21035 ;      22998 [ =          6
      5 : [      22998 ;      24961 ] =          2

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

 Histogram of the number of interior faces per rank:

    minimum value =             328957
    maximum value =             608206

      1 : [     328957 ;     384806 [ =          2
      2 : [     384806 ;     440656 [ =          4
      3 : [     440656 ;     496506 [ =         28
      4 : [     496506 ;     552356 [ =         11
      5 : [     552356 ;     608206 ] =          3

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

 Histogram of the number of boundary faces per rank:

    minimum value =                  0
    maximum value =               4051

      1 : [          0 ;        810 [ =         15
      2 : [        810 ;       1620 [ =         13
      3 : [       1620 ;       2430 [ =         13
      4 : [       2430 ;       3240 [ =          4
      5 : [       3240 ;       4051 ] =          3

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

 Histogram of the number of interior faces per cell:

    minimum value =                  4
    maximum value =                  6

      1 : [          4 ;          5 [ =         50
      2 : [          5 ;          6 ] =    7509100

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

 Histogram of the number of neighboring domains per rank:

    minimum value =                  3
    maximum value =                  8

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

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

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

 Renumbering mesh:

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

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

   renumbering for boundary faces:
     numbering:                           none

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

 Histogram of volume mesh matrix bandwidth per rank:

    minimum value =             119866
    maximum value =             217907

      1 : [     119866 ;     139474 [ =          2
      2 : [     139474 ;     159082 [ =          5
      3 : [     159082 ;     178690 [ =         28
      4 : [     178690 ;     198298 [ =          9
      5 : [     198298 ;     217907 ] =          4

 Histogram of volume mesh matrix profile/lines per rank:

    minimum value =               9809
    maximum value =              17547

      1 : [       9809 ;      11356 [ =          5
      2 : [      11356 ;      12904 [ =         20
      3 : [      12904 ;      14451 [ =         16
      4 : [      14451 ;      15999 [ =          6
      5 : [      15999 ;      17547 ] =          1

 Mesh coordinates:               minimum    and maximum
                       X : -1.8299983e+00  3.0500002e+00
                       Y : -2.4400001e+00  2.4400001e+00
                       Z : -3.3825634e-09  6.1000004e-02
 Mesh
     Number of cells:          7509150
     Number of interior faces: 22793391
     Number of boundary faces: 68850
     Number of vertices:       7845318

 Periodic faces (which are also interior faces):
     Periodicity  1:          300366 face couples

 Groups:
    "INLET"
       boundary faces:        37925
    "OUTLET"
       boundary faces:         9900
    "SIDE_1"
       interior faces:       300366
    "SIDE_2"
       interior faces:       300366
    "WALL"
       boundary faces:        21025

 --- Information on the volumes
       Minimum control volume      =  1.7442262e-10
       Maximum control volume      =  5.9230663e-07
       Total volume for the domain =  1.2950347e+00

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

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

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

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

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

 Computing geometric quantities (0.598 s)
                                                             
       ALMAX  =    0.10899E+01 (Characteristic length       )
       ALMAX is the cubic root of the domain volume.

       ALMAX is the length used to initialize the turbulence.

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

   ** READING MAIN AND AUXILIARY RESTART FILES
      ----------------------------------------


      READING THE MAIN RESTART FILE

 Start reading
 Reading dimensions complete
  Reading the previous time step number (restarting computation)  NTPABS =        10923
  Reading the previous time step number (restarting computation)  TTPABS =   0.5534E-02
 Reading options complete
  Read variables from restart: restart/main
 Reading complete

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


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


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

 -----------------------------------------
  Variable          Min. value  Max. value
 -----------------------------------------                   
  Velocity          0.2509E+01  0.3122E+03
  Velocity         -0.9260E+02  0.1363E+03
  Velocity         -0.5854E-02  0.5855E-02
  Pressure          0.7482E+05  0.1393E+06
  nu_tilda          0.1736E-03  0.3145E-01
  TotEner           0.2191E+06  0.2678E+06
  temperature       0.3000E+03  0.3000E+03
 ---------------------------------

                dt  0.5000E-06  0.5000E-06
 ---------------------------------



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

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



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



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


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




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


 INSTANT    0.553450000E-02   TIME STEP NUMBER           10924
 ============================================================= 


 -----------------------------------------
 Property           Min. value  Max. value
 -----------------------------------------
  density           0.8688E+00  0.1618E+01
  molecular_viscos  0.1830E-04  0.1830E-04
  turbulent_viscos  0.1771E-03  0.3694E-01
 -----------------------------------------

 --- Diffusivity:
 -----------------------------------------------
 Scalar           Number  Min. value  Max. value
 -----------------------------------------------
 TotEner               1  0.3417E-04  0.3417E-04
 temperature           2  0.2495E-01  0.2495E-01
 -----------------------------------------------

