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

Arguments

Type IntentOptional Attributes Name
type(mus_scheme_type), intent(inout) :: scheme

scheme information
integer, intent(in) :: iLevel

current level

Calls

proc~~update_connectivity~~CallsGraph proc~update_connectivity update_connectivity cxdirinv cxdirinv proc~update_connectivity->cxdirinv nghelems nghelems proc~update_connectivity->nghelems property property proc~update_connectivity->property state_varpos state_varpos proc~update_connectivity->state_varpos val val proc~update_connectivity->val

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