lax_friedrichs_split_2d Subroutine

private pure subroutine lax_friedrichs_split_2d(q, dir, gam, fp, fm)

Lax-Friedrichs directional splitting: F± = ½(F ± α Q), α = |vel_dir| + c. Guarantees fp + fm == F exactly (same α applied to both halves).

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=wp),	intent(in)			::	q(neq2d)
integer,	intent(in)			::	dir
real(kind=wp),	intent(in)			::	gam
real(kind=wp),	intent(out)			::	fp(neq2d)
real(kind=wp),	intent(out)			::	fm(neq2d)