fillFinerGhostsFromMe_quadraticCompact3d Subroutine

public subroutine fillFinerGhostsFromMe_quadraticCompact3d(method, fieldProp, tLevelDesc, level, sState, snSize, tState, tnSize, tAuxField, layout, nTargets, targetList, physics, time, varSys, derVarPos)

Interpolation with momentum gradients from the stress tensor source points to the target position located at targetCoord. The source points are arranged in a square from (0,0,0)x(1,1,1) The order of the source points are according to the morton curve See tem_param for the explicit ordering

Deactivated because actually not in use Determine the coefficients of quadratic polynomial from momentum and strain rate of source elements. Evaluate the result of polynomial at targetCoord. The source points are arranged in a square from (0,0)x(1,1) The order of the source points are according to the morton curve

  1      2      3      4
 (0,0); (1,0); (0,1); (1,1)

This routine solves the equation system:

 A * x = b by x = A^(-1) * b

Deactivated, because actually not in use No comment yet!

This subroutine's interface must match the abstract interface definition intpRoutine in intp/mus_interpolate_header_module.f90 in order to be callable via do_intp function pointer.

Arguments

TypeIntentOptionalAttributesName
class(mus_interpolation_method_type), intent(in) :: method
type(mus_field_prop_type), intent(in), target:: fieldProp(:)

Array of field properties (fluid or species)

type(tem_levelDesc_type), intent(in) :: tLevelDesc

level descriptor on target level

integer, intent(in) :: level

my refinement level

real(kind=rk), intent(in) :: sState(:)

State vector of SOURCE FLUID elements

integer, intent(in) :: snSize
real(kind=rk), intent(inout) :: tState(:)

State vector of TARGET GHOST elements

integer, intent(in) :: tnSize
real(kind=rk), intent(inout) :: tAuxField(:)

AuxField variable to fill on target GHOST elements

type(mus_scheme_layout_type), intent(in) :: layout

the layout used

integer, intent(in) :: nTargets

List of target elements ( their position in depSource list )

integer, intent(in) :: targetList(nTargets)
type(mus_physics_type), intent(in) :: physics

physics type to convert lattice to physics SI unit and vice versa

type(tem_time_type), intent(in) :: time

time required to compute viscosity on target element barycenter

type(tem_varSys_type), intent(in) :: varSys

scheme variable system

type(mus_derVarPos_type), intent(in) :: derVarPos(:)

position of all derive variable in varSys for all fields


Calls

proc~~fillfinerghostsfromme_quadraticcompact3d~~CallsGraph proc~fillfinerghostsfromme_quadraticcompact3d fillFinerGhostsFromMe_quadraticCompact3d interface~tem_intp_trilinearreduced tem_intp_trilinearReduced proc~fillfinerghostsfromme_quadraticcompact3d->interface~tem_intp_trilinearreduced proc~get_moment get_moment proc~fillfinerghostsfromme_quadraticcompact3d->proc~get_moment proc~get_momentvector get_momentVector proc~get_moment->proc~get_momentvector proc~mus_imomvector mus_iMomVector proc~get_momentvector->proc~mus_imomvector

Contents


Variables

TypeVisibilityAttributesNameInitial
integer, private :: sourceLevel
integer, private :: sourceElem
integer, private :: targetLevel
integer, private :: targetElem
integer, private :: iVal
integer, private :: iElem
integer, private :: indElem
integer, private :: iSourceElem
integer, private :: nSourceElems
real(kind=rk), private :: f(layout%fStencil%QQ)
real(kind=rk), private :: fEq(layout%fStencil%QQ)
real(kind=rk), private :: sourceMom(layout%fStencil%QQ,method%nMaxSources)
real(kind=rk), private :: targetMom(layout%fStencil%QQ)
real(kind=rk), private :: targetCoord(3)
real(kind=rk), private, parameter:: targetCoord_trafo(3) =[div1_4, div1_4, div1_4]
real(kind=rk), private :: omfac
real(kind=rk), private :: tMom(layout%fStencil%QQ)
integer, private :: iCase
integer, private :: iDir
integer, private :: QQ
real(kind=rk), private :: velocity(3)
type(mus_fluid_type), private, pointer:: fluid
integer, private :: nScalars
real(kind=rk), private, allocatable:: momBuf(:,:)