get_cmax Function

private pure function get_cmax(u, x, flux_dim) result(wC)

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: u(nw_phys)
real(kind=dp), intent(in) :: x(1:ndim)
integer, intent(in) :: flux_dim

Return Value real(kind=dp)


Called by

proc~~get_cmax~8~~CalledByGraph proc~get_cmax~8 get_cmax proc~setdt~4 setdt proc~setdt~4->proc~get_cmax~8

Source Code

pure real(dp) function get_cmax(u, x, flux_dim) result(wC)
  
  real(dp), intent(in)  :: u(nw_phys)
  real(dp), intent(in)  :: x(1:ndim)
  integer, intent(in)   :: flux_dim
  real(dp)              :: mag

  mag = u(iw_b1-1+flux_dim)

  wC = dsqrt(phys_gamma*u(iw_e)/u(iw_rho)) + abs(u(iw_mom(1))*mag)

end function get_cmax