adam_flume_cpu_euler_kernels

ADAM, FLUME CPU kernels of the Euler model.

The model-agnostic loop bodies (adam_flume_cpu_{face,aux}_kernels_agnostic.INC) instantiated on the Euler physics: the local arrays are sized by the Euler constants NV_K = NV_EULER, NV_AUX_K = NV_AUX (issue #41, section 4), plus the Euler signal-speed loop.

Source: src/app/flume/cpu/adam_flume_cpu_euler_kernels.F90

Dependencies

Contents

• compute_face_fluxes
• compute_q_aux
• compute_lambda_max
• face_split_fluxes

Subroutines

compute_face_fluxes

Compute the WENO face fluxes of direction d: face (i,j,k) lies between cells (i,j,k) and (i+di,j+dj,k+dk), so the face array starts at index 0 along d.

Per face: gather the stencil m = 1-S ... S into constant-bound locals, project and split it, reconstruct every field (the local v is the packed stencil of the WENO primitive), back-project.

subroutine compute_face_fluxes(d, di, dj, dk, ni, nj, nk, ngc, blocks_number, gamma, ch, is_characteristic, weno, q, q_aux, fl)

Arguments

Name	Type	Intent	Attributes	Description
d	integer(kind=I4P)	in		Direction, 1=x, 2=y, 3=z.
di	integer(kind=I4P)	in		Unit step along d.
dj	integer(kind=I4P)	in		Unit step along d.
dk	integer(kind=I4P)	in		Unit step along d.
ni	integer(kind=I4P)	in		Grid dimensions.
nj	integer(kind=I4P)	in		Grid dimensions.
nk	integer(kind=I4P)	in		Grid dimensions.
ngc	integer(kind=I4P)	in		Grid dimensions.
blocks_number	integer(kind=I4P)	in		Actual blocks number.
gamma	real(kind=R8P)	in		Specific heats ratio.
ch	real(kind=R8P)	in		GLM cleaning speed (models with GLM only).
is_characteristic	logical	in		Characteristic (or conservative) variables.
weno	type(weno_object)	in		WENO coefficients.
q	real(kind=R8P)	in		Conservative variables.
q_aux	real(kind=R8P)	in		Auxiliary variables.
fl	real(kind=R8P)	inout		Face fluxes of direction d.

Call graph

compute_q_aux

Compute the auxiliary variables on every cell, ghost cells included.

subroutine compute_q_aux(ni, nj, nk, ngc, blocks_number, gamma, R, q, q_aux)

Arguments

Name	Type	Intent	Attributes	Description
ni	integer(kind=I4P)	in		Grid dimensions.
nj	integer(kind=I4P)	in		Grid dimensions.
nk	integer(kind=I4P)	in		Grid dimensions.
ngc	integer(kind=I4P)	in		Grid dimensions.
blocks_number	integer(kind=I4P)	in		Actual blocks number.
gamma	real(kind=R8P)	in		Specific heats ratio.
R	real(kind=R8P)	in		Gas constant.
q	real(kind=R8P)	in		Conservative variables.
q_aux	real(kind=R8P)	inout		Auxiliary variables.

Call graph

compute_lambda_max

Compute max(sum_d (|u_d| + a) / dx_d) over the interior cells (null directions excluded).

subroutine compute_lambda_max(ni, nj, nk, ngc, blocks_number, gamma, R, dxyz, is_null, q, lambda_max)

Arguments

Name	Type	Intent	Attributes	Description
ni	integer(kind=I4P)	in		Grid dimensions.
nj	integer(kind=I4P)	in		Grid dimensions.
nk	integer(kind=I4P)	in		Grid dimensions.
ngc	integer(kind=I4P)	in		Grid dimensions.
blocks_number	integer(kind=I4P)	in		Actual blocks number.
gamma	real(kind=R8P)	in		Specific heats ratio.
R	real(kind=R8P)	in		Gas constant.
dxyz	real(kind=R8P)	in		Blocks space steps [3, nb].
is_null	logical	in		Null directions.
q	real(kind=R8P)	in		Conservative variables.
lambda_max	real(kind=R8P)	out		Maximum of sum_d (|u_d| + a) / dx_d.

Call graph

face_split_fluxes

Split adapter of the shared face kernel (issue #41, M2-P3): the Euler split, ch unused.

Attributes: pure

subroutine face_split_fluxes(gamma, ch, d, S, is_characteristic, qs, qas, fsplit, er)

Arguments

Name	Type	Intent	Attributes	Description
gamma	real(kind=R8P)	in		Specific heats ratio.
ch	real(kind=R8P)	in		GLM cleaning speed.
d	integer(kind=I4P)	in		Direction, 1=x, 2=y, 3=z.
S	integer(kind=I4P)	in		WENO stencil half-width, S <= S_MAX.
is_characteristic	logical	in		Characteristic (or conservative) variables.
qs	real(kind=R8P)	in		Stencil conservative variables.
qas	real(kind=R8P)	in		Stencil auxiliary variables.
fsplit	real(kind=R8P)	out		Split fields in the WENO upwind layout.
er	real(kind=R8P)	out		Right eigenvectors.

Call graph