function [ U, V ] = spiral_initial_condition ( X, Y, alpha )

%*****************************************************************************80
%
%% spiral_initial_condition(): initial values for spiral reaction-diffusion equation.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    16 December 2020
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    real X(nr,nc), Y(nr,nc): the node coordinates.
%
%    real alpha: a problem parameter.
%
%  Output:
%
%    real U(nr,nc), V(nr,nc): the node velocities.
%
  U = ( 40.0 <= X );

  V = 0.5 * alpha * ( 40.0 <= Y );

  return
end

