/*============================================================================
 * Code_Saturne documentation page
 *============================================================================*/

/*
  This file is part of Code_Saturne, a general-purpose CFD tool.

  Copyright (C) 1998-2021 EDF S.A.

  This program 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.

  This program 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.

  You should have received a copy of the GNU General Public License along with
  this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
  Street, Fifth Floor, Boston, MA 02110-1301, USA.
*/

/*----------------------------------------------------------------------------*/

/*!
  \page cs_var_dico Variables and structures reference (C and Fortran)

  \section cs_var_dico_intro Introduction

  This page is meant to help users find their way through when implementing
  user C functions or Fortran subroutines or even developing inside the code
  kernel. It provides cross-reference tables containing the names of Fortran
  variables and their C counterparts as well as guidelines about
  how to manage the mesh entities and the fields (variables, properties, ...).
  In addition, some naming conventions are described.

  Note: variables with the same name in Fortran and in C are most of the time not listed.

  \section cs_var_dico_namingcontent Contents

  Cross-reference tables and guidelines are organized as follows:

   - \subpage mesh
   - \subpage field
   - \subpage local
   - \subpage field_kw

*/
// _____________________________________________________________________________
/*!

  \page mesh How to access and manage the main  mesh entities and mesh quantities ?

  \section cs_var_dico_mesh_vars Mesh variables

  These variables are defined in the files \ref mesh.f90 and \ref cs_mesh.h.

  - The Fortran variables are global in the code, accessible directly by their
    name.
  - Members of the global C structure \c cs_glob_mesh are accessed as: \n
    <tt>cs_glob_mesh->name</tt>, \n
    \e i.e. adding <tt> cs_glob_mesh-> </tt> in front of the name of the
    variable.
  - For a more concise syntax, defining a local \c m or \c mesh variable is recommended:
    \code{.c}
    const cs_mesh_t *m = cs_glob_mesh;
    \endcode
    so that the shorter name of the local variable may be used.
  - In functions which have a <tt> domain </tt> argument,
    using <tt> domain->mesh </tt> is recommended, and may be combined
    with the use of a local argument (as above) for conciseness.
  - Note that the number of elements given here are *local* to a given
    rank when running in parallel. Global sizes should only be used in specific
    instances, though they are available in \ref cs_mesh_t structures.

  Fortran code  | C code                                | Description
  ------------- | ------------------------------------- | ------------
  <tt> ndim     | cs_glob_mesh->dim  </tt>              | Space dimension (always 3)
  <tt> ncelet   | cs_glob_mesh->n_cells_with_ghosts </tt> | Total number of cells on the local rank \n (n_cells + n_ghost_cells)
  <tt> ncel     | cs_glob_mesh->n_cells </tt>           | Number of cells
  <tt> nfac     | cs_glob_mesh->n_i_faces </tt>         | Number of interior faces
  <tt> nfabor   | cs_glob_mesh->n_b_faces </tt>         | Number of boundary faces
  <tt> nnod     | cs_glob_mesh->n_vertices </tt>        | Number of vertices
  -             | <tt> cs_glob_mesh->n_b_cells </tt>    | Number of boundary cells
  <tt> lndfac   | cs_glob_mesh->i_face_vtx_connect_size </tt> | Size of the connectivity \n interior faces -> vertices
  <tt> lndfbr   | cs_glob_mesh->b_face_vtx_connect_size  </tt>| Size of the connectivity \n boundary faces -> vertices
  <tt> ifacel   | cs_glob_mesh->i_face_cells </tt>      | Interior faces -> cells connectivity
  <tt> ifabor   | cs_glob_mesh->b_face_cells </tt>      | Boundary faces -> cells connectivity
  <tt> ipnfac   | cs_glob_mesh->i_face_vtx_idx </tt>    | Interior faces -> vertices index
  <tt> nodfac   | cs_glob_mesh->i_face_vtx_lst </tt>    | Interior faces -> vertices connectivity
  <tt> ipnfbr   | cs_glob_mesh->b_face_vtx_idx </tt>    | Boundary faces -> vertices index
  <tt> nodfbr   | cs_glob_mesh->b_face_vtx_lst </tt>    | Boundary faces -> vertices connectivity
  -             | <tt> cs_glob_mesh->b_cells </tt>      | Boundary cell list
  <tt> xyznod   | cs_glob_mesh->vtx_coord </tt>         | Vertex coordinates

  \section cs_var_dico_mesh_q_vars Mesh quantity variables

  These variables are defined in the files \ref mesh.f90 and
  \ref cs_mesh_quantities.h.

  - The Fortran variables are global in the code, accessible directly by their
    name.
  - Members of the global C structure \c cs_glob_mesh_quantities are accessed
    as: \n
    <tt>cs_glob_mesh_quantities->name</tt>, \n
    \e i.e. adding <tt> cs_glob_mesh_quantities-> </tt> in front of the name of
    the variable.
  - For a more concise syntax, defining a local \c mq variable is recommended:
    \code{.c}
    const cs_mesh_quantities_t *mq = cs_glob_mesh_quantities;
    \endcode
    so that the shorter name of the local variable may be used.
  - In functions which have a <tt> domain </tt> argument,
    using <tt> domain->mesh_quantities </tt> is recommended, and may be combined
    with the use of a local argument (as above) for conciseness.

  Fortran code  | C code                                        | Description
  ------------- | --------------------------------------------- |-------------
  <tt> isympa   |  cs_glob_mesh_quantities->b_sym_flag </tt>    | Symmetry flag for boundary faces
  <tt> xyzcen   |  cs_glob_mesh_quantities->cell_cen </tt>      | Cell center coordinates
  <tt> surfac   |  cs_glob_mesh_quantities->i_face_normal </tt> | Surface normal of interior faces
  <tt> surfbo   |  cs_glob_mesh_quantities->b_face_normal </tt> | Surface normal of border faces
  <tt> cdgfac   |  cs_glob_mesh_quantities->i_face_cog </tt>    | Center of gravity of interior faces
  <tt> cdgfbo   |  cs_glob_mesh_quantities->b_face_cog </tt>    | Center of gravity of border faces
  <tt> volume   |  cs_glob_mesh_quantities->cell_vol </tt>      | Cell volume
  <tt> surfan   |  cs_glob_mesh_quantities->i_face_surf </tt>   | Surface of interior faces
  <tt> surfbn   |  cs_glob_mesh_quantities->b_face_surf </tt>   | Surface of boundary faces
  <tt> dist     |  cs_glob_mesh_quantities->i_dist </tt>        | Distance between the cell center and \n the center of gravity of interior faces
  <tt> distb    |  cs_glob_mesh_quantities->b_dist </tt>        | Distance between the cell center and \n the center of gravity of border faces
  <tt> pond     |  cs_glob_mesh_quantities->weight </tt>        | Interior faces weighting factor

*/
// _____________________________________________________________________________
/*!

  \page field How to access and manage variables and properties using the cs_field API?

  \ref cs_var_dico_vars "Variables" and \ref cs_var_dico_props "properties" can be accessed both in Fortran and in C using the \ref cs_field.h "cs_field" API.

  First, let us specify a few conventions used in code_saturne:

  - The *dynamic variables* designation includes velocity, pressure, and
    variables related to the turbulence model (<em>k</em>, <em>ε</em>,
    <em>R<sub>ij</sub></em>, <em>ω</em>, <em>ϕ</em>, \f$ \overline{f} \f$,
    <em>α</em>, <em>ν<sub>t</sub></em>). \n

  - The designation “scalar” refers to (usually scalar) variables which are
    solution of an advection equation, apart from the *dynamic* variables listed
    above; for instance the temperature, scalars which may be passive or not,
    user-defined or model-defined.

    The mean value of the square of the fluctuations of a “scalar” is a
    “scalar”, too.

    Scalars may be divided into two groups: **user-defined** scalars
    and **model-defined** (sometimes referred to as “specific physics”) scalars.
    In Fortran, there are \c nscaus user-defined scalars and \c nscapp
    “specific physics” (i.e. model-defined) scalars, with
    <tt>nscal = nscaus + nscapp</tt>.

  - Depending on the thermal model used; the solved thermal scalar variable may be
    the temperature, the enthalpy, or in the case of the compressible module,
    the total energy. When the temperature is not the solved thermal variable,
    it is usually available as a property.

  - When a user scalar represents the mean of the
      square of the fluctuations of another scalar (<em>i.e.</em> the average
      \f$ \overline{\varphi^\prime\varphi^\prime} \f$ for a fluctuating scalar
      <em>ϕ</em>). This can be made either {\em via} the
      interface or by declaring that scalar using
      \texttt{cs\_parameters\_add\_variable\_variance} in
      \texttt{cs\_user\_parameters.f90} (if the scalar in question is not a ``user''
      scalar, the selection is made automatically). For instance, if \texttt{j}
      and \texttt{k} are ``user'' scalars, the variable $\varphi$ corresponding
      to \texttt{k} is the variable number
      \texttt{isca(k)=isca(iscavr(j))}.
\end{list}


  \par Accessing variables and properties in Fortran:

    - The most general way of accessing a field is to use its name: \n\n

      <tt>call \ref field::field_get_val_s_by_name "field_get_val_s_by_name"("pressure", cvar_pr)
          \n pres =  cvar_pr(iel)</tt> \n\n

    - Both \ref cs_var_dico_vars "variables" and \ref cs_var_dico_props "properties" can be accessed via the
      \ref field.f90 "cs_field" API, using the global field indices and indirections arrays, as in the following examples:

      - For one-dimensional arrays:\n\n
      <tt>call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref numvar::ipr "ipr"), cvar_pr)
          \n pres =  cvar_pr(iel)</tt>, \n\n
      <tt>call \ref field::field_get_val_s "field_get_val_s"(\ref cstphy::icp "icp", cpro_cp) \n
          cp = cpro_cp(iel)</tt> \n\n
          The values of scalar variable can be accessed as follows:\n\n
       <tt>call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref isca "isca"(iscalt)), cvar_scalt) \n
          temp = cvar_scalt(iel)</tt> \n\n

      - For interleaved multidimensional arrays:\n\n
            <tt>call \ref field::field_get_val_v "field_get_val_v"(ivarfl(\ref numvar::iu "iu"), cvar_vel)
          \n ux = cvar_vel(1,iel)</tt> \n\n

       \ref numvar::ipr "ipr", \ref numvar::iu "iu" are here variable indices and the array ivarfl allows to get the corresponding field indices.\n
       \ref isca "isca"(iscalt) is also a variable index (\ref optcal::iscalt "iscalt" is the scalar index of the thermal scalar).\n
       \ref cstphy::icp "icp" is directly a field index (there is no equivalent to the array ivarfl for field of type properties).\n\n

    - the solved variable number for scalar <tt>j</tt>, where
      <tt>1 <= j <= nscal</tt>, is given by the <tt>iscal</tt> array.
      So to access the associated field, use:\n
      <tt>ivarfl(iscal(j))</tt>). \n\n

    - For model-defined scalars, the total scalar number for model scalar <tt>j</tt>, where
      <tt>1 <= j <= nscapp</tt>, is given by the <tt>iscapp</tt> array.
      So to access the associated field, use:\n
      <tt>ivarfl(iscal(iscapp(j)))</tt>). \n\n

    - The thermal scalar can be accessed using the <tt>iscalt</tt> index in the list of the
      scalars. It's variable number is thus <tt> isca(iscalt) </tt>.\n
      if there is no thermal scalar, <tt>iscalt</tt> is equal to -1. \n\n

  \par Accessing variables and properties in C:

    - The most general way of accessing a field is to use its name: \n\n

      <tt>cs_real_t *cvar_pr = cs_field_by_name "field_get_val_s_by_name"("pressure")->val; </tt>\n
      then: \n
      <tt> cvar_pr(cell_id) </tt> \n\n

    - The main \ref cs_var_dico_vars "variables" and \ref cs_var_dico_props "properties" can be accessed using the \ref CS_F_ macro:\n\n

     - For one-dimensional arrays:\n
      <tt>press = CS_F_(p)->val[cell_id]</tt>, \n
      <tt>cp = CS_F_(cp)->val[cell_id]</tt> \n
      <tt>temp = CS_F_(t)->val[cell_id]</tt> \n\n

     - For multidimensional arrays:\n
      <tt>uz = CS_F_(vel)->val[3*cell_id + 2]</tt>\n\n
      These arrays can also be cast as follows (for a 3-D array):\n
      <tt>\ref cs_real_3_t *cvar_vel = (\ref cs_real_3_t *)CS_F_(vel)->val</tt> \n\n
      The cell value can then be accessed as \n
      <tt>ux = cvar_vel[cell_id][0]</tt>\n\n

      <b>\c vel, \c p, or \c cp </b> are defined in \ref cs_field_pointer.h. \n\n

    - Indexed variables (such as user scalars) and indexed properties
      are accessed as: \n
      <tt>\ref CS_FI_(name,ii-1)->val[cell_id]</tt>. \n\n

    - The thermal scalar can be accessed using the \ref cs_thermal_model_field
      function: \n
      <tt>field_t *tf = cs_thermal_model_field();</tt>\n
      then: \n
      <tt>tf->val[cell_id]</tt>. \n
      if there is no thermal scalar, cs_thermal_model_field returns NULL. \n\n

    - In any case, all variables can be accessed using the function \ref cs_field_by_name:\n
      <tt>\ref cs_real_t *cvar_pr = \ref cs_field_by_name("pressure")->val</tt> \n\n

    - A field's \ref scalar_id key-word should be ≥ 0 if it is a scalar. \n\n

    - When a field is user-defined (rather than model-defined),
      its \ref cs_field_t::type "type" flag
      should also match the \ref CS_FIELD_USER mask: \n
      <tt> if (f->type & CS_FIELD_USER) ... </tt>

  \remark Note that indexes in C begin at 0, while indexes in Fortran begin at 1.
          Thus, Fortran and C loop counters are related in the following as:\n
          <tt>cell_id = iel-1</tt>.

  Cross-reference tables are available for the variables and properties of the
  standard solver and the specific physics features:
    - \ref cs_var_dico_vars
    - \ref cs_var_dico_props
    - \ref cs_var_dico_part
    - \ref cs_var_dico_atmo
    - \ref cs_var_dico_comb
    - \ref cs_var_dico_cfbl
    - \ref cs_var_dico_elec
    - \ref cs_var_dico_cogz
    - \ref cs_var_dico_rayt

  \section cs_var_dico_vars Variables

  The Fortran variables indexes are defined in the files \ref numvar.f90 (with
  the exception of \c ihm and \c iscal, which are respectively defined in \ref
  ppincl.f90 and \ref optcal.f90) and the C variables names are defined in
  \ref cs_field_pointer.h. \n Note that \c dt is accessible by its name
  (using the \ref cs_field_by_name family of functions), but not through
  a permanent index..

  Fortran code                             | C code                       | Description
  ------------------------------------------------ | ---------------------------- | ------------
  <tt> call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ipr), cvar_pr)  | CS_F_(p)->val       | Pressure
  call \ref field::field_get_val_v "field_get_val_v"(ivarfl(\ref iu), cvar_vel)       | CS_F_(vel)->val     | Velocity
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ivolf2), cvar_voidf) | CS_F_(void_f)->val  | Void fraction for Volume of Fluid model
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ik  ), cvar_k  )     | CS_F_(k)->val       | Turbulent kinetic energy \f$ k \f$
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref iep ), cvar_eps)     | CS_F_(eps)->val     | Turbulent dissipation \f$ \varepsilon \f$
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ir11), cvar_r11)     | CS_F_(r11)->val     | Reynolds stress component \f$ R_{xx} \f$
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ir22), cvar_r22)     | CS_F_(r22)->val     | Reynolds stress component \f$ R_{yy} \f$
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ir33), cvar_r33)     | CS_F_(r33)->val     | Reynolds stress component \f$ R_{zz} \f$
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ir12), cvar_r12)     | CS_F_(r12)->val     | Reynolds stress component \f$ R_{xy} \f$
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ir23), cvar_r23)     | CS_F_(r23)->val     | Reynolds stress component \f$ R_{yz} \f$
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ir13), cvar_r13)     | CS_F_(r13)->val     | Reynolds stress component \f$ R_{xz} \f$
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref iphi), cvar_phi)     | CS_F_(phi)->val     | \f$ \phi \f$ for \f$ \phi-f_b \f$ model
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ifb ), cvar_fb )     | CS_F_(f_bar)->val   | \f$ f_b \f$ for \f$ \phi-f_b \f$ model
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref ial ), cvar_al )     | CS_F_(alp_bl)->val  | \f$ \alpha \f$ for \f$ Bl-v^2-k \f$ \n or EBRSM model
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref iomg), cvar_omg)     | CS_F_(omg)->val     | \f$ \omega \f$ for \f$ k-\omega \f$ SST model
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref inusa), cvar_nusa)   | CS_F_(nusa)->val    | \f$ \widetilde{\nu}_T \f$ for Spalart-Allmaras
  call \ref field::field_get_val_v "field_get_val_v"(ivarfl(\ref iuma), cvar_mesh_v)  | CS_F_(mesh_u)->val  | Mesh velocity
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(isca(\ref ppincl::ihm "ihm")), cvar_hm) | CS_F_(h)->val | Enthalpy
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(isca(\ref optcal::iscalt "iscalt")), cvar_scalt) | CS_F_(t)->val | Temperature </tt>


  \section cs_var_dico_props Properties

  These properties are defined in the files \ref numvar.f90 and
  \ref cs_field_pointer.h.

  Fortran code                                                                                    |C code                             | Description
  ------------------------------------------------------------------------------------------------|---------------------------------- | ------------
  <tt> dt                                                                                         |CS_F_(dt)->val                     | Local time step
  call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::iviscl "iviscl", cpro_viscl)    |CS_F_(mu)->val                     | Molecular viscosity
  call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::ivisct "ivisct", cpro_visct)    |CS_F_(mu_t)->val                   | Turbulent dynamic viscosity
  call \ref field::field_get_val_s "field_get_val_s"(\ref cstphy::icp "icp", cpro_cp)             |CS_F_(cp)->val                     | Specific heat
  call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::icrom "icrom", cpro_crom)       |CS_F_(rho)->val                    | Density (at cells)
  call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::ibrom "ibrom", bpro_rho)        |CS_F_(rho_b)->val[face_id]         | Density (at boundary faces)
  call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::ismago "ismago", cpro_smago)    |\ref cs_real_t *cpro_smago = \ref cs_field_by_name("smagorinsky_constant^2")->val| Field id of the anisotropic turbulent viscosity
  call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::icour "icour", cpro_cour)       |\ref cs_real_t *cpro_cour = \ref cs_field_by_name("courant_number")->val | Courant number
  call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::ifour "ifour", cpro_four)       |\ref cs_real_t *cpro_four = \ref cs_field_by_name("fourier_number")->val | Fourier number
  call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::iprtot "iprtot", cpro_prtot)    |\ref cs_real_t *cpro_prtot = \ref cs_field_by_name("total_pressure")->val | Total pressure at cell centers
  call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::ivisma "ivisma", cpro_visma_s)  |\ref cs_real_t *cpro_visma_s = \ref cs_field_by_name("mesh_viscosity")->val | Mesh velocity viscosity (scalar) for the ALE module
  call \ref field::field_get_val_v "field_get_val_v"(\ref numvar::ivisma "ivisma", cpro_visma_s)  |\ref cs_real_t *cpro_visma_v = \ref cs_field_by_name("mesh_viscosity")->val | Mesh velocity viscosity (vector) for the ALE module
  call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::itsrho "itsrho"), cpro_tsrho )   |\ref cs_real_t *cpro_tsrho = \ref cs_field_by_name("dila_st")->val        | Global dilatation source terms
  call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::ibeta "ibeta"), cpro_beta  )     |\ref cs_real_t *cpro_beta = \ref cs_field_by_name("thermal_expansion")->val | Thermal expansion coefficient
  call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::ipori "ipori", cpro_ipori)       |CS_F_(poro)->val                   | Porosity
  call \ref field::field_get_val_v "field_get_val_v"(\ref numvar::iporf "iporf", cpro_iporf)       |CS_F_(t_poro)->val                 | Tensorial porosity
  call \ref field::field_get_val_v "field_get_val_v"(\ref numvar::iforbr "iforbr", bpro_forbr)     |\ref cs_real_t *bpro_forbr = \ref cs_field_by_name("boundary_forces")->val| Field id of the stresses at boundary
  call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::iyplbr "iyplbr", bpro_yplus)     |\ref cs_real_t *bpro_yplus = \ref cs_field_by_name("yplus")->val          | Field id of \f$y^+\f$ at boundary
  call \ref field::field_get_val_v "field_get_val_v"(\ref numvar::idtten "idtten", dttens)         |\ref cs_real_t *dttens = \ref cs_field_by_name("dttens")->val         | Field id for the dttens tensor
  call \ref field::field_get_val_s "field_get_val_s"(\ref numvar::itempb "itempb", t_b)            |CS_F_(t_b)->val                   | Boundary temperature </tt>

Some physical properties such as specific heat or diffusivity are often
constant (depending on the model or user parameters).
In that case, these properties are stored as a simple real numbers
rather than in a field over all mesh cells.

- This is the case for the specific heat <em>C<sub>p</sub></em>.

  - If <em>C<sub>p</sub></em> is constant, it is based on the reference
    value \ref cs_glob_fluid_properties->cp0 (<tt>cp0</tt> in Fortran).\n
    When this is the case, \ref cs_glob_fluid_properties->icp (mapped as <tt>icp</tt> in
    Fortran) should remain set to 0.\n
    When <tt>icp</tt> is initialized to 1 (by the GUI, or in
    \ref cs_user_parameters.c, it is automatically reset to the id of
    the  cell-based property field referenced in the above table.

- This is the same for the diffusivity <em>K</em> of each scalar.

  - When <em>K</em> is constant, its value is based on the field's reference
    diffusivity, accessible through the scalar field's \ref diffusivity_ref
    keyword.

    When it is variable, the matching field can be specified and accessed using
    the base scalar field's \ref diffusivity_id key (accessible using
    \ref cs_field_key_id("diffusivity_id") in C, or <tt>kivisl</tt> in Fortran).\n
    For example, for a scalar field <tt>f</tt>:\n
    <tt> int k_f_id = \ref cs_field_get_key_int(f, \ref cs_field_key_id("diffusivity_id")); \n
    cs_field_t *kf = \ref cs_field_by_id(k_f_id); </tt>


  \section cs_var_dico_part Specific physical models

  \subsection cs_var_dico_atmo Atmospheric

  Defined in \ref optcal.f90, \ref atincl.f90, \ref atvarp.f90 and
  \ref cs_field_pointer.h.

  Fortran code | C code | Description
  ------------ | ------ | ------------
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(isca(iscalt)), cvar_scalt) | <tt> CS_F_(pot_t)->val </tt>  | Potential temperature
  -            | <tt> CS_F_(totwt)->val </tt> | Total water content
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref isca(\ref atincl::intdrp "intdrp")), cvar_intdrp) | <tt> CS_F_(ntdrp)->val </tt> | Total number of droplets
  call \ref field::field_get_val_s "field_get_val_s"(ivarfl(\ref isca(\ref atchem::isca_chem "isca_chem"(iesp))), cvar_sc) | <tt> CS_FI_(chemistry,iesp-1)->val </tt>| Chemistry species (indexed)

  \subsection cs_var_dico_comb Coal combustion

  Defined in \ref ppincl.f90, \ref ppcpfu.f90 and \ref cs_field_pointer.h.

  Fortran code                                                                    | C code                           | Description
  ------------------------------------------------------------------------------- | -------------------------------- | ------------
  <tt> call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::inp "inp"(iesp)), cvar_inpcl)  | \ref CS_FI_(np,iesp-1)->val  | Particles per kg for coal class
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ixch "ixch"(iesp)), cvar_xchcl)     | \ref CS_FI_(xch,iesp-1)->val | Reactive coal mass fraction for coal class
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ixck "ixck"(iesp)), cvar_xckcl)     | \ref CS_FI_(xck,iesp-1)->val | Coke mass fraction for coal class
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ixwt "ixwt"(iesp)), cvar_xwtcl)     | \ref CS_FI_(xwt,iesp-1)->val | Water mass fraction for coal class
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ih2 "ih2"(iesp)), cvar_h2cl)        | \ref CS_FI_(h2,iesp-1)->val  | Mass enthalpy for coal class (permeatic case)
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if1m "if1m"(iesp)), cvar_f1mcl)     | \ref CS_FI_(f1m,iesp-1)->val | Mean value light volatiles for coal class
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if2m "if2m"(iesp)), cvar_f2mcl)     | \ref CS_FI_(f2m,iesp-1)->val | Mean value heavy volatiles for coal class
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if4m "if4m"), cvar_f4m)             | CS_F_(f4m)->val         | Oxydant 2 mass fraction
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if5m "if5m"), cvar_f5m))            | CS_F_(f5m)->val         | Oxydant 3 mass fraction
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if6m "if6m"), cvar_f6m))            | CS_F_(f6m)->val         | Water from coal drying mass fraction
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if7m "if7m"), cvar_f7m))            | CS_F_(f7m)->val         | Carbon from coal oxidyzed by O2 mass fraction
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if8m "if8m"), cvar_f8m))            | CS_F_(f8m)->val         | Carbon from coal gasified by CO2 mass fraction
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::if9m "if9m"), cvar_f9m))            | CS_F_(f9m)->val         | Carbon from coal gasified by H2O mass fraction
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ifvp2m "ifvp2m"), cvar_fvp2m)       | CS_F_(fvp2m)->val       | f1f2 variance
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppcpfu::iyco2 "iyco2"), cvar_yco2)          | CS_F_(yco2)->val        | CO2 fraction
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppcpfu::iyhcn "iyhcn"), cvar_yhnc)          | CS_F_(yhcn)->val        | HCN fraction
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppcpfu::iyno "iyno"), cvar, yno)            | CS_F_(yno)->val         | NO fraction
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppcpfu::iynh3 "iynh3"), cvar_ynh3)          | CS_F_(ynh3)->val        | NH3 enthalpy
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppcpfu::ihox "ihox"), cvar_hox)             | CS_F_(hox)->val         | Ox enthalpy </tt>


  \subsection cs_var_dico_cfbl Compressible

  Defined in \ref ppincl.f90 and \ref cs_field_pointer.h.

  Fortran code                                                                   | C code                        | Description
  ------------------------------------------------------------------------------ | ----------------------------- | ------------
  <tt> call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ienerg "ienerg"), cvar_energ) | CS_F_(e_tot)->val   | Total energy
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::itempk "itempk"), cvar_tempk)      | CS_F_(t_kelvin)->val | Temperature, in Kelvin </tt>


  \subsection cs_var_dico_elec Electric arcs

  Defined in \ref ppincl.f90 and \ref cs_field_pointer.h.

  Fortran code                                                                  | C code                             | Description
  ----------------------------------------------------------------------------- | ---------------------------------- | ------------
  <tt> call \ref field::field_get_val_s_by_name "field_get_val_s_by_name"("elec_pot_r", cvar_potr)   | CS_F_(potr)->val          | Electric potential, real part
  call \ref field::field_get_val_s_by_name "field_get_val_s_by_name"("elec_pot_i", cvar_poti)        | CS_F_(poti)->val          | Electric potential, imaginary part
  call \ref field::field_get_val_v_by_name "field_get_val_v_by_name"("vec_potential", cvar_potva)    | CS_F_(potva)->val         | Vector potential
  call \ref field::field_get_val_s_by_name "field_get_val_s_by_name"("esl_fraction_01", cvar_ycoel_01) | \ref CS_FI_(ycoel,iesp-1)->val | Constituent mass fraction </tt>


  \subsection cs_var_dico_cogz Gas combustion

  Defined in \ref ppincl.f90 and \ref cs_field_pointer.h.

  Fortran code                                                              | C code                     | Description
  ------------------------------------------------------------------------- | -------------------------- | ------------
  <tt> call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ifm "ifm"), cvar_fm)     | CS_F_(fm)->val    | Mixture fraction
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ifp2m "ifp2m"), cvar_fp2m)    | CS_F_(fp2m)->val  | Mixture fraction variance
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::ifsm "ifsm"), cvar_fsm)       | CS_F_(fsm)->val   | Soot mass fraction
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::inpm "inpm"), cvar_npm)       | CS_F_(npm)->val   | Soot precursor number
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::iygfm "iygfm"), cvar_ygfm)    | CS_F_(ygfm)->val  | Fresh gas fraction
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::iyfm "iyfm"), cvar_yfm)       | CS_F_(yfm)->val   | Mass fraction
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::iyfp2m "iyfp2m"), cvar_yfp2m) | CS_F_(yfp2m)->val | Mass fraction variance
  call \ref field::field_get_val_s "field_get_val_s"(\ref isca(\ref ppincl::icoyfp "icoyfp"), cvar_coyfp) | CS_F_(coyfp)->val | Mass fraction covariance </tt>


  \subsection cs_var_dico_rayt Radiative transfer

  Defined in \ref cs_field_pointer.h.

  C code                                  | Description
  --------------------------------------- | ------------
  <tt> CS_F_(rad_lumin)->val </tt>        | Radiative luminance
  <tt> CS_F_(rad_q)->val </tt>            | Radiative flux
  <tt> CS_FI_(rad_ets,iesp-1)->val </tt>  | Radiative flux explicit source term
  <tt> CS_FI_(rad_its,iesp-1)->val </tt>  | Radiative flux implicit source term
  <tt> CS_FI_(rad_abs,iesp-1)->val </tt>  | Radiative absorption
  <tt> CS_FI_(rad_emi,iesp-1)->val </tt>  | Radiative emission
  <tt> CS_FI_(rad_cak,iesp-1)->val </tt>  | Radiative absorption coefficient
  <tt> CS_F_(qinci)->val </tt>            | Radiative incident radiative flux density
  <tt> CS_F_(xlam)->val </tt>             | Wall thermal conductivity
  <tt> CS_F_(epa)->val </tt>              | Wall thickness
  <tt> CS_F_(emissivity)->val </tt>       | Wall emissivity
  <tt> CS_F_(fnet)->val </tt>             | Boundary radiative flux
  <tt> CS_F_(fconv)->val </tt>            | Boundary radiative convective flux
  <tt> CS_F_(hconv)->val </tt>            | Radiative exchange coefficient


  \subsection cs_var_dico_multiphase Eulerian-Eulerian multiphase flows

  Defined in \ref cs_field_pointer.h.

  C code                                  | Description
  --------------------------------------- | ------------
  <tt> CS_FI_(yf_ncond,inc)->val </tt>    | Non-condensable gas mass fractions
  <tt> CS_FI_(qp,ip)->val </tt>           | Particles turbulent kinetic energy Q2
  <tt> CS_FI_(qfp,ip)->val </tt>          | Covariance of the turbulent Q12
  <tt> CS_FI_(qfpxx,ip)->val </tt>        | XX component of qfp
  <tt> CS_FI_(qfpxy,ip)->val </tt>        | XY component of qfp
  <tt> CS_FI_(qfpxz,ip)->val </tt>        | XZ component of qfp
  <tt> CS_FI_(qfpyx,ip)->val </tt>        | YX component of qfp
  <tt> CS_FI_(qfpyy,ip)->val </tt>        | YY component of qfp
  <tt> CS_FI_(qfpyz,ip)->val </tt>        | YZ component of qfp
  <tt> CS_FI_(qfpzx,ip)->val </tt>        | ZX component of qfp
  <tt> CS_FI_(qfpzy,ip)->val </tt>        | ZY component of qfp
  <tt> CS_FI_(qfpzz,ip)->val </tt>        | ZZ component of qfp
  <tt> CS_FI_(gamma,ip)->val </tt>        | Interfacial mass transfer
  <tt> CS_FI_(ia,ip)->val </tt>           | Interfacial area
  <tt> CS_FI_(x2,ip)->val </tt>           | Droplets x2
  <tt> CS_FI_(d32,ip)->val </tt>          | Droplets Sauter mean diameter
  <tt> CS_FI_(drag,ipcpl)->val </tt>      | Drag between phases
  <tt> CS_FI_(ad_mass,ip)->val </tt>      | Added mass
  <tt> CS_FI_(th_diff,ip)->val </tt>      | Thermal diffusivity
  <tt> CS_FI_(th_diff_t,ip)->val </tt>    | Turbulent thermal diffusivity
  <tt> CS_FI_(drho_dp,ip)->val </tt>      | dRho over dP
  <tt> CS_FI_(drho_dh,ip)->val </tt>      | dRho over dH
  <tt> CS_FI_(tau12_t,ip)->val </tt>      | Turbulent tau12 for particles
  <tt> CS_FI_(lift,ip)->val </tt>         | Particles lift
  <tt> CS_FI_(disp_t,ip)->val </tt>       | Particles turbulent dispersion
  <tt> CS_FI_(drift_vel,ip)->val </tt>    | Particles drift velocity

*/
// _____________________________________________________________________________
/*!

  \page local How to name common local variables ?

  The following table provides a non-exhaustive list of local variables which
  are used in the code in a recurring manner.

  Fortran code       | C code     | Description
  ------------------ | ---------- | ------------
  <tt> iel           | cell_id    </tt> | Cell index
  <tt> ifac          | face_id    </tt> | Face index
  <tt> ig            | g_id       </tt> | Interior face number of associated groups (OpenMP)
  <tt> it            | t_id       </tt> | Interior face number of threads (OpenMP)
  <tt> idimtr        | tr_dim     </tt> | Indicator for tensor perodicity of rotation
  <tt> flumas        | i_massflux </tt> | Mass flux at interior faces
  <tt> flumab        | b_massflux </tt> | Mass flux at boundary faces
  <tt> viscf         | i_visc     </tt> | \f$ \mu_\fij \dfrac{S_\fij}{\ipf \jpf} \f$ \n  at interior faces for the r.h.s.
  <tt> viscb         | b_visc     </tt> | \f$ \mu_\fib \dfrac{S_\fib}{\ipf \centf} \f$ \n  at border faces for the r.h.s.
  <tt> smbrp         | rhs        </tt> | Right hand side \f$ \vect{Rhs} \f$

  \section cs_var_dico_conv Local naming convention for fields (Fortran and C)

  Rules have been stablished for local names denoting fields, depending on their nature. The convention, applying both in Fortran and in C, is as follows:

  - The first letter of the name indicates the location at which the field values are defined:
    - \b c for values at the cell centers.
    - \b i for values at the interior faces.
    - \b b for values at the boundary faces.
  - The next three letters indicate if the field is a variable (at the current time step or the previous time step) or a property:
    - \b var for variables at the current time step.
    - \b vara for variables at the previous time step.
    - \b pro for properties.
  - An underscore \b _ follows.
  - Finally, the <b> short name </b> of the variable/property is specified. This short name is built from the variable/property Fortran index, removing the \c i at the beginning of the word.

  The following examples ilustrate this convention:

  \c cvar_pr: Values of the variable pressure field defined at the cell centers, at the current time step. \n
  \c cvara_pr: Values of the variable pressure field defined at the cell centers, at the previous time step. \n
  \c cpro_cp: Values of the property specific heat defined field at the cell centers. \n

*/
// _____________________________________________________________________________
/*!

  \page field_kw List of main field keys

  A non-exhaustive \ref field_keywords "list of field keys" which are used in the code is provided.

*/
