This function returns the sigma for the spongelayer from treeids
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(tem_spongeLayer_box_type) | :: | me |
Spacetime function to evaluate |
|||
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) | :: | n |
Number of arrays to return |
return value
function spongelayer_box_scalar_for_treeIDs(me, treeids, tree, n) &
& result(res)
! --------------------------------------------------------------------------
!> Spacetime function to evaluate
type(tem_spongeLayer_box_type) :: me
!> Number of arrays to return
integer, intent(in) :: n
!> global treelm mesh
type( treelmesh_type ), intent(in) ::tree
!> treeIds of elements in given level
integer(kind=long_k), intent(in) :: treeIds(n)
!> return value
real(kind=rk) :: res(n)
! --------------------------------------------------------------------------
select case (trim(me%dampProfile))
case ('linear', 'exponential')
res(:) = spongeLayer_box_expon_for_treeIDs(me, treeIDs, tree, n)
case ('polynomial_n5')
if (me%rounded_corner) then
res(:) = spongeLayer_box_roundCornerPolyn5_for_treeIDs(me, treeIDs, &
& tree, n)
else
res(:) = spongeLayer_box_sharpCornerPolyn5_for_treeIDs(me, treeIDs, &
& tree, n)
end if
case ('polynomial_n6')
if (me%rounded_corner) then
res(:) = spongeLayer_box_roundCornerPolyn6_for_treeIDs(me, treeIDs, &
& tree, n)
else
res(:) = spongeLayer_box_sharpCornerPolyn6_for_treeIDs(me, treeIDs, &
& tree, n)
end if
end select
end function spongeLayer_box_scalar_for_treeIDs