This function returns the sigma for the 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_box_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_box_expon_for_coord(me, coord, n) case ('polynomial_n5') if (me%rounded_corner) then res(:) = spongeLayer_box_roundCornerPolyn5_for_coord(me, coord, n) else res(:) = spongeLayer_box_sharpCornerPolyn5_for_coord(me, coord, n) end if case ('polynomial_n6') if (me%rounded_corner) then res(:) = spongeLayer_box_roundCornerPolyn6_for_coord(me, coord, n) else res(:) = spongeLayer_box_sharpCornerPolyn6_for_coord(me, coord, n) end if end select end function spongelayer_box_scalar_for_coord