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

TypeIntentOptionalAttributesName
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.!


Calls

proc~~compute_rhs_cubes_modg_1d~~CallsGraph proc~compute_rhs_cubes_modg_1d compute_rhs_cubes_modg_1d proc~tem_stoptimer tem_stopTimer proc~compute_rhs_cubes_modg_1d->proc~tem_stoptimer proc~tem_abort tem_abort proc~compute_rhs_cubes_modg_1d->proc~tem_abort proc~atl_modg_1d_advection_numflux atl_modg_1d_advection_numFlux proc~compute_rhs_cubes_modg_1d->proc~atl_modg_1d_advection_numflux proc~atl_modg_1d_heat_numflux atl_modg_1d_heat_numFlux proc~compute_rhs_cubes_modg_1d->proc~atl_modg_1d_heat_numflux proc~atl_modg_1d_euler_numflux atl_modg_1d_euler_numFlux proc~compute_rhs_cubes_modg_1d->proc~atl_modg_1d_euler_numflux proc~modg_1d_compute_project_physflux modg_1d_compute_project_physFlux proc~compute_rhs_cubes_modg_1d->proc~modg_1d_compute_project_physflux proc~atl_modg_1d_project_testfunc atl_modg_1d_project_testFunc proc~compute_rhs_cubes_modg_1d->proc~atl_modg_1d_project_testfunc proc~tem_starttimer tem_startTimer proc~compute_rhs_cubes_modg_1d->proc~tem_starttimer proc~atl_modg_1d_finetocoarseface atl_modg_1d_fineToCoarseFace proc~compute_rhs_cubes_modg_1d->proc~atl_modg_1d_finetocoarseface proc~modg_1d_advection_onedim_numflux modg_1d_advection_oneDim_numFlux proc~atl_modg_1d_advection_numflux->proc~modg_1d_advection_onedim_numflux proc~modg_1d_heat_onedim_numflux modg_1d_heat_oneDim_numFlux proc~atl_modg_1d_heat_numflux->proc~modg_1d_heat_onedim_numflux proc~modg_1d_euler_onedim_numflux modg_1d_euler_oneDim_numFlux proc~atl_modg_1d_euler_numflux->proc~modg_1d_euler_onedim_numflux proc~modg_1d_compute_project_physflux->proc~tem_stoptimer proc~modg_1d_compute_project_physflux->proc~tem_abort proc~modg_1d_compute_project_physflux->proc~tem_starttimer proc~atl_modg_1d_euler_physflux_nonconst atl_modg_1d_euler_physFlux_nonConst proc~modg_1d_compute_project_physflux->proc~atl_modg_1d_euler_physflux_nonconst proc~atl_modg_1d_heat_physflux atl_modg_1d_heat_physFlux proc~modg_1d_compute_project_physflux->proc~atl_modg_1d_heat_physflux proc~atl_modg_1d_euler_penalization_const atl_modg_1d_euler_penalization_Const proc~modg_1d_compute_project_physflux->proc~atl_modg_1d_euler_penalization_const proc~atl_physfluxeuler_1d atl_physFluxEuler_1d proc~modg_1d_compute_project_physflux->proc~atl_physfluxeuler_1d interface~ply_poly_project_m2n ply_poly_project_m2n proc~modg_1d_compute_project_physflux->interface~ply_poly_project_m2n interface~ply_poly_project_n2m ply_poly_project_n2m proc~modg_1d_compute_project_physflux->interface~ply_poly_project_n2m proc~ply_convert2oversample ply_convert2oversample proc~modg_1d_compute_project_physflux->proc~ply_convert2oversample proc~atl_modg_1d_loclineuler_physflux atl_modg_1d_LoclinEuler_physFlux proc~modg_1d_compute_project_physflux->proc~atl_modg_1d_loclineuler_physflux proc~atl_modg_1d_advection_physflux atl_modg_1d_advection_physFlux proc~modg_1d_compute_project_physflux->proc~atl_modg_1d_advection_physflux proc~atl_modg_1d_euler_physflux_const atl_modg_1d_euler_physFlux_const proc~modg_1d_compute_project_physflux->proc~atl_modg_1d_euler_physflux_const proc~atl_modg_1d_euler_penalization_nonconst atl_modg_1d_euler_penalization_NonConst proc~modg_1d_compute_project_physflux->proc~atl_modg_1d_euler_penalization_nonconst proc~atl_initialize_state_der atl_initialize_state_der proc~modg_1d_compute_project_physflux->proc~atl_initialize_state_der proc~atl_modg_1d_project_physflux_testfunc atl_modg_1d_project_PhysFlux_testFunc proc~modg_1d_compute_project_physflux->proc~atl_modg_1d_project_physflux_testfunc proc~ply_convertfromoversample ply_convertFromOversample proc~modg_1d_compute_project_physflux->proc~ply_convertfromoversample proc~modg_1d_project_testfunc_q modg_1d_project_testFunc_Q proc~atl_modg_1d_project_testfunc->proc~modg_1d_project_testfunc_q

Called by

proc~~compute_rhs_cubes_modg_1d~~CalledByGraph proc~compute_rhs_cubes_modg_1d compute_rhs_cubes_modg_1d proc~compute_rhs_cubes compute_rhs_cubes proc~compute_rhs_cubes->proc~compute_rhs_cubes_modg_1d interface~atl_compute_rhs atl_compute_rhs interface~atl_compute_rhs->proc~compute_rhs_cubes proc~rktaylor_substep rktaylor_substep proc~rktaylor_substep->interface~atl_compute_rhs proc~rktaylor_substep->proc~rktaylor_substep proc~rk2_substep rk2_substep proc~rk2_substep->interface~atl_compute_rhs proc~rk2_substep->proc~rk2_substep proc~rk4_substep rk4_substep proc~rk4_substep->interface~atl_compute_rhs proc~rk4_substep->proc~rk4_substep proc~imexrk_substep imexrk_substep proc~imexrk_substep->interface~atl_compute_rhs proc~imexrk_substep->proc~imexrk_substep proc~mesh_timestep_euler mesh_timestep_euler proc~mesh_timestep_euler->interface~atl_compute_rhs proc~mesh_timestep_euler->proc~mesh_timestep_euler proc~global_corrector_substep global_corrector_substep proc~global_corrector_substep->interface~atl_compute_rhs proc~mesh_timestep_rktaylor mesh_timestep_rktaylor proc~mesh_timestep_rktaylor->proc~rktaylor_substep proc~mesh_timestep_predcor_cerk4 mesh_timestep_predcor_cerk4 proc~mesh_timestep_predcor_cerk4->proc~global_corrector_substep none~substep substep none~substep->proc~rk4_substep proc~mesh_timestep_imexrk mesh_timestep_imexrk proc~mesh_timestep_imexrk->proc~imexrk_substep proc~mesh_timestep_ssprk2 mesh_timestep_ssprk2 proc~mesh_timestep_ssprk2->proc~rk2_substep

Contents


Variables

TypeVisibilityAttributesNameInitial
logical, private :: usePenalization
integer, private :: iDir
integer, private :: iFace