/* SRC/cs_user_les_inflow.cpp */

#include "base/cs_defs.h"
#include <assert.h>
#include <math.h>
#include <string.h>

#if defined(HAVE_MPI)
#include <mpi.h>
#endif

#include <ple_coupling.h>

#include "cs_headers.h"


BEGIN_C_DECLS

void
cs_user_les_inflow_define(void)
{
  cs_les_inflow_set_restart(false,   /* allow_read */
                    false);  /* allow_write */

  /* First synthetic turbulence inlet: the Batten Method is used
     for boundary faces of zone "INLET" */

  {
    /* Number of turbulence structures */
    int n_entities = 50;

    /* Velocity, turbulent kinetic energy and dissipation rate */
    cs_real_t vel_r[3] = {2.932153, 0, 0};
    cs_real_t k_r = 0.0322;
    cs_real_t eps_r = 6.6;

    cs_les_inflow_add_inlet(CS_INFLOW_SEM,
                            false,
                            cs_boundary_zone_by_name("INLET"),
                            n_entities,
                            1,  /* verbosity */
                            vel_r,
                            k_r,
                            eps_r);
  }


}

END_C_DECLS