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.
Type | Intent | Optional | 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
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