entropy_func_ECAdvection2D_t Function

  1. SELF_ECAdvection2D_t.f90
  2. SELF_ECAdvection2D_t
  3. entropy_func_ECAdvection2D_t

public pure function entropy_func_ECAdvection2D_t(this, s) result(e)

Quadratic entropy: eta(u) = u^2 / 2

Arguments

TypeIntentOptionalAttributesName
class(ECAdvection2D_t), intent(in) :: this
real(kind=prec), intent(in) :: s(1:this%nvar)

Return Value real(kind=prec)


Contents

Source Code

entropy_func_ECAdvection2D_t

Source Code

  pure function entropy_func_ECAdvection2D_t(this,s) result(e)
    !! Quadratic entropy: eta(u) = u^2 / 2
    class(ECAdvection2D_t),intent(in) :: this
    real(prec),intent(in) :: s(1:this%nvar)
    real(prec) :: e

    e = 0.5_prec*s(1)*s(1)
    if(.false.) e = e+this%u ! suppress unused-dummy-argument warning

  endfunction entropy_func_ECAdvection2D_t