twopointflux2d_ECDGModel2D_t Function

  1. SELF_ECDGModel2D_t.f90
  2. SELF_ECDGModel2D_t
  3. twopointflux2d_ECDGModel2D_t

public pure function twopointflux2d_ECDGModel2D_t(this, sL, sR) result(flux)

Entropy-conserving two-point flux function.

Returns the physical-space two-point EC flux for the state pair (sL, sR). flux(ivar, d) is the d-th physical component (d=1: x, d=2: y) for variable ivar.

Override this procedure in a concrete model to supply an EC flux. The stub returns zero (no flux), which is safe but trivial.

Arguments

TypeIntentOptionalAttributesName
class(ECDGModel2D_t), intent(in) :: this
real(kind=prec), intent(in) :: sL(1:this%nvar)
real(kind=prec), intent(in) :: sR(1:this%nvar)

Return Value real(kind=prec)(1:this%nvar,1:2)


Contents

Source Code

twopointflux2d_ECDGModel2D_t

Source Code

  pure function twopointflux2d_ECDGModel2D_t(this,sL,sR) result(flux)
    !! Entropy-conserving two-point flux function.
    !!
    !! Returns the physical-space two-point EC flux for the state pair (sL, sR).
    !! flux(ivar, d) is the d-th physical component (d=1: x, d=2: y) for variable ivar.
    !!
    !! Override this procedure in a concrete model to supply an EC flux.
    !! The stub returns zero (no flux), which is safe but trivial.
    implicit none
    class(ECDGModel2D_t),intent(in) :: this
    real(prec),intent(in) :: sL(1:this%nvar)
    real(prec),intent(in) :: sR(1:this%nvar)
    real(prec) :: flux(1:this%nvar,1:2)

    flux = 0.0_prec
    if(.false.) flux(1,1) = flux(1,1)+sL(1)+sR(1) ! suppress unused-argument warning

  endfunction twopointflux2d_ECDGModel2D_t