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 | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(mus_scheme_type), | intent(inout) | :: | scheme |
scheme information |
||
| integer, | intent(in) | :: | iLevel |
current level |