command: 
 ./cs_solver --param dhc_Ax1_kwsst_htles.xml

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

                                  (R)
                      Code_Saturne

                      Version 5.1.6


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

  revision 5.1.6-fff2a0c24
  build Wed 03 Apr 2019 17:18:15 CEST
  MPI version 3.1 (MPICH 3.3b3)


  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:                Fri 17 Jul 2020 12:15:03 CEST
  System:              Linux 3.10.0-957.27.2.el7.x86_64
  Machine:             login02
  Processor:           model name	: Intel(R) Xeon(R) Gold 6130 CPU @ 2.10GHz
  Memory:              192619 MB
  User:                pbmreddy (Puneeth Bikkanahally Muni Reddy)
  Directory:           /gpfs/home/pbmreddy/CS_directory/SIMULATIONS/htles_sims/htles_base/grid_convergence/param_z/r105_z40_fs_alpha/CASE1/RESU/20200717-1214
  OpenMP threads:      1
  Processors/node:     64
  I/O read method:     standard input and output, serial access
  I/O write method:    standard input and output, serial access

  External libraries for partitioning:
    PT-SCOTCH 6.0.6

 Reading metadata from file: "mesh_input"

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

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


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



 ISCHTP =            2
@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@
@ @@   WARNING :      WHEN READING INPUT DATA
@    =========
@   CHOICE OF TIME-SCHEME
@
@     TIME-SCHEME FOR VELOCITY IS SECOND ORDER
@       (THETAV =   0.50E+00)
@     CERTAIN TERMS 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.
@
@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@


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


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

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

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


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

 --- Physics
       NVAR   =              9 (Nb variables                )
       NSCAL  =              1 (Nb scalars                  )
       NSCAUS =              0 (Nb user scalars             )
       NSCAPP =              1 (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
 omega                 1  0.1000E+01
 Alpha                 0  0.1000E+01
 TempC                 1  0.1000E+01
 Wall distance         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      )
                                                             

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


 ** STOKES
    ------

  -- Phase continue :

       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                   )


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


 ** WALL DISTANCE COMPUTATION
    -------------------------

       ICDPAR =              1 ( 1: std, reread if restart
                               (-1: std, recomputed if restrt
                               ( 2: old, reread if restart
                               (-2: old, recomputed if restrt


       NITMAY =          10000 (Nb iter for iter resolution )
       NSWRSY =              1 (Nb iter for rhs reconstr.   )
       NSWRGY =            100 (Nb iter for grad. reconstr. )
       IMLIGY =             -1 (Gradient limitation method  )
       IRCFLY =              1 (Conv. Diff. flow reconstr.  )
       ISCHCY =              1 (Convective scheme           )
       ISSTPY =              0 (Slope tet use               )
       IWARNY =              0 (Verbosity level             )
       NTCMXY =           1000 (Nb iter for steady convect. )

       BLENCY =    0.00000E+00 (2nd order conv. scheme prop.)
       EPSILY =    0.10000E-07 (Iterative solver precision  )
       EPSRSY =    0.10000E-04 (rhs reconstruction precision)
       EPSRGY =    0.10000E-04 (Gradient reconstr. precision)
       CLIMGY =    0.15000E+01 (Coeff. for grad. limitation )
       EXTRAY =    0.00000E+00 (Coeff. for grad. extrapolat.)
       COUMXY =    0.50000E+04 (Max CFL for convection      )
       EPSCVY =    0.10000E-07 (Precision for steady conv.  )
       YPLMXY =    0.20000E+03 (y+ max w. damping influence )


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


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

------------------------------------
 Variable              EPSILO IDIRCL
------------------------------------
 Velocity          0.1000E-04      1
 Pressure          0.1000E-04      1
 k                 0.1000E-07      1
 omega             0.1000E-07      1
 Alpha             0.1000E-07      0
 TempC             0.1000E-04      1
 Wall distance     0.1000E-07      1
------------------------------------

       EPSILO =                (resolution precision)
       IDIRCL = 0 ou 1         (shift diagonal if
                                ISTAT=0 and no Dirichlet)


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


 ** SCALARS
    -------

       ITBRRB =              0 (T or H reconstruction at bdy)

-------------------------------------------------------------
 Variable         Number ISCACP  ITURT     VISLS0      SIGMAS
-------------------------------------------------------------
 TempC                 1      1      0  0.3162E-05  0.1000E+01
-------------------------------------------------------------

------------------------------------
 Variable         Number      RVARFL
------------------------------------
 TempC                 1  0.8000E+00
-------------------------------------------

-------------------------------------------------------
 Variable         Number ICLVFL      SCAMIN      SCAMAX
-------------------------------------------------------
 TempC                 1     -1 -0.2731E+03  0.1000E+13
-------------------------------------------------------

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

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

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

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

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


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


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

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


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

       INPDT0 =              0 (1: 0 time step calcuation   )
       NTMABS =         800000 (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 =           1000 (Output frequency     )

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

 Velocity                       0
 Pressure                       0
 k                              0
 omega                          0
 Alpha                          0
 TempC                          0
 TempC                       -999
 TurbVisc                    -999
 Hybrid blending             -999
 HybInd1_field               -999
 HybInd2_field               -999
 HybInd3_field               -999
 HybInd4_field               -999
 HybInd5_field               -999
 Energy total                -999
 Modelised Energy            -999
 Mean Specific Di            -999
 Mean Dissipation            -999
 Ratio energy                -999
 Mean velocity               -999
 Mean Ux                     -999
 Mean Uy                     -999
 Mean Uz                     -999
 Mean UxUx                   -999
 Mean UyUy                   -999
 Mean UzUz                   -999
 F1SST_field                 -999
 taum_field                  -999
 tauij_field                 -999
 CourantNb                   -999
 FourierNb                   -999
 total_pressure              -999
 grad_temp_C                 -999
 grad_U_C                    -999
 grad_V_C                    -999
 grad_W_C                    -999
 grad_temp_FORTRA            -999
 grad_U_FORTRAN              -999
 grad_V_FORTRAN              -999
 grad_W_FORTRAN              -999
 Wall distance                  0
 Yplus                       -999
 U_mean                      -999
 V_mean                      -999
 W_mean                      -999
 UU_mean                     -999
 VV_mean                     -999
 WW_mean                     -999
 UV_mean                     -999
 VW_mean                     -999
 UW_mean                     -999
 Kmod_mean                   -999
 Omg_mean                    -999
 Pr_mean                     -999
 Nut_mean                    -999
 UP_mean                     -999
 VP_mean                     -999
 WP_mean                     -999
 T_mean                      -999
 UT_mean                     -999
 VT_mean                     -999
 WT_mean                     -999
 TT_mean                     -999
 Yplus_mean                  -999
 hybrid_blend_mea            -999
 HybInd1_mean                -999
 HybInd2_mean                -999
 HybInd3_mean                -999
 HybInd4_mean                -999
 HybInd5_mean                -999
 F1SST_mean                  -999
 TAU_M_mean                  -999
 TIJ_M_mean                  -999
 grad_U_mean_C               -999
 grad_V_mean_C               -999
 grad_W_mean_C               -999
 grad_temp_mean_C            -999
 grad_U_mean_FORT            -999
 grad_V_mean_FORT            -999
 grad_W_mean_FORT            -999
 grad_temp_mean_F            -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 =              1 (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: 
      time dependency: fixed mesh
      output: at calculation end

  -5: name: 
      directory: monitoring
      format: time_plot
      options: 
      time dependency: fixed mesh
      output: every 1 time steps

  -6: name: 
      directory: profiles
      format: plot
      options: 
      time dependency: fixed mesh
      output: at calculation end

  -7: name: histograms
      directory: histograms
      format: histogram
      options: txt
      time dependency: fixed mesh
      output: at calculation end


 Reading file:        mesh_input
 Finished reading:    mesh_input


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

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

  Selection criteria: "5 or 6"

  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:          1102240
     Number of interior faces: 3265884
     Number of boundary faces: 81672
     Number of vertices:       1143449

  Global number of boundary faces selected for joining:      55112

  Element selection successfully done.
  Determination of possible face intersections:

    bounding-box tree layout: 3D

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

  Joining operation is conforming.

  Edge intersections and vertex creation done.

  Pre-merge for 27889 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:      82668

  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:          1102240
     Number of interior faces: 3320996
     Number of boundary faces: 26560
     Number of vertices:       1143449


  Joining  1 completed (11.7 s).

  All joining operations successfully finished:

    Wall clock time:                  11.7


 ----------------------------------------------------------
 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
 Number of cells:                                      1102240

 ----------------------------------------------------------
 Number of cells + halo cells:                         1157352


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

 Number of ghost cells:                               55112

 ----------------------------------------------------------
 Number of interior faces:                             3320996

 ----------------------------------------------------------
 Number of boundary faces:                             26560

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

 Histogram of the number of interior faces per cell:

    minimum value =                  4
    maximum value =                  6

      1 : [          4 ;          5 [ =        160
      2 : [          5 ;          6 ] =    1102080

 ----------------------------------------------------------
Synchronizing cell families

 Writing file:        mesh_output
 Finished writing:    mesh_output

 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

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

 Matrix bandwidth for volume mesh :          1129796
 Matrix profile/lines for volume mesh :      35241

 Mesh coordinates:               minimum    and maximum
                       X :  0.0000000e+00  1.0000000e+00
                       Y :  0.0000000e+00  1.0000000e+00
                       Z :  0.0000000e+00  1.5000000e-01
 Mesh
     Number of cells:          1102240
     Number of interior faces: 3320996
     Number of boundary faces: 26560
     Number of vertices:       1143449

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

 Groups:
    "1"
       cells:               1102240
       boundary faces:         6640
    "2"
       boundary faces:         6640
    "3"
       boundary faces:         6640
    "4"
       boundary faces:         6640
    "5"
       interior faces:        27556
    "6"
       interior faces:        27556

 --- Information on the volumes
       Minimum control volume      =  7.3749709e-10
       Maximum control volume      =  2.2019394e-06
       Total volume for the domain =  1.5000000e-01

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

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

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

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

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

 Computing geometric quantities (5.75 s)

 Probe set: "probes":
   maximum distance between cell centers and requested coordinates: 1.656e-04
                                                             
       ALMAX  =    0.53133E+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-05  0.6000E-05
  omega             0.4610E-02  0.4610E-02
  Alpha             0.0000E+00  0.0000E+00
  TempC             0.0000E+00  0.0000E+00
  Wall distance     0.0000E+00  0.0000E+00
 ---------------------------------

                dt  0.1000E-01  0.1000E-01
 ---------------------------------



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

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



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



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


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




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


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


 -----------------------------------------
 Property           Min. value  Max. value
 -----------------------------------------
  density           0.1000E+01  0.1000E+01
  molecular_viscos  0.2245E-05  0.2245E-05
  turbulent_viscos  0.1301E-02  0.1301E-02
 -----------------------------------------

 --- Diffusivity:
 -----------------------------------------------
 Scalar           Number  Min. value  Max. value
 -----------------------------------------------
 TempC                 1  0.3162E-05  0.3162E-05
 -----------------------------------------------



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

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



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

---------------------------------------------------------------
Boundary type          Code    Nb faces           Mass flow
---------------------------------------------------------------
Inlet                     2           0         0.000000000E+00
Smooth wall               5       26560         0.000000000E+00
Rough wall                6           0         0.000000000E+00
Symmetry                  4           0         0.000000000E+00
Free outlet               3           0         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
---------------------------------------------------------------



/gpfs/softs/contrib/apps/code_saturne/5.1.6/sources/src/base/cs_field.c:1123: Fatal error.

Field "alpha"
 does not have associated BC coefficients.


Call stack:
   1: 0x7f0956bd0053 <bft_error+0xae>                 (libsaturne.so.5)
   2: 0x7f0955f39a9e <cs_f_field_bc_coeffs_ptr_by_id+0xf2> (libsaturne.so.5)
   3: 0x7f0956b67fb7 <__field_MOD_field_get_coefa_s+0x42> (libsaturne.so.5)
   4: 0x4db59c     <clptur_+0x1e7a>                 (cs_solver)
   5: 0x487b08     <condli_+0x7067>                 (cs_solver)
   6: 0x7f09560d475f <tridim_+0x8eef>                 (libsaturne.so.5)
   7: 0x7f0955e100a5 <caltri_+0x3afd>                 (libsaturne.so.5)
   8: 0x7f0955dd283a <cs_run+0x4b0>                   (libsaturne.so.5)
   9: 0x7f0955dd2ad7 <main+0x15b>                     (libsaturne.so.5)
  10: 0x7f09500603d5 <__libc_start_main+0xf5>         (libc.so.6)
  11: 0x465ba9     <>                               (cs_solver)
End of stack

