command: 
 ./cs_solver

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

                                  (R)
                      Code_Saturne

                      Version 4.0.6


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

  revision 4.0.6
  build Sun 15 Jan 2017 07:22:30 PM MYT
  MPI version 3.0 (Open MPI 1.10.2)


  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 05 Mar 2017 06:05:17 PM MYT
  System:              Linux 4.4.0-62-generic
  Machine:             oscae01
  Processor:           model name	: Intel(R) Core(TM) i7-3770 CPU @ 3.40GHz
  Memory:              11976 MB
  User:                theuser (The User)
  Directory:           /home/theuser/ugp/test_bline_chimney/test1/RESU/20170305-1805
  I/O read method:     standard input and output, serial access
  I/O write method:    standard input and output, serial access

  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   =              7 (Nb variables                )
       NSCAL  =              1 (Nb scalars                  )
       NSCAUS =              0 (Nb user scalars             )
       NSCAPP =              1 (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.98100E+01 (Gravity z component         )

       ICORIO =          0 (Coriolis source terms       )

  -- Continuous phase:

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

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

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

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


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


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

 --- Continuous phase:

   - Commons
       ITURB  =             20 (Turbulence model            )
       IDEUCH =              1 (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)     )

   - k-epsilon           (ITURB = 20)
       ALMAX  =    0.50000E+00 (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         )


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

   - k-epsilon           (ITURB = 20)
       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
    -------------

    STEADY ALGORITHM

 --- Global parameters

       IDTVAR =             -1 (-1: steady algorithm        )
       RELXST =    0.70000E+00 (Reference relaxation coeff. )


 --- Frozen velocity field

       ICCVFG =              0 (1: frozen velocity field    )

 --- Per variable relaxation coefficient

-----------------------------
 Variable          RELAXV
-----------------------------
 Velocity          0.7000E+00
 Pressure          0.3000E+00
 Turb Kinetic Ene  0.7000E+00
 Turb Dissipation  0.7000E+00
 Temperature       0.7000E+00
----------------------------

       RELAXV =  [0.,1.]       (relaxation coefficient      )


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


 ** 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
 Turb Kinetic Ene      1      1      1      1      0 0.00E+00 0.10E+01
 Turb Dissipation      1      1      1      1      0 0.00E+00 0.10E+01
 Temperature           1      1      1      1      0 0.10E+01 0.10E+01
-------------------------------------------------------------

       ICONV  =  0 ou  1       (1 for convection active     )
       IDIFF  =  0 ou  1       (1 for total diffusion active)
       IDIFFT =  0 ou  1       (1 for turbulent diff. active)
       ISCHCV =  0 ou  1       (SOLU or CD                  )
       ISSTPC =  0 ou  1       (1: no slope test            )
       BLENCV =  [0.;1.]       (1-upwind proportion         )
       THETAV =  [0.;1.]       (0.5 Crank-Nicolson/AB       )
                               (theta for convection-       )
                               (diffusion terms uses        )
                               ((1-theta).old+theta.new     )


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


 ** STOKES
    ------
       IDILAT =              0 (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)

       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
 Turb Kinetic Ene    100      1  0.1000E-04  0.1000E-06  0.0000E+00
 Turb Dissipation    100      1  0.1000E-04  0.1000E-06  0.0000E+00
 Temperature         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
 Turb Kinetic Ene      1     -1  0.1500E+01
 Turb Dissipation      1     -1  0.1500E+01
 Temperature           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
 Turb Kinetic Ene  0.1000E-07      1
 Turb Dissipation  0.1000E-07      1
 Temperature       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    )


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


 ** SCALARS
    -------

       ITBRRB =              0 (T or H reconstruction at bdy)

-------------------------------------------------------------
 Variable         Number ISCACP ITURT      VISLS0      SIGMAS
-------------------------------------------------------------
 Temperature           1      1      0  0.1800E-04  0.1000E+01
---------------------------------------------------------------------

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

-------------------------------------------------------
 Variable         Number ICLVFL      SCAMIN      SCAMAX
-------------------------------------------------------
 Temperature           1     -1 -0.1000E+13  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.0000E+00         0  0.0000E+00
------------------------------------------------------

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


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


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

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


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

       INPDT0 =              0 (1: 0 time step calcuation   )
       NTMABS =           2000 (Final time step required    )

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


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


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

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

 --- Post-processed variables


Velocity
Pressure
Turb Kinetic Energy
Turb Dissipation
Temperature
Density
Boundary Density
Laminar Viscosity
Turb Viscosity
Total Pressure
         --           --

 --- 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 Velocity[X]                   -1
         3 Velocity[Y]                   -1
         4 Velocity[Z]                   -1
         5 Pressure                      -1
         6 Turb Kinetic Ene              -1
         7 Turb Dissipation              -1
         8 Temperature                   -1
        11 Turb Viscosity                -1
         --           --                --

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

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

         2 Velocity                       0
         5 Pressure                       0
         6 Turb Kinetic Ene               0
         7 Turb Dissipation               0
         8 Temperature                    0
         9 Density                     -999
        10 Laminar Viscosit            -999
        11 Turb Viscosity              -999
        14 Total Pressure              -999
         --           --                --

 --- Additional post-processing variables (ipstdv)
       ipstfo =              0 (Force exerted by the
                                     fluid on the boundary)
       ipstft =              0 (Thermal flux   at boundary)
       ipsttb =              0 (Temperature    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: 
      time dependency: fixed mesh
      output: at calculation end


 Reading file:        mesh_input
 Finished reading:    mesh_input

 Number of cells:                                      31250

 ----------------------------------------------------------
 Number of cells + halo cells:                         31250


 ----------------------------------------------------------
 Number of interior faces:                             90625

 ----------------------------------------------------------
 Number of boundary faces:                             6250

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

 Histogram of the number of interior faces per cell:

    minimum value =                  3
    maximum value =                  6

      1 : [          3 ;          4 [ =          8
      2 : [          4 ;          5 [ =        376
      3 : [          5 ;          6 ] =      30866

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

 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 :          1250
 Matrix profile/lines for volume mesh :      1250

 Mesh coordinates:               minimum    and maximum
                       X :  0.0000000e+00  1.0000000e+00
                       Y :  0.0000000e+00  3.1700000e-01
                       Z :  0.0000000e+00  4.2500000e-01
 Mesh
     Number of cells:          31250
     Number of interior faces: 90625
     Number of boundary faces: 6250
     Number of vertices:       34476

 Groups:
    "Back"
       boundary faces:          625
    "Bottom"
       boundary faces:         1250
    "Front"
       boundary faces:          625
    "Left"
       boundary faces:         1250
    "Right"
       boundary faces:         1250
    "Top"
       boundary faces:         1250

 --- Information on the volumes
       Minimum control volume      =  4.3112000e-06
       Maximum control volume      =  4.3112000e-06
       Total volume for the domain =  1.3472500e-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 (0.0521 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
  Turb Kinetic Ene  0.6000E-03  0.6000E-03
  Turb Dissipation  0.2645E-05  0.2645E-05
  Temperature       0.3000E+03  0.3000E+03
 ---------------------------------



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

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



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



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


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




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


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




/opt/edf-BUILD/code_saturne_build/code_saturne-4.0.6/src/base/cs_field.c:987: Fatal error.

Fortran pointer of rank 1 requested for values of field "velocity",
which have rank 2.


Call stack:
   1: 0x7fd369b7bfb8 <cs_f_field_var_ptr_by_id+0x88>  (libsaturne.so.0)
   2: 0x7fd36a0184b5 <__field_MOD_field_get_val_s+0x21> (libsaturne.so.0)
   3: 0x403866     <usphyv_+0xe2>                   (cs_solver)
   4: 0x7fd369c2508e <phyvar_+0x1cc>                  (libsaturne.so.0)
   5: 0x7fd369c490ec <tridim_+0xe75>                  (libsaturne.so.0)
   6: 0x7fd369b30354 <caltri_+0x2cd1>                 (libsaturne.so.0)
   7: 0x7fd369b08e92 <cs_run+0x412>                   (libsaturne.so.0)
   8: 0x7fd369b08961 <main+0x111>                     (libsaturne.so.0)
   9: 0x7fd368f42830 <__libc_start_main+0xf0>         (libc.so.6)
  10: 0x403119     <_start+0x29>                    (cs_solver)
End of stack

