spongelayer_radial_expon_for_treeIDs Function

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

This function returns the sigma for the radial viscosity spongelayer for 2D and 3D from treeIDs

Arguments

Type IntentOptional Attributes Name
type(tem_spongeLayer_radial_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) :: nDim

Dimension
integer, intent(in) :: n

Number of arrays to return

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

return value


Calls

proc~~spongelayer_radial_expon_for_treeids~~CallsGraph proc~spongelayer_radial_expon_for_treeids spongelayer_radial_expon_for_treeIDs proc~tem_baryofid tem_BaryOfId proc~spongelayer_radial_expon_for_treeids->proc~tem_baryofid proc~tem_coordofid tem_CoordOfId proc~tem_baryofid->proc~tem_coordofid proc~tem_elemsizelevel tem_ElemSizeLevel proc~tem_baryofid->proc~tem_elemsizelevel proc~tem_levelof tem_LevelOf proc~tem_coordofid->proc~tem_levelof

Called by

proc~~spongelayer_radial_expon_for_treeids~~CalledByGraph proc~spongelayer_radial_expon_for_treeids spongelayer_radial_expon_for_treeIDs proc~spongelayer_radial_scalar_for_treeids spongelayer_radial_scalar_for_treeIDs proc~spongelayer_radial_scalar_for_treeids->proc~spongelayer_radial_expon_for_treeids interface~tem_spongelayer_radial_for tem_spongeLayer_radial_for interface~tem_spongelayer_radial_for->proc~spongelayer_radial_scalar_for_treeids proc~spongelayer_radial_vector_for_treeids spongelayer_radial_vector_for_treeIDs interface~tem_spongelayer_radial_for->proc~spongelayer_radial_vector_for_treeids proc~spongelayer_radial_vector_for_treeids->proc~spongelayer_radial_scalar_for_treeids 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 proc~tem_spatial_scalar_for_index->proc~tem_spatial_for_coord proc~tem_spatial_vector_for_index tem_spatial_vector_for_index proc~tem_spatial_vector_for_index->proc~tem_spatial_vector_for_coord

Source Code

  function spongelayer_radial_expon_for_treeIDs(me, 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
    !> 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)
    ! --------------------------------------------------------------------------
    integer :: i
    real(kind=rk) :: sigma, rad, origin(3), outerRadius, vec(3), coord(3)
    ! --------------------------------------------------------------------------
    origin(:) = me%origin
    outerRadius = me%radius + me%thickness
    vec = 0.0_rk

    do i = 1, n
      !barycentric coordinate
      coord = tem_BaryOfId( tree, treeIds(i) )
      vec(1:nDim) = coord(1:nDim) - origin(1:nDim)
      rad = sqrt( vec(1)**2 + vec(2)**2 + vec(3)**2 )

      if ( rad > me%radius .and. rad < outerRadius ) then
        sigma = me%dampFactor                                        & 
          &   * ( (rad-me%radius)/me%thickness )**me%dampExponent
      else
        sigma = 0.0_rk
      end if

      res(i) = sigma

    end do

  end function spongelayer_radial_expon_for_treeIDs