function [ U0, V0 ] = spiral_initial_condition ( X, Y )

%*****************************************************************************80
%
%% spiral_initial_condition(): initial values for spiral reaction-diffusion PDE.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    05 February 2021
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    real X(nr,nc), Y(nr,nc): the node coordinates.
%
%  Output:
%
%    real U(nr,nc), V(nr,nc): the node velocities.
%
  [ alpha, beta, delta, epsilon ] = spiral_parameters ( );

  U0 = ( 40.0 <= X );

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

  return
end

