This function calculates the sigma for the radial viscosity spongelayer for 2D and 3D from treeIDs and fills up rest with target_state. This function is currectly used to define viscosity sponge in musubi.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(tem_spongeLayer_radial_type) | :: | me |
Spacetime function to evaluate |
|||
integer, | intent(in) | :: | nComp |
Number of entrys in each array |
||
integer(kind=long_k), | intent(in) | :: | treeIds(n) |
treeIds of elements in given level |
||
type(treelmesh_type), | intent(in) | :: | tree |
global treelm mesh |
||
integer, | intent(in) | :: | nDim |
Dimension |
||
integer, | intent(in) | :: | n |
Number of arrays to return |
return value
function spongelayer_radial_vector_for_treeIDs(me, nComp, treeids, tree, &
& nDim, n) result(res)
! --------------------------------------------------------------------------
!> Spacetime function to evaluate
type(tem_spongeLayer_radial_type) :: me
!> Number of arrays to return
integer, intent(in) :: n
!> Number of entrys in each array
integer, intent(in) :: nComp
!> global treelm mesh
type( treelmesh_type ), intent(in) ::tree
!> treeIds of elements in given level
integer(kind=long_k), intent(in) :: treeIds(n)
!> Dimension
integer, intent(in) :: nDim
!> return value
real(kind=rk) :: res(n, nComp)
! --------------------------------------------------------------------------
integer :: i
! --------------------------------------------------------------------------
res(:, 1) = spongeLayer_radial_scalar_for_treeIDs(me, treeids, tree, &
& nDim, n)
if (ncomp > 1) then
do i = 1,n
res(i,2:) = me%targetState(:)
end do
end if
end function spongelayer_radial_vector_for_treeIDs