spongelayer_radial_vector_for_treeIDs Function

private function spongelayer_radial_vector_for_treeIDs(me, nComp, treeIds, tree, nDim, n) result(res)

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.

Arguments

Type IntentOptional 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 real(kind=rk), (n,nComp)

return value


Source Code

  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