ic_tgv_uy_for Function

  pure function ic_tgv_uy_for( me, coord, n ) result( uy )
    ! ---------------------------------------------------------------------------
    !> global gauss pulse data
    type(ic_tgv_type), intent(in) :: me
    !> number of return values
    integer, intent(in) :: n
    !> coordinate of an element
    real(kind=rk), intent(in) :: coord(n, 3)
    !> return value which is sent to state variable
    real(kind=rk) :: uy(n)
    ! ---------------------------------------------------------------------------
    integer :: iElem
    ! ---------------------------------------------------------------------------

    do iElem = 1, n
      uy(iElem) =  me%u0(2) * ( sin(coord(iElem,1)-me%x0(1)) &
        &                     * cos(coord(iElem,2)-me%x0(2)) )
    end do
  end function ic_tgv_uy_for