command: 
 /home/butsb/Desktop/HEATED-RECTANGULAR-CAVITY/PURE-CONVECTION/RESU/20161123-1106/cs_solver --param PURE-CONVECTION.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:                Wed 23 Nov 2016 11:06:24 AM PKT
  System:              Linux 3.2.0-60-generic
  Machine:             butsb-PowerEdge-R815
  Processor:           model name	: AMD Opteron(tm) Processor 6348
  Memory:              257938 MB
  User:                butsb (butsb)
  Directory:           /home/butsb/Desktop/HEATED-RECTANGULAR-CAVITY/PURE-CONVECTION/RESU/20161123-1106
  MPI ranks:           12 (appnum attribute: 0)
  I/O read method:     collective MPI-IO (explicit offsets)
  I/O write method:    collective MPI-IO (explicit offsets)
  I/O rank step:        1

  External libraries for partitioning:
    ParMETIS 3.1
    PT-SCOTCH 1.12.2147483631

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

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

 --- Physics
       NVAR   =              5 (Nb variables                )
       NSCAL  =              1 (Nb scalars                  )
       NSCAUS =              1 (Nb user scalars             )
       NSCAPP =              0 (Nb specific physics scalars )
       NPROCE =              6 (Nb cell properties          )
       NPROFA =              1 (Nb internal face properties )
       NPROFB =              2 (Nb boundary face properties )


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


 ** PHYSICAL PROPERTIES
    -------------------

       GX     =    0.00000E+00 (Gravity x component         )
       GY     =   -0.98100E+01 (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.11650E+01 (Reference density           )
       VISCL0 =    0.18500E-04 (Ref. molecular dyn. visc.   )
       CP0    =    0.10090E+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 =              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

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


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

 --- Continuous phase:

   - Commons
       ITURB  =              0 (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.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)     )
       ISCALT =              1 (Temperature salar number    )

 --- 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 =              2 (0 cst; 1,2 var (t, t-space  )
       IPTLRO =              0 (1: rho-related DT clipping  )
       COUMAX =    0.80000E+01 (Maximum target CFL          )
       FOUMAX =    0.10000E+02 (Maximum target Fourier      )
       VARRDT =    0.10000E+00 (For var. DT, max. increase  )
       DTMIN  =    0.10000E-01 (Minimum time step           )
       DTMAX  =    0.10000E+03 (Maximum time step           )
       DTREF  =    0.10000E+00 (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
 TempK                 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      1 0.10E+01 0.10E+01
 VelocityY             1      1      1      1      1 0.10E+01 0.10E+01
 VelocityZ             1      1      1      1      1 0.10E+01 0.10E+01
 TempK                 1      1      1      1      1 0.10E+01 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
    ------
       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 =              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
 TempK               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
 TempK                 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
 TempK                -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
 TempK                 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

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


 ** SCALARS
    -------

       ITBRRB =              0 (T or H reconstruction at bdy)

---------------------------------------------------------------------
 Variable         Number ISCSTH IVISLS  ITURT      VISLS0      SIGMAS
---------------------------------------------------------------------
 TempK                 1      1      0      0  0.2660E-01  0.1000E+01
---------------------------------------------------------------------

-------------------------------------------
 Variable         Number ISCAVR      RVARFL
-------------------------------------------
 TempK                 1      0  0.8000E+00
-------------------------------------------

-------------------------------------------------------
 Variable         Number ICLVFL      SCAMIN      SCAMAX
-------------------------------------------------------
 TempK                 1     -1  0.0000E+00  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.

       ISCSTH = -1,0, 1 ou 2   (T (C), Passive, T (K) or H  )
       IVISLS = 0 ou >0        (Viscosity: constant or not  )
       VISLS0 = >0             (Reference viscosity         )
       SIGMAS = >0             (Schmidt                     )
       ISCAVR = 0 ou >0        (Associat. scalar if variance)
       RVARFL = >0             (Rf, cf variance dissipation )
       ICLVFL = 0, 1 ou 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
------------------------------------------------------

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


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


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

 --- Restarted calculation
       ISUITE =              0 (1: restarted calculation    )
       ILEAUX =              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 =           1500 (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 TempK           
         7 Density         
         9 CourantNb       
        10 FourierNb       
        11 total_pressure  
        12 LocalTime       
         --           --

 --- 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 TempK                         -1
         7 Density                       -1
         8 LamVisc                       -1
         9 CourantNb                     -1
        10 FourierNb                     -1
        11 total_pressure                -1
        12 LocalTime                     -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 TempK                          0
         7 Density                     -999
         9 CourantNb                   -999
        10 FourierNb                   -999
        11 total_pressure              -999
        12 LocalTime                   -999
         --           --                --

 --- Additional post-processing variables (ipstdv)
       ipstfo =              0 (Force exerted  at boundary)
       ipstyp =              0 (Yplus          at boundary)
       ipstft =              1 (Thermal flow   at boundary)
       ipsttb =              1 (Temperature    at boundary)
       ipstnu =              1 (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_12" file available;

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

 Partitioning by space-filling curve: Morton (in bounding box).
  Number of cells per domain (histogramm):
    [      12994 ;      12995 [ =          1
    [      12995 ;      12996 [ =          0
    [      12996 ;      12997 [ =          0
    [      12997 ;      12998 [ =          0
    [      12998 ;      12999 [ =          0
    [      12999 ;      13000 [ =          0
    [      13000 ;      13001 [ =          0
    [      13001 ;      13002 ] =         11

 Partitioning finished (1.45 s)


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

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


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

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

  Selection criteria: ""L" OR "R""

  Parameters for the joining operation:
    Shortest incident edge fraction:           0.10000
    Maximum angle between joined face planes: 25.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:          156016
     Number of interior faces: 311049
     Number of boundary faces: 313998
     Number of vertices:       314000

  Global number of boundary faces selected for joining:        398

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

 Glob. Num. |  Tolerance  |              Coordinates

          1 |    0.000000 | -5.0000000000e-01   0.0000000000e+00   0.0000000000e+00 | ORI
       1674 |    0.000103 | -5.0000000000e-01   2.6016054923e-02   0.0000000000e+00 | ORI
  Determination of possible face intersections:

    bounding-box tree layout: 2D

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

  Joining operation is conforming.

  Edge intersections and vertex creation done.

  Pre-merge for 866 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:        597

  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:          156016
     Number of interior faces: 311447
     Number of boundary faces: 313600
     Number of vertices:       314000


  Joining  1 completed (0.0473 s).

  All joining operations successfully finished:

    Wall clock time:                0.0474


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

 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 =              12994
    maximum value =              13002

      1 : [      12994 ;      12995 [ =          1
      2 : [      12995 ;      12997 [ =          0
      3 : [      12997 ;      12998 [ =          0
      4 : [      12998 ;      13000 [ =          0
      5 : [      13000 ;      13002 ] =         11

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

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

    minimum value =              13326
    maximum value =              13647

      1 : [      13326 ;      13390 [ =          4
      2 : [      13390 ;      13454 [ =          0
      3 : [      13454 ;      13518 [ =          2
      4 : [      13518 ;      13582 [ =          0
      5 : [      13582 ;      13647 ] =          6

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

 Histogram of the number of ghost cells per rank:

    minimum value =                332
    maximum value =                645

      1 : [        332 ;        394 [ =          4
      2 : [        394 ;        457 [ =          0
      3 : [        457 ;        519 [ =          2
      4 : [        519 ;        582 [ =          0
      5 : [        582 ;        645 ] =          6

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

 Histogram of the number of interior faces per rank:

    minimum value =              26074
    maximum value =              26269

      1 : [      26074 ;      26113 [ =          4
      2 : [      26113 ;      26152 [ =          0
      3 : [      26152 ;      26191 [ =          0
      4 : [      26191 ;      26230 [ =          2
      5 : [      26230 ;      26269 ] =          6

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

 Histogram of the number of boundary faces per rank:

    minimum value =              26102
    maximum value =              26170

      1 : [      26102 ;      26115 [ =          4
      2 : [      26115 ;      26129 [ =          0
      3 : [      26129 ;      26142 [ =          4
      4 : [      26142 ;      26156 [ =          1
      5 : [      26156 ;      26170 ] =          3

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

 Histogram of the number of interior faces per cell:

    minimum value =                  3
    maximum value =                  4

      1 : [          3 ;          4 ] =     156016

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

 Histogram of the number of neighboring domains per rank:

    minimum value =                  2
    maximum value =                  7

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

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

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

 Writing file:        mesh_output
