This function returns the sigma for the 2d box shape spongelayer
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(tem_spongeLayer_box_type) | :: | me |
Spacetime function to evaluate |
|||
real(kind=rk), | intent(in) | :: | coord(n,3) |
barycentric Ids of an elements. 1st index goes over number of elements and 2nd index goes over x,y,z coordinates |
||
integer, | intent(in) | :: | n |
Number of arrays to return |
return value
function spongelayer_box2d_scalar_for_coord(me, coord, n) &
& result(res)
! --------------------------------------------------------------------------
!> Spacetime function to evaluate
type(tem_spongeLayer_box_type) :: me
!> Number of arrays to return
integer, intent(in) :: n
!> barycentric Ids of an elements.
!! 1st index goes over number of elements and
!! 2nd index goes over x,y,z coordinates
real(kind=rk), intent( in ) :: coord(n,3)
!> return value
real(kind=rk) :: res(n)
! --------------------------------------------------------------------------
select case (trim(me%dampProfile))
case ('linear', 'exponential')
res(:) = spongeLayer_box2d_expon_for_coord(me, coord, n)
case ('polynomial_n5')
if (me%rounded_corner) then
res(:) = spongeLayer_box2d_roundCornerPolyn5_for_coord(me, coord, n)
else
res(:) = spongeLayer_box2d_sharpCornerPolyn5_for_coord(me, coord, n)
end if
case ('polynomial_n6')
if (me%rounded_corner) then
res(:) = spongeLayer_box2d_roundCornerPolyn6_for_coord(me, coord, n)
else
res(:) = spongeLayer_box2d_sharpCornerPolyn6_for_coord(me, coord, n)
end if
end select
end function spongelayer_box2d_scalar_for_coord