command: 
 /home/zhang/cfdexe/codeStaure/libexec/code_saturne/cs_solver --param setup.xml --log 0 --logp 1 --mpi

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

                                  (R)
                      Code_Saturne

                      Version 5.0-beta


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

  revision 5.0-beta-r9986
  build Mon 08 May 2017 12:31:43 PM CST
  MPI version 2.2 (MPICH2 1.4.1p1)


  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.

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


 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 )


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


 ** VOF METHOD                              
    ----------------------------------------

       IVOFMT =             -1 ( -1: inactive                )
                               (  0: active                  )


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

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


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


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

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

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


 ** STOKES
    ------

  -- Phase continue :

       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                   )


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


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

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

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

 Velocity                       0
 Pressure                       0
 k                              0
 epsilon                        0
 TurbVisc                    -999
 CourantNb                   -999
 FourierNb                   -999
 total_pressure              -999
 Local Time Step             -999
 Yplus                       -999
         --           --                --

 --- Additional post-processing variables (ipstdv)
       ipstfo =              1 (Force exerted by the
                                     fluid on the boundary)
       ipstyp =              1 (y+ at boundary)
       ipsttp =              0 (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)


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


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

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

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

 Partitioning finished (0.144 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 =              10880
    maximum value =              10880

      1 : [      10880 ;      10880 ] =          5

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

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

    minimum value =              11301
    maximum value =              12620

      1 : [      11301 ;      11564 [ =          1
      2 : [      11564 ;      11828 [ =          0
      3 : [      11828 ;      12092 [ =          1
      4 : [      12092 ;      12356 [ =          0
      5 : [      12356 ;      12620 ] =          3

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

 Histogram of the number of ghost cells per rank:

    minimum value =                421
    maximum value =               1740

      1 : [        421 ;        684 [ =          1
      2 : [        684 ;        948 [ =          0
      3 : [        948 ;       1212 [ =          1
      4 : [       1212 ;       1476 [ =          0
      5 : [       1476 ;       1740 ] =          3

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

 Histogram of the number of interior faces per rank:

    minimum value =              21812
    maximum value =              22886

      1 : [      21812 ;      22026 [ =          1
      2 : [      22026 ;      22241 [ =          0
      3 : [      22241 ;      22456 [ =          1
      4 : [      22456 ;      22671 [ =          1
      5 : [      22671 ;      22886 ] =          2

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

 Histogram of the number of boundary faces per rank:

    minimum value =              21841
    maximum value =              22119

      1 : [      21841 ;      21896 [ =          2
      2 : [      21896 ;      21952 [ =          1
      3 : [      21952 ;      22007 [ =          0
      4 : [      22007 ;      22063 [ =          1
      5 : [      22063 ;      22119 ] =          1

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

 Histogram of the number of interior faces per cell:

    minimum value =                  2
    maximum value =                  4

      1 : [          2 ;          3 [ =          4
      2 : [          3 ;          4 ] =      54396

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

 Histogram of the number of neighboring domains per rank:

    minimum value =                  3
    maximum value =                  4

      1 : [          3 ;          4 ] =          5

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

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

 Renumbering mesh:

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

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

   renumbering for boundary faces:
     numbering:                           none

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

 Histogram of volume mesh matrix bandwidth per rank:

    minimum value =              10292
    maximum value =              11594

      1 : [      10292 ;      10552 [ =          1
      2 : [      10552 ;      10812 [ =          0
      3 : [      10812 ;      11073 [ =          3
      4 : [      11073 ;      11333 [ =          0
      5 : [      11333 ;      11594 ] =          1

 Histogram of volume mesh matrix profile/lines per rank:

    minimum value =                442
    maximum value =               1434

      1 : [        442 ;        640 [ =          2
      2 : [        640 ;        838 [ =          1
      3 : [        838 ;       1037 [ =          0
      4 : [       1037 ;       1235 [ =          1
      5 : [       1235 ;       1434 ] =          1

 Mesh coordinates:               minimum    and maximum
                       X : -6.9000000e+02  1.6100000e+03
                       Y : -3.8985001e+02  3.8985001e+02
                       Z : -5.9799999e+01  2.3871429e-14
 Mesh
     Number of cells:          54400
     Number of interior faces: 108312
     Number of boundary faces: 109776
     Number of vertices:       109776

 Groups:
    "DOWN"
       boundary faces:          220
    "FAR"
       boundary faces:          696
    "FOIL"
       boundary faces:          280
    "INLET"
       boundary faces:          128
    "OUT"
       boundary faces:          128
    "SYM"
       boundary faces:        54400
    "SYM1"
       boundary faces:        54400
    "UP"
       boundary faces:          220

 --- Information on the volumes
       Minimum control volume      =  1.4091005e-03
       Maximum control volume      =  3.6252512e+04
       Total volume for the domain =  1.0699132e+08

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

  Criterion 2: Offset:
    Number of bad cells detected: 148 -->   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 %

 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 (0.0122 s)
                                                             

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


 ** 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.1162E+01  0.1162E+01
  epsilon           0.1127E+00  0.1127E+00
 ---------------------------------

                dt  0.1000E+00  0.1000E+00
 ---------------------------------



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

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



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



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


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




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


 INSTANT    0.100000000E+00   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.1270E+01  0.1270E+01
 -----------------------------------------



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

-------------------------------------------------------------------------
Boundary type          Code    Nb faces
-------------------------------------------------------------------------
Inlet                     2         128
Smooth wall               5         280
Rough wall                6           0
Symmetry                  4      109240
Free outlet               3         128
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.16100E+04  -0.12923E+01  -0.29900E+02



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

---------------------------------------------------------------
Boundary type          Code    Nb faces           Mass flow
---------------------------------------------------------------
Inlet                     2         128         0.000000000E+00
Smooth wall               5         280         0.000000000E+00
Rough wall                6           0         0.000000000E+00
Symmetry                  4      109240         0.000000000E+00
Free outlet               3         128         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
   ---------------------------------------
------------------------------------------------------------
                                         Minimum     Maximum
------------------------------------------------------------
   Rel velocity at the wall uiptn :  0.00000E+00 0.00000E+00
   Friction velocity        uet   :  0.00000E+00 0.00000E+00
   Friction velocity        uk    :  0.55797E+00 0.58215E+00
   Dimensionless distance   yplus :  0.12748E+02 0.21351E+02
   ------------------------------------------------------
   Nb of reversal of the velocity at the wall   :          0
   Nb of faces within the viscous sub-layer     :          0
   Total number of wall faces                   :        280
------------------------------------------------------------



   ** INFORMATION ON CONVERGENCE
      --------------------------
----------------------------------------------------------------------------
   Variable     Rhs norm      N_iter  Norm. residual   Drift   Time residual
----------------------------------------------------------------------------
c  Velocity             NaN       1   0.00000E+00           NaN         NaN 
c  Velocity[X]                                              NaN 
c  Velocity[Y]                                              NaN 
c  Velocity[Z]                                              NaN 
c  Pressure             NaN       0   0.00000E+00           NaN 0.00000E+00 
c  k                    NaN       1   0.00000E+00           NaN         NaN 
c  epsilon              NaN       1   0.00000E+00           NaN         NaN 
----------------------------------------------------------------------------
  
  
