getGradU_1D Function

private pure function getGradU_1D(auxField, gradData, velPos, nAuxScalars, nDims, nSolve, elemOffset) result(gradU)

This function computes gradient of velocity from gradient and veleocity data. Gradient is computed using central difference. if an element has an boundary then neighbor refers to current element then forward difference is used

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in)			::	auxField(:)	auxField
type(mus_gradData_type),	intent(in)			::	gradData	gradient data
integer,	intent(in)			::	velPos(3)	Position of velocity field in auxField
integer,	intent(in)			::	nAuxScalars	Number of scalars in auxField array
integer,	intent(in)			::	nDims	Dimensions
integer,	intent(in)			::	nSolve	Number of element to solve in this level
integer,	intent(in)			::	elemOffset	Offset for elements when computing chunkwise

Return Value real(kind=rk)(nDims,nDims,nSolve)

output: gradient of velocity

Contents

Variables

iElem elempos leftngh rightngh leftvel rightvel

Variables

Type	Visibility	Attributes		Name		Initial
integer,	private		::	iElem
integer,	private		::	elempos
integer,	private		::	leftngh(1)
integer,	private		::	rightngh(1)
real(kind=rk),	private		::	leftvel(1,1)
real(kind=rk),	private		::	rightvel(1,1)