spongelayer_radial_vector_for_coord Function

private function spongelayer_radial_vector_for_coord(me, nComp, coord, nDim, n) result(res)

This function calculates the sigma for the radial viscosity spongelayer for 2D and 3D 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

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) :: nDim

Dimension

integer, intent(in) :: n

Number of arrays to return

Return Value real(kind=rk), (n,nComp)

return value


Calls

proc~~spongelayer_radial_vector_for_coord~~CallsGraph proc~spongelayer_radial_vector_for_coord spongelayer_radial_vector_for_coord proc~spongelayer_radial_scalar_for_coord spongelayer_radial_scalar_for_coord proc~spongelayer_radial_vector_for_coord->proc~spongelayer_radial_scalar_for_coord proc~spongelayer_radial_expon_for_coord spongelayer_radial_expon_for_coord proc~spongelayer_radial_scalar_for_coord->proc~spongelayer_radial_expon_for_coord proc~spongelayer_radial_polyn5_for_coord spongeLayer_radial_polyn5_for_coord proc~spongelayer_radial_scalar_for_coord->proc~spongelayer_radial_polyn5_for_coord proc~spongelayer_radial_polyn6_for_coord spongeLayer_radial_polyn6_for_coord proc~spongelayer_radial_scalar_for_coord->proc~spongelayer_radial_polyn6_for_coord

Called by

proc~~spongelayer_radial_vector_for_coord~~CalledByGraph proc~spongelayer_radial_vector_for_coord spongelayer_radial_vector_for_coord interface~tem_spongelayer_radial_for tem_spongeLayer_radial_for interface~tem_spongelayer_radial_for->proc~spongelayer_radial_vector_for_coord proc~tem_spatial_for_coord tem_spatial_for_coord proc~tem_spatial_for_coord->interface~tem_spongelayer_radial_for proc~tem_spatial_for_treeids tem_spatial_for_treeIDs proc~tem_spatial_for_treeids->interface~tem_spongelayer_radial_for proc~tem_spatial_vector_for_coord tem_spatial_vector_for_coord proc~tem_spatial_vector_for_coord->interface~tem_spongelayer_radial_for proc~tem_spatial_vector_for_treeids tem_spatial_vector_for_treeIDs proc~tem_spatial_vector_for_treeids->interface~tem_spongelayer_radial_for interface~tem_spatial_for tem_spatial_for interface~tem_spatial_for->proc~tem_spatial_for_coord interface~tem_spatial_for->proc~tem_spatial_for_treeids interface~tem_spatial_for->proc~tem_spatial_vector_for_coord interface~tem_spatial_for->proc~tem_spatial_vector_for_treeids proc~tem_spatial_scalar_for_index tem_spatial_scalar_for_index interface~tem_spatial_for->proc~tem_spatial_scalar_for_index proc~tem_spatial_vector_for_index tem_spatial_vector_for_index interface~tem_spatial_for->proc~tem_spatial_vector_for_index proc~tem_spatial_scalar_for_index->proc~tem_spatial_for_coord proc~tem_spatial_vector_for_index->proc~tem_spatial_vector_for_coord proc~tem_spacetime_for_coord tem_spacetime_for_coord proc~tem_spacetime_for_coord->interface~tem_spatial_for proc~tem_spacetime_for_treeids tem_spacetime_for_treeIDs proc~tem_spacetime_for_treeids->interface~tem_spatial_for proc~tem_spacetime_scalar_for_index tem_spacetime_scalar_for_index proc~tem_spacetime_scalar_for_index->interface~tem_spatial_for proc~tem_spacetime_vector_for_coord tem_spacetime_vector_for_coord proc~tem_spacetime_vector_for_coord->interface~tem_spatial_for proc~tem_spacetime_vector_for_index tem_spacetime_vector_for_index proc~tem_spacetime_vector_for_index->interface~tem_spatial_for proc~tem_spacetime_vector_for_treeids tem_spacetime_vector_for_treeIDs proc~tem_spacetime_vector_for_treeids->interface~tem_spatial_for proc~tem_spatial_scalar_storeval tem_spatial_scalar_storeVal proc~tem_spatial_scalar_storeval->interface~tem_spatial_for proc~tem_spatial_vector_storeval tem_spatial_vector_storeVal proc~tem_spatial_vector_storeval->interface~tem_spatial_for

Source Code

  function spongelayer_radial_vector_for_coord(me, nComp, coord, 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
    !> 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)
    !> Dimension
    integer, intent(in) :: nDim
    !> return value
    real(kind=rk) :: res(n, nComp)
    ! --------------------------------------------------------------------------
    integer :: i
    ! --------------------------------------------------------------------------
    res(:, 1) = spongeLayer_radial_scalar_for_coord(me, coord, 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_coord