mms_diffusion_m Module

Method of Manufactured Solutions for diffusion model

Uses

- circular_equilibrium_m
- variable_m
- mms_diffusion_dommaschk_m
- error_handling_grillix_m
- mms_diffusion_slab_m
- screen_io_m
- parallelisation_setup_m
- status_codes_grillix_m
- dommaschk_equilibrium_m
- slab_equilibrium_m
- precision_grillix_m
- mpi
- equilibrium_m
- parallel_map_m
- mesh_cart_m
- mms_diffusion_circular_m
- commons_misc_m
- comm_handler_m

Contents

Functions

mms_sol_diffusion mms_source_diffusion

Subroutines

mms_diagnostics_diffusion

Functions

public function mms_sol_diffusion(equi, x, y, z, t)

MMS solution

Arguments

Type	Intent	Optional	Attributes		Name
class(equilibrium_t)				::	equi	Equilibrium
real(kind=GP),	intent(in)			::	x	x-coordinate
real(kind=GP),	intent(in)			::	y	y-coordinate
real(kind=GP),	intent(in)			::	z	Axial coordinate z (represented by toroidal coordinate)
real(kind=GP),	intent(in)			::	t	Time

Return Value real(kind=gp)

public function mms_source_diffusion(equi, x, y, z, t, charfun, penval)

Source corresponding to MMS solution

Arguments

Type	Intent	Optional	Attributes		Name
class(equilibrium_t),	intent(inout)			::	equi	Equilibrium
real(kind=GP),	intent(in)			::	x	x-coordinate
real(kind=GP),	intent(in)			::	y	y-coordinate
real(kind=GP),	intent(in)			::	z	Axial coordinate z (represented by toroidal coordinate)
real(kind=GP),	intent(in)			::	t	Time
real(kind=GP),	intent(in)			::	charfun	Value of characteristic function (penalisation)
real(kind=GP),	intent(in)			::	penval	Value that solution shall be penalised to

Return Value real(kind=gp)

Subroutines

public subroutine mms_diagnostics_diffusion(comm_handler, equi, mesh, map, tau, num_sol)

Prints information on numerical errors of numerical solution

Arguments

Type	Intent	Optional	Attributes		Name
type(comm_handler_t),	intent(in)			::	comm_handler
class(equilibrium_t),	intent(inout)			::	equi	Equilibrium
type(mesh_cart_t),	intent(in)			::	mesh	Mesh within poloidal plane
type(parallel_map_t),	intent(in)			::	map	Mesh within poloidal plane
real(kind=GP),	intent(in)			::	tau	Time
type(variable_t),	intent(in)			::	num_sol	Numerical solution