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
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 |
output: gradient of velocity
Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|
integer, | private | :: | iElem | ||||
integer, | private | :: | elempos | ||||
integer, | private | :: | leftngh(3) | ||||
integer, | private | :: | rightngh(3) | ||||
real(kind=rk), | private | :: | leftvel(3,3) | ||||
real(kind=rk), | private | :: | rightvel(3,3) |