command: 
 /home/fafan/Code_Saturne/4.3.0/code_saturne-4.3.0/arch/Linux_x86_64/libexec/code_saturne/cs_solver --param pipe.xml --mpi

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

                                  (R)
                      Code_Saturne

                      Version 4.3.0


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

  revision 4.3.0
  build Fri 16 Sep 2016 09:23:06 PM EDT
  MPI version 3.0 (MVAPICH2 2.0, based on MPICH 3.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:                Thu 29 Sep 2016 11:01:54 PM EDT
  System:              Linux 2.6.32-504.16.2.el6.x86_64
  Machine:             compute-0-3.local
  Processor:           model name	: Intel(R) Xeon(R) CPU E5-2640 0 @ 2.50GHz
  Memory:              96859 MB
  User:                fafan (Fa-Gung Fan - fafan)
  Directory:           /home/fafan/Code_Saturne/Saturne/pipe/pipe_v4.3.0/CASE1/RESU/20160927-0059
  MPI ranks:           12 (appnum attribute: 0)
  OpenMP threads:      12
  Processors/node:     1
  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
    PT-SCOTCH 6.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   =              6 (Nb variables                )
       NSCAL  =              0 (Nb scalars                  )
       NSCAUS =              0 (Nb user scalars             )
       NSCAPP =              0 (Nb specific physics scalars )
       NPROCE =              6 (Nb cell properties          )


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


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

       GX     =    0.00000E+00 (Gravity x component         )
       GY     =    0.00000E+00 (Gravity y component         )
       GZ     =    0.00000E+00 (Gravity z component         )

       ICORIO =          0 (Coriolis source terms       )

  -- Continuous phase:

       RO0    =    0.11786E+01 (Reference density           )
       VISCL0 =    0.18300E-04 (Ref. molecular dyn. visc.   )
       CP0    =    0.10172E+04 (Ref. specific heat          )
       ICP    =              0 (> 0: variable CP (usphyv)   )
       P0     =    0.10132E+06 (Ref. total pressure         )
       PRED0  =    0.00000E+00 (Ref. reduced pressure       )
       T0     =    0.29315E+03 (Ref. temperature            )

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

       Initial reference point for pressure
       XYZP0  =    0.00000E+00   0.00000E+00   0.00000E+00

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


 ** 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 =              0 (0: no thermal model         )
                               (1: temperature              )
                               (2: enthalpy                 )
                               (3: total energy             )
       ITPSCL =              0 (0: none                     )
                               (1: temperature in Kelvin    )
                               (2: temperature in Celsius   )
       ISCALT =             -1 (Thermal scalar number       )


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

 --- Continuous phase:

   - Commons
       ITURB  =             21 (Turbulence model            )
       IWALLF =              3 (wall function               )
                               (0: disabled                 )
                               (1: one scale power law
                                   (forbidden for k-epsilon))
                               (2: one scale log law        )
                               (3: two scales log law       )
                               (4: scalable wall function   )
                               (5: two scales V. Driest     )
                               (6: two scales smooth/rough  )
       IWALLT =              0 (Exch. coeff. correlation    )
                               (0: not activated            )
                               (1: activated                )
       YPLULI =    0.23810E+01 (Limit Y+                    )
       IGRHOK =              0 (1: computed Grad(rho k)     )

   - Linear production k-epsilon (ITURB = 21)
       ALMAX  =   -0.99900E+03 (Characteristic length       )
       UREF   =    0.10000E+01 (Characteristic velocity     )
       ICLKEP =              0 (k-epsilon clipping model    )
       IKECOU =              0 (k-epsilon coupling mode     )
       IGRAKE =              1 (Account for gravity         )
       RELAXV =    0.10000E+01 for k        (Relaxation)
       RELAXV =    0.10000E+01 for epsilon  (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))

   - Linear production k-epsilon (ITURB = 21)
       Ce1    =    0.14400E+01 (Cepsilon 1: production coef.)
       CE2    =    0.19200E+01 (Cepsilon 2: dissipat.  coef.)
       SIGMAK =    0.10000E+01 (Prandtl relative to k       )
       SIGMAE =    0.13000E+01 (Prandtl relative to epsilon )
       CMU    =    0.90000E-01 (Cmu constant                )


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


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

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


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


 ** TIME STEPPING
    -------------

    UNSTEADY ALGORITHM

 --- Time step parameters

       IDTVAR =              2 (0 cst; 1,2 var (t, t-space  )
       IPTLRO =              0 (1: rho-related DT clipping  )
       COUMAX =    0.10000E+01 (Maximum target CFL          )
       FOUMAX =    0.10000E+02 (Maximum target Fourier      )
       VARRDT =    0.10000E+00 (For var. DT, max. increase  )
       DTMIN  =    0.10000E-03 (Minimum time step           )
       DTMAX  =    0.10000E+01 (Maximum time step           )
       DTREF  =    0.10000E-02 (Reference time step         )

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

 --- Frozen velocity field

       ICCVFG =              0 (1: frozen velocity field    )

 --- Per-variable properties

------------------------------------
 Variable          ISTAT      CDTVAR
------------------------------------
 Velocity              1  0.1000E+01
 Pressure              0  0.1000E+01
 k                     1  0.1000E+01
 epsilon               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      1      0 0.10E+01 0.10E+01
 Pressure              0      1      1      1      0 0.00E+00 0.10E+01
 k                     1      1      1      1      0 0.00E+00 0.10E+01
 epsilon               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,1,2,3       (0: slope test
                                1: no slope test
                                2: min/max limiter
                                3: Roe Sweby limiter        )
       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 =              0 (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-06  0.0000E+00
 Pressure            100      2  0.1000E-04  0.1000E-06  0.0000E+00
 k                   100      1  0.1000E-04  0.1000E-06  0.0000E+00
 epsilon             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
 k                     1     -1  0.1500E+01
 epsilon               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                )


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


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

------------------------------------
 Variable              EPSILO IDIRCL
------------------------------------
 Velocity          0.1000E-07      1
 Pressure          0.1000E-07      1
 k                 0.1000E-07      1
 epsilon           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    )


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


 ** 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 =            300 (Final time step required    )

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


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


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

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

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

       Number Name                  Nb. probes (-1: all)
         2 Velocity[X]                   -1
         3 Velocity[Y]                   -1
         4 Velocity[Z]                   -1
         5 Pressure                      -1
         6 k                             -1
         7 epsilon                       -1
         9 LamVisc                       -1
        10 TurbVisc                      -1
        11 CourantNb                     -1
        12 FourierNb                     -1
        13 total_pressure                -1
        14 Local Time Step               -1
         --           --                --

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

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

         2 Velocity                       0
         5 Pressure                       0
         6 k                              0
         7 epsilon                        0
        10 TurbVisc                    -999
        11 CourantNb                   -999
        12 FourierNb                   -999
        13 total_pressure              -999
        14 Local Time Step             -999
         --           --                --

 --- Additional post-processing variables (ipstdv)
       ipstfo =              1 (Force exerted by the
                                     fluid on the boundary)
       ipstyp =              1 (y+ at boundary)
       ipsttp =              1 (T+ at boundary)
       ipstft =              1 (Thermal flux   at boundary)
       ipstnu =              0 (Dimensionless thermal
                                          flux at boundary)


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


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

       IALE   =              0 (1: activated                )
       NALINF =              0 (Fluid initialization
                                                  iterations)
       IFLXMW =              0 (ALE mass flux computation
                                0: thanks to vertices
                                1: thanks to mesh velocity)


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


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

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

  -3: name: 
      directory: monitoring
      format: time_plot
      options: 
      time dependency: fixed mesh
      output: every 5 time steps and at calculation end


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

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

 Partitioning 66704 cells to 12 domains on 12 ranks
   (SCOTCH_dgraphPart).

  wall-clock time: 0.215577 s

  Number of cells per domain (histogramm):
    [       5339 ;       5381 [ =          1
    [       5381 ;       5423 [ =          0
    [       5423 ;       5466 [ =          1
    [       5466 ;       5508 [ =          2
    [       5508 ;       5551 [ =          2
    [       5551 ;       5593 [ =          2
    [       5593 ;       5635 [ =          1
    [       5635 ;       5678 [ =          0
    [       5678 ;       5720 [ =          2
    [       5720 ;       5763 ] =          1

 Writing file:        partition_output/domain_number_12
 Finished writing:    partition_output/domain_number_12

 Partitioning finished (0.255 s)


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

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

 Face interfaces creation
 Vertex interfaces creation
 Halo creation
 Halo definition
    Local halo definition
    Distant halo creation
    Updating the faces -> cells connectivity

 Histogram of the number of cells per rank:

    minimum value =               5339
    maximum value =               5763

      1 : [       5339 ;       5423 [ =          1
      2 : [       5423 ;       5508 [ =          3
      3 : [       5508 ;       5593 [ =          4
      4 : [       5593 ;       5678 [ =          1
      5 : [       5678 ;       5763 ] =          3

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

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

    minimum value =               5606
    maximum value =               6272

      1 : [       5606 ;       5739 [ =          2
      2 : [       5739 ;       5872 [ =          1
      3 : [       5872 ;       6005 [ =          2
      4 : [       6005 ;       6138 [ =          4
      5 : [       6138 ;       6272 ] =          3

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

 Histogram of the number of ghost cells per rank:

    minimum value =                244
    maximum value =                574

      1 : [        244 ;        310 [ =          3
      2 : [        310 ;        376 [ =          0
      3 : [        376 ;        442 [ =          0
      4 : [        442 ;        508 [ =          3
      5 : [        508 ;        574 ] =          6

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

 Histogram of the number of interior faces per rank:

    minimum value =              10107
    maximum value =              11395

      1 : [      10107 ;      10364 [ =          1
      2 : [      10364 ;      10622 [ =          1
      3 : [      10622 ;      10879 [ =          3
      4 : [      10879 ;      11137 [ =          4
      5 : [      11137 ;      11395 ] =          3

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

 Histogram of the number of boundary faces per rank:

    minimum value =                743
    maximum value =               1409

      1 : [        743 ;        876 [ =          8
      2 : [        876 ;       1009 [ =          1
      3 : [       1009 ;       1142 [ =          0
      4 : [       1142 ;       1275 [ =          2
      5 : [       1275 ;       1409 ] =          1

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

 Histogram of the number of interior faces per cell:

    minimum value =                  2
    maximum value =                  4

      1 : [          2 ;          3 [ =        108
      2 : [          3 ;          4 ] =      66596

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

 Histogram of the number of neighboring domains per rank:

    minimum value =                  1
    maximum value =                  3

      1 : [          1 ;          2 [ =          2
      2 : [          2 ;          3 ] =         10

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

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

   renumbering for boundary faces:
     numbering:                           none

 Histogram of thread pools size for interior faces per rank:

    minimum value =                 12
    maximum value =                 12

      1 : [         12 ;         12 ] =         12

 Histogram of threading groups count for interior faces per rank:

    minimum value =                  8
    maximum value =                  8

      1 : [          8 ;          8 ] =         12

 Histogram of thread imbalance for interior faces per rank:

    minimum value =         5.12489e-01
    maximum value =         6.78004e-01

      1 : [ 5.12489e-01 ; 6.78004e-01 ] =         12

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

 Histogram of volume mesh matrix bandwidth per rank:

    minimum value =               5533
    maximum value =               6202

      1 : [       5533 ;       5666 [ =          2
      2 : [       5666 ;       5800 [ =          1
      3 : [       5800 ;       5934 [ =          2
      4 : [       5934 ;       6068 [ =          4
      5 : [       6068 ;       6202 ] =          3

 Histogram of volume mesh matrix profile/lines per rank:

    minimum value =               2826
    maximum value =               3132

      1 : [       2826 ;       2887 [ =          2
      2 : [       2887 ;       2948 [ =          3
      3 : [       2948 ;       3009 [ =          2
      4 : [       3009 ;       3070 [ =          2
      5 : [       3070 ;       3132 ] =          3

 Mesh coordinates:               minimum    and maximum
                       X :  0.0000000e+00  1.1000000e+00
                       Y : -9.9999727e-02  6.0000000e-01
                       Z : -1.0000000e-01  1.0000000e-01
 Mesh
     Number of cells:          66704
     Number of interior faces: 127906
     Number of boundary faces: 11004
     Number of vertices:       14142

 Groups:
    "F1"
       boundary faces:          309
    "F2"
       boundary faces:           80
    "F3"
       boundary faces:          451
    "Walls"
       boundary faces:        10164

 --- Information on the volumes
       Minimum control volume      =  1.3175985e-07
       Maximum control volume      =  2.3125474e-06
       Total volume for the domain =  4.2498043e-02

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

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

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

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

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

 Warning:
 --------
    Mesh quality issue has been detected

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

    The calculation will run but the solution quality may be degraded...

 Computing geometric quantities (2.3 s)
                                                             
       ALMAX  =    0.34897E+00 (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
 -----------------------------------------                   
  Velocity          0.0000E+00  0.0000E+00
  Velocity          0.0000E+00  0.0000E+00
  Velocity          0.0000E+00  0.0000E+00
  Pressure          0.0000E+00  0.0000E+00
  k                 0.6000E-03  0.6000E-03
  epsilon           0.3790E-05  0.3790E-05
 ---------------------------------

                dt  0.1000E-02  0.1000E-02
 ---------------------------------



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

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



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



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


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




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


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


 -----------------------------------------
 Property           Min. value  Max. value
 -----------------------------------------
  density           0.1179E+01  0.1179E+01
  molecular_viscos  0.1830E-04  0.1830E-04
  turbulent_viscos  0.1007E-01  0.1007E-01
 -----------------------------------------



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

-------------------------------------------------------------------------
Boundary type          Code    Nb faces
-------------------------------------------------------------------------
Inlet                     2         389
Smooth wall               5       10164
Rough wall                6           0
Symmetry                  4           0
Free outlet               3         451
Free inlet               14           0
Convective inlet         16           0
Free surface             15           0
Undefined                 1           0
-------------------------------------------------------------------------


Boundary faces with free inlet/outlet detected
Update of reference point for total pressure
 XYZP0 =    0.70403E+00   0.60000E+00  -0.30520E-02



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

---------------------------------------------------------------
Boundary type          Code    Nb faces           Mass flow
---------------------------------------------------------------
Inlet                     2         389         0.000000000E+00
Smooth wall               5       10164         0.000000000E+00
Rough wall                6           0         0.000000000E+00
Symmetry                  4           0         0.000000000E+00
Free outlet               3         451         0.000000000E+00
Free inlet               14           0         0.000000000E+00
Convective inlet         16           0         0.000000000E+00
Free surface             15           0         0.000000000E+00
Undefined                 1           0         0.000000000E+00
---------------------------------------------------------------









   ** BOUNDARY CONDITIONS FOR SMOOTH WALLS
 INSTANT    0.300000000E+00   TIME STEP NUMBER             300
 ============================================================= 




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

---------------------------------------------------------------
Boundary type          Code    Nb faces           Mass flow
---------------------------------------------------------------
Inlet                     2         389        -0.202798560E+00
Smooth wall               5       10164        -0.191441957E-18
Rough wall                6           0         0.000000000E+00
Symmetry                  4           0         0.000000000E+00
Free outlet               3         451         0.202798560E+00
Free inlet               14           0         0.000000000E+00
Convective inlet         16           0         0.000000000E+00
Free surface             15           0         0.000000000E+00
Undefined                 1           0         0.000000000E+00
---------------------------------------------------------------


   ** BOUNDARY CONDITIONS FOR SMOOTH WALLS
   ---------------------------------------
------------------------------------------------------------
                                         Minimum     Maximum
------------------------------------------------------------
   Rel velocity at the wall uiptn :  0.27565E-01 0.58680E+01
   Friction velocity        uet   :  0.30197E-02 0.58235E+00
   Friction velocity        uk    :  0.51474E-01 0.54217E+00
   Dimensionless distance   yplus :  0.10536E+02 0.13486E+03
   ------------------------------------------------------
   Nb of reversal of the velocity at the wall   :          0
   Nb of faces within the viscous sub-layer     :          0
   Total number of wall faces                   :      10164
------------------------------------------------------------


@                                                            
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@                                                            
@ @@ WARNING: MESH TOO REFINED AT THE WALL                   
@    ========                                                
@    The mesh seems to be too refined at the wall to use     
@      a high-Reynolds turbulence model.                     
@                                                            
@    The last time step at which too small values for the    
@      dimensionless distance to the wall (yplus) have been  
@      observed is the time step          1
@                                                            
@    The minimum value for yplus must be greater than the    
@      limit value YPLULI =    0.23810E+01
@                                                            
@    Have a look at the distribution of yplus at the wall    
@      (with EnSight or ParaView for example) to conclude on 
@      the way the results quality might be affected.        
@                                                            
@    This warning is only printed at the first two           
@      occurences of the problem and at the last time step   
@      of the calculation. The vanishing of the message does 
@      not necessarily mean the vanishing of the problem.    
@                                                            
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@                                                            



 Write final restart files
   checkpoint at iteration        300, Physical time    0.30000E+00


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

 Start writing
  Wrote field names and types to checkpoint: checkpoint/main
   End writing the options                              
  Wrote main variables to checkpoint: checkpoint/main
 End writing

   ** Writing the auxiliary restart file
      ----------------------------------

 Start writing
   End writing the dimensions and options               
   End writing the physical properties                  
   End writing the time step                            
  Wrote boundary condition coefficients to checkpoint: checkpoint/auxiliary
 End writing

   ** INFORMATION ON CONVERGENCE
      --------------------------
-----------------------------------------------------------------------------
   Variable    Rhs norm      N_iter  Norm. residual      Drift  Time residual
-----------------------------------------------------------------------------
c  Velocity     0.12947E+01      25   0.90507E-02   0.49766E+00 0.46610E+01 
c  Velocity[X]                                      0.20695E+00 
c  Velocity[Y]                                      0.21392E+00 
c  Velocity[Z]                                      0.76792E-01 
c  Pressure     0.54565E-03      23   0.28868E-02   0.12526E+00 0.10695E+02 
c  k            0.35459E-01      25   0.14648E-03   0.19716E-02 0.76935E+01 
c  epsilon      0.10163E+01      23   0.17973E-03   0.39222E+01 0.96520E+01 
-----------------------------------------------------------------------------
  
  

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

   field               initial min     initial max  clips to min  clips to max
-  ---------------  --------------  --------------  ------------  ------------
a  k              
a  epsilon        
a  Local Time Step


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

   field                   minimum         maximum        set mean    spatial mean
-  ---------------  --------------  --------------  --------------  --------------
v  Velocity[X]             -4.7814          6.9843          2.8105          2.8339
v  Velocity[Y]             -1.3915          7.5937          2.4045          2.4239
v  Velocity[Z]             -2.0599          2.1306      -0.0014094      0.00061584
v  ǁVelocityǁ             0.081588          8.0053          4.9385          4.9614
v  Pressure                -22.857          19.209          7.0899          7.0586
v  k                     0.0093848          1.3783         0.14199         0.13991
v  epsilon                  0.1272          228.42          3.7078          3.4189
v  TurbVisc             7.0421e-05       0.0038102      0.00087861      0.00088688
v  CourantNb               0.48739               1         0.99994         0.99993
v  FourierNb              0.014351          6.5271         0.10674         0.10626
v  total_pressure        1.013e+05      1.0134e+05      1.0133e+05      1.0133e+05
v  Local Time Step      0.00023456        0.084304       0.0010604       0.0010985

  ** Computed fields on boundary_faces
     ---------------------------------

   field         minimum         maximum        set mean    spatial mean
-  -----  --------------  --------------  --------------  --------------
v  Yplus               0          134.86          45.084          45.324


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

   criterion               minimum         maximum        set mean    spatial mean
   ---------------  --------------  --------------  --------------  --------------
   Courant/Fourier          1.0144          7.5271          1.1067          1.1062



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



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


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





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




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


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

 Destroying structures and ending computation

