function [ X, Y ] = xy_shear ( nr, nc )

%*****************************************************************************80
%
%% xy_shear returns node coordinates for the shear flow.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    16 July 2020
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer nr, nc: the number of rows and columns.
%
%  Output:
%
%    real X(nr,nc), Y(nr,nc): the node coordinates.
%
  x = linspace (  0.0, +1.0, nc );
  y = linspace (  0.0, +1.0, nr );

  [ X, Y ] = ndgrid ( x, y );

  X = rot90 ( X );
  Y = rot90 ( Y );

  return
end

