entropy_func_ECAdvection3D_t Function

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

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

Arguments

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

Return Value real(kind=prec)


Contents


Source Code

  pure function entropy_func_ECAdvection3D_t(this,s) result(e)
    !! Quadratic entropy: eta(u) = u^2 / 2
    class(ECAdvection3D_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_ECAdvection3D_t