Taking a dot product between Cell-center velocity and unit normal
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=wp), | intent(in), | dimension(-2:dims%imx,-2:dims%jmx,-2:dims%kmx,-2:dims%n_var) | :: | qp | ||
type(extent), | intent(in) | :: | dims |
subroutine find_CCVn(qp, dims)
!< Taking a dot product between Cell-center velocity and unit normal
implicit none
type(extent), intent(in) :: dims
real(wp), dimension(-2:dims%imx,-2:dims%jmx,-2:dims%kmx,-2:dims%n_var), intent(in) :: qp
CCVn = CCnormalX*qp(:,:,:,2) + CCnormalY*qp(:,:,:,3) + CCnormalZ*qp(:,:,:,4) ! (nx,ny,nz).(u,v,w)
end subroutine find_CCVn