Settings soils and tracers for the groundwater flow module with CDO schemes (cs_user_gwf.c)

The add of a new soil or a new tracer takes place in cs_user_model (see Settings related to the groundwater flow module with CDO/HHO schemes for more details).

Soils

Example for a saturated soil defined by an isotropic permeability.

  {
    cs_gwf_soil_t  *s1 = cs_gwf_soil_by_name("soil1");
    cs_gwf_set_iso_saturated_soil(s1,
                                  k1,     /* saturated permeability */
                                  1.0,    /* saturated moisture */
                                  1.0);   /* bulk density (useless) */

    cs_gwf_soil_t  *s2 = cs_gwf_soil_by_name("soil2");
    cs_gwf_set_iso_saturated_soil(s2,
                                  k2,     /* saturated permeability */
                                  1.0,    /* saturated moisture */
                                  1.0);   /* bulk density (useless) */
  }

cs_gwf_set_aniso_saturated_soil sets a saturated soil defined by an anisotropic permeability. Soils which behave according to a Van Genuchten model can be specified using cs_gwf_set_aniso_genuchten_soil or cs_gwf_set_iso_genuchten_soil. More advanced definition using a user-defined model is also possible using cs_gwf_set_user_soil.

Tracers

Here is an example for a standard tracer.

  {
    cs_gwf_tracer_t *tr = cs_gwf_tracer_by_name("Tracer_01");
    cs_gwf_set_main_tracer_param(tr,
                                 NULL,     /* soil name or NULL for all */
                                 0.,       /* water molecular diffusivity */
                                 0., 0.,   /* alpha (longi. and transvesal) */
                                 0.,       /* distribution coef. */
                                 0.);      /* 1st order decay coef. */
  }