6.3
general documentation
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.

{
k1, /* saturated permeability */
1.0, /* saturated moisture */
1.0); /* bulk density (useless) */
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");
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. */
}
cs_gwf_tracer_by_name
cs_gwf_tracer_t * cs_gwf_tracer_by_name(const char *eq_name)
Retrieve the pointer to the cs_gwf_tracer_t structure associated to the name given as parameter.
Definition: cs_gwf.c:909
cs_gwf_set_main_tracer_param
void cs_gwf_set_main_tracer_param(cs_gwf_tracer_t *tracer, const char *soil_name, double wmd, double alpha_l, double alpha_t, double distrib_coef, double reaction_rate)
Set a tracer for a specified soil when the tracer is attached to the default model.
Definition: cs_gwf_tracer.c:687
cs_gwf_soil_t
Definition: cs_gwf_soil.h:167
cs_gwf_set_iso_saturated_soil
void cs_gwf_set_iso_saturated_soil(cs_gwf_soil_t *soil, double k_s, double theta_s, double rho)
Set a soil defined by a saturated hydraulic model and attached to an isotropic permeability.
Definition: cs_gwf_soil.c:599
cs_gwf_soil_by_name
cs_gwf_soil_t * cs_gwf_soil_by_name(const char *name)
Retrieve a soil structure from its name.
Definition: cs_gwf_soil.c:521