update_connectivity Subroutine

private subroutine update_connectivity(scheme, iLevel)

Construct the propagation stencil for each element

the bounce back rules have to be applied here

     ! Set all the pdf links of target element to an infinitesimally
     ! small value. This will prevent density from becoming NaN if pdfs
     ! are set to zero and besides changed elements will have a very
     ! small density which will help in visualization and differentiate
     ! them from other fluid elements
     if (btest(elemProp, prp_fluidify) ) then
        scheme%state(iLevel)%val(                                         &
          &   ( ielem-1)* scheme%varsys%nscalars+idx_idir,:) =  &
          &   0.000000001_rk
     end if


type(mus_scheme_type), intent(inout) :: scheme

scheme information

integer, intent(in) :: iLevel

current level

Called by

proc~~update_connectivity~~CalledByGraph proc~update_connectivity update_connectivity proc~mus_geomincr mus_geomIncr proc~mus_geomincr->proc~update_connectivity proc~check_flow_status check_flow_status proc~check_flow_status->proc~mus_geomincr proc~mus_solve mus_solve proc~mus_solve->proc~check_flow_status program~musubi musubi program~musubi->proc~mus_solve



integer, private :: iElem
integer, private :: iField
integer, private :: idx_idir
integer, private :: QQ
integer, private :: iDir
integer, private :: invDir
integer, private :: GetFromDir
integer, private :: GetFromPos
integer, private :: neighPos
integer(kind=long_k), private :: neighProp
integer(kind=long_k), private :: elemProp
integer, private :: nElems
integer, private :: nSize