fillFinerGhostsFromMe_quadraticD2Q9Incomp Subroutine

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

Fill fine ghost from coarse fluid by quadratic interpolation for D2Q9 stencil. 1. Compute moments for all source elements, save in momBuf 2. For each target, interpolate moments (den, vel, tau) (10 moments for 3D and 6 moments for 2D) 3. calculate fEq and use it to calculate high order moments 4. convert moments to PDF This routine is used by acoustic quadratic interpolation for incompressible.

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.


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


proc~~fillfinerghostsfromme_quadraticd2q9incomp~~CallsGraph proc~fillfinerghostsfromme_quadraticd2q9incomp fillFinerGhostsFromMe_quadraticD2Q9Incomp proc~mus_interpolate_quad2d_leastsq mus_interpolate_quad2D_leastSq proc~fillfinerghostsfromme_quadraticd2q9incomp->proc~mus_interpolate_quad2d_leastsq proc~mus_intp_eq_d2q9_a mus_intp_eq_d2q9_a proc~fillfinerghostsfromme_quadraticd2q9incomp->proc~mus_intp_eq_d2q9_a proc~getnoneqfac getNonEqFac proc~fillfinerghostsfromme_quadraticd2q9incomp->proc~getnoneqfac



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 :: sourceMom(6,method%nMaxSources)
real(kind=rk), private :: fac
real(kind=rk), private :: tMom(layout%fStencil%QQ)
integer, private :: QQ
integer, private :: posInIntpMatLSF
type(mus_fluid_type), private, pointer:: fluid
integer, private :: nScalars
real(kind=rk), private, allocatable:: momBuf(:,:)