command: 
 /gpfs5r/workgpfs/rech/iam/riam642/ETUDES/taycoupoi301/hybride/RESU/20131117-1739/cs_solver --mpi

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

                                  (R)
                      Code_Saturne

                      Version 3.0.1


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

  build Wed Nov 13 21:54:46 2013
  MPI version 2.2 (MPICH2 1.5)


  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:                Sun Nov 17 17:40:17 2013
  System:              CNK 2.6.32-279.14.1.bgq.el6_V1R2M0_36.ppc64
  Machine:             bgqion09-ib
  Processor:           model		: ibm,bluegeneq

  Memory:              16384 MB
  Directory:           /gpfs5r/workgpfs/rech/iam/riam642/ETUDES/taycoupoi301/hybride/RESU/20131117-1739
  MPI ranks:           128
  Hardware threads:    32
  Block shape:         <2,2,4,2,2>
  Torus links enabled: <1,1,1,1,1>
  I/O read method:     collective MPI-IO (explicit offsets)
  I/O write method:    collective MPI-IO (explicit offsets)
  I/O rank step:        1

 Reading metadata from file: "mesh_input"

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

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

@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@
@ @@ ATTENTION : ENTREE DES DONNEES
@    =========
@    REFERENCE VELOCITY  UREF WAS NOT DEFINED
@    or is ill defined  (NEGATIVE value ? ).
@   It IS EQUAL  -0.10000E+14
@
@  calculation can only run if turbulence is defined
@  via a restart file or the interface or the
@  cs_user_initialization subroutine.
@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@

@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@
@ @@   WARNING :      WHILE READING INPUT DATA
@    =========
@    DEFECTS CAN APPEAR WHEN USING A COMBINATION OF
@     ANGULAR PERIODICITE (ROTATION) AND RSTM RIJ-EPSILON.
@
@   COMMANDE_PERIO was defined in the run-case script
@
@    and IPERIO =          1)
@    and some periodic boundaries involve rotation
@      Flag for turb ITURB is =         32
@
@    Job can run.
@
@    The defects are related to the turbulent transport terms
@    in the Re stress equations (equivalent to an anisotropic
@    diffusion tensor), but these terms are generaly small
@
@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@

@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@
@ @@   WARNING :      WHILE READING INPUT DATA
@    =========
@
@   DEFECTS CAN APPEAR WHEN USING A COMBINATION OF
@     ANGULAR PERIODICITE (ROTATION) AND SECOND ORDER SCHEME
@     IN TIME FOR VELOCITY
@
@   COMMANDE_PERIO was defined in the run-case script
@
@    and IPERIO =          1)
@  Flags  THETAV for 3 velocity components      Ux, Uy, Uz
@    of velocity
@    are selected
@    with the following values :
@    THETAV(IU)  THETAV(IV)  THETAV(IW)
@           0.50000E+00        0.50000E+00        0.50000E+00
@
@  Job can run.
@     defects are due to the computation of the velocity
@    gradient tensor at rotation perodicity boundaries
@    This uses an explicit approach which reduces the order
@    of the time discretization scheme
@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@


 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   =              3 (Nb families                 )

 --- Physics
       NVAR   =             12 (Nb variables                )
       NSCAL  =              0 (Nb scalars                  )
       NSCAUS =              0 (Nb user scalars             )
       NSCAPP =              0 (Nb specific physics scalars )
       NPROCE =             27 (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.10000E+01 (Reference density           )
       VISCL0 =    0.26702E-03 (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  =             32 (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    )

   - Rij-epsilon EBRSM     (ITURB = 32)
       ALMAX  =   -0.99900E+03 (Characteristic length       )
       UREF   =   -0.10000E+14 (Characteristic velocity     )
       IRIJNU =              0 (Matrix stabilization        )
       IRIJRB =              0 (Reconstruct at boundaries   )
       IGRARI =              1 (Account for gravity         )
       ICLSYR =              1 (Symmetry implicitation      )
       ICLPTR =              0 (Wall implicitation          )

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

   - EBRSM Rij-epsilon     (ITURB = 32)
       CEBMS1 =    0.17000E+01 (Cs1 coeff.                  )
       CEBMR1 =    0.90000E+00 (Cr1 coeff.                  )
       CEBMR2 =    0.80000E+00 (Cr2 coeff.                  )
       CEBMR3 =    0.65000E+00 (Cr3 coeff.                  )
       CEBMR4 =    0.62500E+00 (Cr4 coeff.                  )
       CEBMR5 =    0.20000E+00 (Cr5 coeff.                  )
       CSEBM  =    0.21000E+00 (Rij Cs diffusion coeff.     )
       CEBMR6 =    0.60000E+00 (Gravity term coeff.         )
       CEBME2 =    0.18300E+01 (Coef Ceps2                  )
       Ce1    =    0.14400E+01 (Coef Ceps1                  )
       SIGEBM =    0.11500E+01 (Coef sigma_eps              )
       XA1    =    0.10000E+00 (Coef A1                     )
       SIGMAK =    0.10000E+01 (Coef sigma_k                )
       XCETA  =    0.80000E+02 (Coef Ceta                   )
       XCT    =    0.60000E+01 (Coef CT                     )


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


 ** 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.14044E-02 (Minimum time step           )
       DTMAX  =    0.14044E+02 (Maximum time step           )
       DTREF  =    0.14044E-01 (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
 R11                   1  0.1000E+01
 R22                   1  0.1000E+01
 R33                   1  0.1000E+01
 R12                   1  0.1000E+01
 R13                   1  0.1000E+01
 R23                   1  0.1000E+01
 Turb Dissipation      1  0.1000E+01
 Alphap                0  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      0 0.10E+01 0.50E+00
 VelocityY             1      1      1      1      0 0.10E+01 0.50E+00
 VelocityZ             1      1      1      1      0 0.10E+01 0.50E+00
 R11                   1      1      1      1      0 0.00E+00 0.50E+00
 R22                   1      1      1      1      0 0.00E+00 0.50E+00
 R33                   1      1      1      1      0 0.00E+00 0.50E+00
 R12                   1      1      1      1      0 0.00E+00 0.50E+00
 R13                   1      1      1      1      0 0.00E+00 0.50E+00
 R23                   1      1      1      1      0 0.00E+00 0.50E+00
 Turb Dissipation      1      1      1      1      0 0.00E+00 0.50E+00
 Alphap                0      1      1      1      0 0.00E+00 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 =              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 =              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
 R11                 100      1  0.1000E-04  0.1000E-06  0.0000E+00
 R22                 100      1  0.1000E-04  0.1000E-06  0.0000E+00
 R33                 100      1  0.1000E-04  0.1000E-06  0.0000E+00
 R12                 100      1  0.1000E-04  0.1000E-06  0.0000E+00
 R13                 100      1  0.1000E-04  0.1000E-06  0.0000E+00
 R23                 100      1  0.1000E-04  0.1000E-06  0.0000E+00
 Turb Dissipation    100      1  0.1000E-04  0.1000E-06  0.0000E+00
 Alphap              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
 R11                   1     -1  0.1500E+01
 R22                   1     -1  0.1500E+01
 R33                   1     -1  0.1500E+01
 R12                   1     -1  0.1500E+01
 R13                   1     -1  0.1500E+01
 R23                   1     -1  0.1500E+01
 Turb Dissipation      1     -1  0.1500E+01
 Alphap                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 =             18 (Number of moments           )

------------------------------------------------------
 IMOM IMOOLD NDTMOM IDFMOM
------------------------------------------------------
    1     -1  35604   0.00000E+00      2      0      0      0      0
    2     -1  35604   0.00000E+00      3      0      0      0      0
    3     -1  35604   0.00000E+00      4      0      0      0      0
    4     -1  35604   0.00000E+00      5      0      0      0      0
    5     -1  35604   0.00000E+00      6      0      0      0      0
    6     -1  35604   0.00000E+00      7      0      0      0      0
    7     -1  35604   0.00000E+00      8      0      0      0      0
    8     -1  35604   0.00000E+00      9      0      0      0      0
    9     -1  35604   0.00000E+00     10      0      0      0      0
   10     -1  35604   0.00000E+00      2      2      0      0      0
   11     -1  35604   0.00000E+00      3      3      0      0      0
   12     -1  35604   0.00000E+00      4      4      0      0      0
   13     -1  35604   0.00000E+00      2      3      0      0      0
   14     -1  35604   0.00000E+00      2      4      0      0      0
   15     -1  35604   0.00000E+00      3      4      0      0      0
   16     -1  35604   0.00000E+00     11      0      0      0      0
   17     -1  35604   0.00000E+00     12      0      0      0      0
   18     -1  35604   0.00000E+00      1      0      0      0      0
------------------------------------------------------

       IMOM   = 0 ou > 0       (moment number               )
       IMOOLD =-1 ou > 0       (old moment corresponding    )
                               (  to calculation restart or )
                               (-1 if the moment is         )
                               (  reinitialized             )
       NTDMOM = 0 ou > 0       (moment calculation starting )
                               (time step number            )
       TTDMOM = 0 ou > 0       (moment calculation starting )
                               (time                        )
       IDFMOM = 0 ou > 0       (number of variables of      )
                               (which the moment is composed)


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


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

--------------------------------------------------
 Variable         IRESOL NITMAX      EPSILO IDIRCL
--------------------------------------------------
 Pressure             -1  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
 R11                  -1  10000  0.1000E-07      1
 R22                  -1  10000  0.1000E-07      1
 R33                  -1  10000  0.1000E-07      1
 R12                  -1  10000  0.1000E-07      1
 R13                  -1  10000  0.1000E-07      1
 R23                  -1  10000  0.1000E-07      1
 Turb Dissipation     -1  10000  0.1000E-07      1
 Alphap               -1  10000  0.1000E-07      0
-----------------------------------

       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
 R11                   0    100     10
 R22                   0    100     10
 R33                   0    100     10
 R12                   0    100     10
 R13                   0    100     10
 R23                   0    100     10
 Turb Dissipation      0    100     10
 Alphap                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 =         142411 (Final time step required    )

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


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


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

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

 --- Post-processed variables

       Number Name
         2 Pressure        
         3 VelocityX       
         4 VelocityY       
         5 VelocityZ       
         6 R11             
         7 R22             
         8 R33             
         9 R12             
        10 R13             
        11 R23             
        12 Turb Dissipation
        13 Alphap          
        14 Density         
        20 MoyU            
        21 MoyV            
        22 MoyW            
        23 MoyuuSGS        
        24 MoyvvSGS        
        25 MoywwSGS        
        26 MoyuvSGS        
        27 MoyuwSGS        
        28 MoyvwSGS        
        29 MoyUU           
        30 MoyVV           
        31 MoyWW           
        32 MoyUV           
        33 MoyUW           
        34 MoyVW           
        35 MoyEps          
        36 MoyAlpha        
        37 MoyP            
         --           --

 --- 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 R11                           -1
         7 R22                           -1
         8 R33                           -1
         9 R12                           -1
        10 R13                           -1
        11 R23                           -1
        12 Turb Dissipation              -1
        13 Alphap                        -1
        20 MoyU                          -1
        21 MoyV                          -1
        22 MoyW                          -1
        23 MoyuuSGS                      -1
        24 MoyvvSGS                      -1
        25 MoywwSGS                      -1
        26 MoyuvSGS                      -1
        27 MoyuwSGS                      -1
        28 MoyvwSGS                      -1
        29 MoyUU                         -1
        30 MoyVV                         -1
        31 MoyWW                         -1
        32 MoyUV                         -1
        33 MoyUW                         -1
        34 MoyVW                         -1
        35 MoyEps                        -1
        36 MoyAlpha                      -1
        37 MoyP                          -1
         --           --                --

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

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

         2 Pressure                       0
         3 VelocityX                      0
         4 VelocityY                      0
         5 VelocityZ                      0
         6 R11                            0
         7 R22                            0
         8 R33                            0
         9 R12                            0
        10 R13                            0
        11 R23                            0
        12 Turb Dissipation               0
        13 Alphap                         0
        20 MoyU                        -999
        21 MoyV                        -999
        22 MoyW                        -999
        23 MoyuuSGS                    -999
        24 MoyvvSGS                    -999
        25 MoywwSGS                    -999
        26 MoyuvSGS                    -999
        27 MoyuwSGS                    -999
        28 MoyvwSGS                    -999
        29 MoyUU                       -999
        30 MoyVV                       -999
        31 MoyWW                       -999
        32 MoyUV                       -999
        33 MoyUW                       -999
        34 MoyVW                       -999
        35 MoyEps                      -999
        36 MoyAlpha                    -999
        37 MoyP                        -999
         --           --                --

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


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


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

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


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


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

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

 Partitioning by space-filling curve: Morton (in bounding box).
  Number of cells per domain (histogramm):
    [       1624 ;       1630 [ =          1
    [       1630 ;       1636 [ =          0
    [       1636 ;       1643 [ =          0
    [       1643 ;       1649 [ =          0
    [       1649 ;       1656 [ =          0
    [       1656 ;       1662 [ =          0
    [       1662 ;       1668 [ =          0
    [       1668 ;       1675 [ =          0
    [       1675 ;       1681 [ =          0
    [       1681 ;       1688 ] =        127

 Writing file:        partition_output/domain_number_128
 Finished writing:    partition_output/domain_number_128

 Partitioning finished (8.84 s)


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

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


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

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

  Selection criteria: "Z < 0.01 or Z > 9.99"

  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:          216000
     Number of interior faces: 637200
     Number of boundary faces: 21600
     Number of vertices:       226981

  Global number of boundary faces selected for joining:       7200

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

 Glob. Num. |  Tolerance  |              Coordinates

          3 |    0.000304 |  9.0000000000e+00   0.0000000000e+00   0.0000000000e+00 | ORI
       2571 |    0.004822 |  6.6057437479e+00   5.3492195260e+00   0.0000000000e+00 | ORI
  Determination of possible face intersections:

    bounding-box tree layout: 3D

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

  Joining operation is conforming.

  Edge intersections and vertex creation done.

  Pre-merge for 26009 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:      10800

  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:          216000
     Number of interior faces: 644400
     Number of boundary faces: 14400
     Number of vertices:       226981


  Joining  1 completed (0.724 s).

 -------------------------------------------------------
  Joining number 2:

  Periodicity type: rotation
  Transformation matrix:             0           -1            0            0
                                     1            0            0            0
                                     0            0            1            0

  Selection criteria: "X < 0.001 or Y < 0.001"

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

  Advanced joining parameters:
    Verbosity level:                                 2
    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:          216000
     Number of interior faces: 644400
     Number of boundary faces: 14400
     Number of vertices:       226981

  Global number of boundary faces selected for joining:       7200

  Element selection successfully done.
  Global number of selected vertices:        7442

  Global min/max. tolerance:

 Glob. Num. |  Tolerance  |              Coordinates

          3 |    0.000608 |  9.0000000000e+00   0.0000000000e+00   1.7763568394e-15 | ORI
          2 |    0.009644 |  8.5000000000e+00   0.0000000000e+00   0.0000000000e+00 | ORI
  Determination of possible face intersections:

    bounding-box tree layout: 3D

  Sorted possible intersections between faces.

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

  Joining operation is conforming.

  Edge intersections and vertex creation done.

  Pre-merge for 14702 global element couples.

  No need to merge vertices.

  Merging of equivalent vertices done.

  Number of simplified faces: 0

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

  Number of simplified faces: 0

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

  Merge vertices and mesh update done.

  Global number of faces after splitting:      10800

  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:          216000
     Number of interior faces: 651600
     Number of boundary faces: 7200
     Number of vertices:       226981


  Joining  2 completed (0.436 s).

  All joining operations successfully finished:

    Wall clock time:                  1.16


 ----------------------------------------------------------
 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 =               1624
    maximum value =               1688

      1 : [       1624 ;       1636 [ =          1
      2 : [       1636 ;       1649 [ =          0
      3 : [       1649 ;       1662 [ =          0
      4 : [       1662 ;       1675 [ =          0
      5 : [       1675 ;       1688 ] =        127

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

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

    minimum value =               2819
    maximum value =               4148

      1 : [       2819 ;       3084 [ =         44
      2 : [       3084 ;       3350 [ =         59
      3 : [       3350 ;       3616 [ =         21
      4 : [       3616 ;       3882 [ =          0
      5 : [       3882 ;       4148 ] =          4

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

 Histogram of the number of ghost cells per rank:

    minimum value =               1131
    maximum value =               2460

      1 : [       1131 ;       1396 [ =         44
      2 : [       1396 ;       1662 [ =         59
      3 : [       1662 ;       1928 [ =         21
      4 : [       1928 ;       2194 [ =          0
      5 : [       2194 ;       2460 ] =          4

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

 Histogram of the number of interior faces per rank:

    minimum value =               5478
    maximum value =               6274

      1 : [       5478 ;       5637 [ =          1
      2 : [       5637 ;       5796 [ =         39
      3 : [       5796 ;       5955 [ =         71
      4 : [       5955 ;       6114 [ =         13
      5 : [       6114 ;       6274 ] =          4

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

 Histogram of the number of boundary faces per rank:

    minimum value =                 26
    maximum value =                146

      1 : [         26 ;         50 [ =         46
      2 : [         50 ;         74 [ =         66
      3 : [         74 ;         98 [ =         10
      4 : [         98 ;        122 [ =          4
      5 : [        122 ;        146 ] =          2

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

 Histogram of the number of interior faces per cell:

    minimum value =                  5
    maximum value =                  6

      1 : [          5 ;          6 ] =     216000

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

 Histogram of the number of neighboring domains per rank:

    minimum value =                 11
    maximum value =                 34

      1 : [         11 ;         15 [ =         24
      2 : [         15 ;         20 [ =         81
      3 : [         20 ;         24 [ =         12
      4 : [         24 ;         29 [ =          7
      5 : [         29 ;         34 ] =          4

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

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

 Writing file:        mesh_output
 Finished writing:    mesh_output

 Renumbering mesh:

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


 Histogram of volume mesh matrix bandwidth per rank:

    minimum value =               2097
    maximum value =               3974

      1 : [       2097 ;       2472 [ =          9
      2 : [       2472 ;       2847 [ =         55
      3 : [       2847 ;       3223 [ =         52
      4 : [       3223 ;       3598 [ =          8
      5 : [       3598 ;       3974 ] =          4

 Histogram of volume mesh matrix profile/lines per rank:

    minimum value =                774
    maximum value =               2045

      1 : [        774 ;       1028 [ =         31
      2 : [       1028 ;       1282 [ =         67
      3 : [       1282 ;       1536 [ =         22
      4 : [       1536 ;       1790 [ =          6
      5 : [       1790 ;       2045 ] =          2

 Mesh coordinates:               minimum    and maximum
                       X :  5.3071796e-06  9.0000000e+00
                       Y : -5.9705767e-06  9.0000000e+00
                       Z : -1.7763568e-15  1.0000000e+01
 Mesh
     Number of cells:          216000
     Number of interior faces: 651600
     Number of boundary faces: 7200
     Number of vertices:       226981

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

 Groups:
    "7"
       cells:                216000
       interior faces:        14400
       boundary faces:         7200

 --- Information on the volumes
       Minimum control volume      =  1.0611085e-04
       Maximum control volume      =  1.7932265e-03
       Total volume for the domain =  1.3350244e+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: 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 %

  Criteria 6: User Specific Tag:
    Number of bad cells detected: 216000 --> 100 %
 Warning:
 --------
    Mesh quality issue based on user criteria has been detected

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


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

       ALMAX is the length used to initialize the turbulence.

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


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

 -----------------------------------------
  Variable          Min. value  Max. value
 -----------------------------------------                   
  Pressure          0.0000E+00  0.0000E+00
  VelocityX        -0.4386E+01 -0.1153E-02
  VelocityY        -0.4386E+01  0.5021E+01
  VelocityZ        -0.9465E-01  0.2737E+00
  R11               0.1407E-05  0.1157E+00
  R22               0.1704E-05  0.1168E+00
  R33               0.7581E-03  0.5550E-01
  R12              -0.1050E-01  0.5644E-01
  R13              -0.4915E-01  0.4928E-01
  R23              -0.4941E-01  0.4954E-01
  Turb Dissipation  0.1605E-01  0.2313E+00
  Alphap            0.4825E-01  0.9994E+00
 ---------------------------------

  MoyU              0.0000E+00  0.0000E+00
  MoyV              0.0000E+00  0.0000E+00
  MoyW              0.0000E+00  0.0000E+00
  MoyuuSGS          0.0000E+00  0.0000E+00
  MoyvvSGS          0.0000E+00  0.0000E+00
  MoywwSGS          0.0000E+00  0.0000E+00
  MoyuvSGS          0.0000E+00  0.0000E+00
  MoyuwSGS          0.0000E+00  0.0000E+00
  MoyvwSGS          0.0000E+00  0.0000E+00
  MoyUU             0.0000E+00  0.0000E+00
  MoyVV             0.0000E+00  0.0000E+00
  MoyWW             0.0000E+00  0.0000E+00
  MoyUV             0.0000E+00  0.0000E+00
  MoyUW             0.0000E+00  0.0000E+00
  MoyVW             0.0000E+00  0.0000E+00
  MoyEps            0.0000E+00  0.0000E+00
  MoyAlpha          0.0000E+00  0.0000E+00
  MoyP              0.0000E+00  0.0000E+00
 ---------------------------------

  Local Time Step   0.1404E-01  0.1404E-01
 ---------------------------------

 Time averages (sum over the time-steps)
 ------------------------------------------------------------
   Average  Min. value  Max. value Uniform/Variable in space
 ------------------------------------------------------------
         1  0.0000E+00  0.0000E+00  Uniform
         2  0.0000E+00  0.0000E+00  Uniform
         3  0.0000E+00  0.0000E+00  Uniform
         4  0.0000E+00  0.0000E+00  Uniform
         5  0.0000E+00  0.0000E+00  Uniform
         6  0.0000E+00  0.0000E+00  Uniform
         7  0.0000E+00  0.0000E+00  Uniform
         8  0.0000E+00  0.0000E+00  Uniform
         9  0.0000E+00  0.0000E+00  Uniform
        10  0.0000E+00  0.0000E+00  Uniform
        11  0.0000E+00  0.0000E+00  Uniform
        12  0.0000E+00  0.0000E+00  Uniform
        13  0.0000E+00  0.0000E+00  Uniform
        14  0.0000E+00  0.0000E+00  Uniform
        15  0.0000E+00  0.0000E+00  Uniform
        16  0.0000E+00  0.0000E+00  Uniform
        17  0.0000E+00  0.0000E+00  Uniform
        18  0.0000E+00  0.0000E+00  Uniform
 ------------------------------------------------------------



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

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



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



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


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




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


 INSTANT    0.140438729E-01   TIME STEP NUMBER               1
 ============================================================= 


 -----------------------------------------
 Property           Min. value  Max. value
 -----------------------------------------
  Density           0.1000E+01  0.1000E+01
  Laminar Viscosit  0.2670E-03  0.2670E-03
  Turb Viscosity    0.6384E-06  0.4265E-02
 -----------------------------------------



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

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



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

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


   ** BOUNDARY CONDITIONS FOR SMOOTH WALLS
   ---------------------------------------
------------------------------------------------------------
                                         Minimum     Maximum
------------------------------------------------------------
   Rel velocity at the wall uiptn :  0.00000E+00 0.00000E+00
   Friction velocity        uet   :  0.17798E-01 0.56157E+00
   Friction velocity        uk    :  0.17798E-01 0.56157E+00
   Dimensionless distance   yplus :  0.10127E+00 0.31966E+01
   ------------------------------------------------------
   Nb of reversal of the velocity at the wall   :          0
   Nb of faces within the viscous sub-layer     :       3600
   Total number of wall faces                   :       7200
------------------------------------------------------------


Signal 5 intercepted!
SIGABRT signal (Abort) intercepted.
