compute_physFluxDivCor_nonConst Subroutine

private subroutine compute_physFluxDivCor_nonConst(nDofs, nScalars, nquadpoints, state_der, state, rot, material, poly_proj, modalCoeffs, nodalPhysFlux, pointVal, iElem, nElems)

Compute the physical flux in x direction. For other directions a properly defined variable permutation can be used. This routine covers non-constant material parameters.

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: nDofs

dimensions

integer, intent(in) :: nScalars

dimensions

integer, intent(in) :: nquadpoints

dimensions

real(kind=rk), intent(inout) :: state_der(nDofs,nScalars)

Array to store the fluxes in.

real(kind=rk), intent(in) :: state(nDofs,nScalars)

State to compute the fluxes from.

integer, intent(in) :: rot(8)

Rotationing to index the variables.

real(kind=rk), intent(in) :: material(nElems,nDofs,4)

Material parameters (mu, epsilon) for all elements

type(ply_poly_project_type), intent(inout) :: poly_proj

Data for projection method

real(kind=rk), intent(inout) :: modalCoeffs(:,:)

Working array for modal coefficients of the current element in the loop.

real(kind=rk), intent(inout) :: nodalPhysFlux(:,:)

Working array for nodal representation of the physical flux along the 3 spatial directions.

real(kind=rk), intent(inout) :: pointVal(:,:)

Working array for nodal representation of the polynomial with in each cell.

integer, intent(in) :: iElem
integer, intent(in) :: nElems

Number of elements.