compute_rhs_cubes_modg_1d Subroutine

private subroutine compute_rhs_cubes_modg_1d(minLevel, maxLevel, currentLevel, mesh_list, tree, kerneldata_list, statedata_list, facedata_list, source, penalizationdata_list, scheme_list, poly_proj_pos, poly_proj_list, equation, material_list, general, computePenalization)

Uses

Computes the right hand side for cubical elements and 1D MODG scheme.

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: minLevel

The minimum level of the mesh.
integer, intent(in) :: maxLevel

The maximum level of the mesh.
integer, intent(in) :: currentLevel

The level of the mesh you are computing the rhs for.
type(atl_cube_elem_type), intent(inout) :: mesh_list(minLevel:maxLevel)

List of mesh parts. For each level we have one.
type(treelmesh_type), intent(in) :: tree

treelm mesh
type(atl_kerneldata_type), intent(inout) :: kerneldata_list(minLevel:maxLevel)

List of kerneldatas. For each level we have one
type(atl_statedata_type), intent(inout) :: statedata_list(minLevel:maxLevel)

List of states you want to calc the rhs for. For each level we have one.
type(atl_facedata_type), intent(inout) :: facedata_list(minLevel:maxLevel)

List of facedatas. For each level we have one
type(atl_source_type), intent(inout) :: source

Levelwise list of sources
type(atl_penalizationData_type), intent(inout) :: penalizationdata_list(minLevel:maxLevel)

Levelwise list of penalization data
type(atl_scheme_type), intent(inout) :: scheme_list(minLevel:maxLevel)

List of schemes, for each level.
integer, intent(in) :: poly_proj_pos(minLevel:maxLevel)

List of levelwise position of projection method in unique projection_list
type(ply_poly_project_type), intent(inout) :: poly_proj_list(:)

unique list for projection methods
type(atl_Equations_type), intent(in) :: equation

The equation you are operating with.
type(atl_material_type), intent(in) :: material_list(minlevel:maxlevel)

Material parameter description.
type(tem_general_type), intent(inout) :: general

General treelm settings
logical, intent(in), optional :: computePenalization

Flag to indicate whether penalization terms should be computed or not.

This is used to switch off the application of penalizing terms from the convective computation and instead compute it in an implicit update, see the atl_imexrk_module. Default is .true.!