Atmospheric model (cs_user_atmospheric_model.f90)

Introduction

User subroutine for the atmospheric model.

Atmospheric module

imode corresponds to the number of calls of the usatdv function. Depending on the value of imode , different operations are performed in the following example.

imode = 0

imode = 1

Data Entry for the atmospheric soil model

To activate the model, the user has to set the and to specify the zone id on which the soil model is applied in cs_user_parameters.c, routine cs_user_model:

  at_opt->soil_model = 1; /* Switch on soil model */
 
  /* Set the number of predefined categories (+1 which is the default one)
   * among:
   *  - CS_ATMO_SOIL_5_CAT
   *  - CS_ATMO_SOIL_7_CAT
   * */
  at_opt->soil_cat= CS_ATMO_SOIL_5_CAT; /* Switch on soil model */
 
  /* Specify the boundary zone which is modeled */
  at_opt->soil_zone_id = cs_boundary_zone_by_name("Sol")->id;

Then the user may change default coefficient values for soil categories in cs_user_parameters.f90, routine usipsu

! Example to modify the soil parameters if activated
if (iatsoil.eq.1) then
  ! Example to modify some Soil constants for minerals (from Wangara test case)
  tab_sol(4)%csol = 1.7e-5
  tab_sol(4)%rugthe = 0.0012
  tab_sol(4)%rugdyn = 0.0012
endif
 
! Initializing the soil table of each vertical grid
if (iatsoil.ne.0) then
  do iiv = 1, nvert
    soilvert(iiv)%albedo  = 0.25d0
    soilvert(iiv)%emissi  = 0.965d0
    soilvert(iiv)%ttsoil  = 14.77d0
    soilvert(iiv)%totwat  = 0.0043d0
    soilvert(iiv)%pressure = 1023.d0
    soilvert(iiv)%density = 1.23d0
  enddo
endif

Initialization of atmospheric

The user has to specify the percentage of each categories for all faces of the soil zone in cs_user_initialization.c:

  if (cs_glob_atmo_option->soil_zone_id > -1) {
    const cs_zone_t *z
      = cs_boundary_zone_by_id(cs_glob_atmo_option->soil_zone_id);
 
    cs_field_t *f = cs_field_by_name("atmo_soil_percentages");
 
    for (cs_lnum_t elt_id = 0; elt_id < z->n_elts; elt_id++) {
      for (cs_lnum_t soil_id = 0; soil_id < f->dim; soil_id++)
        f->val[soil_id + f->dim * elt_id] = 0.;
 
      /* 100% of mineral */
      f->val[4 + f->dim * elt_id] = 100.;
    }
 
  }