This function calculates the sigma for the box shape spongelayer fom treeid for polynomial n5 profile
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_box2d_sharpCornerPolyn5_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) ! -------------------------------------------------------------------------- integer :: i real(kind=rk) :: res_i(1), coord(3,1) ! -------------------------------------------------------------------------- do i = 1,n !barycentric coordinate coord(:,1) = tem_BaryOfId( tree, treeIds(i) ) res_i = spongeLayer_box2d_sharpCornerPolyn5_for_coord(me, coord, 1) res(i) = res_i(1) end do end function spongelayer_box2d_sharpCornerPolyn5_for_treeIDs