function grids_display ( input_filename1, input_filename2, rmin, rmax )

%*****************************************************************************80
%
%% grids_display() displays a pair of 1D, 2D or 3D grids.
%
%  Discussion:
%
%    The points of the first grid will be blue, the second red.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    05 February 2014
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'grids_display():\n' );
  fprintf ( 1, '  MATLAB version\n' );
  fprintf ( 1, '  Display a pair of 1D, 2D or 3D grids.\n' );
%
%  First command line argument is file name1.
%
  if ( nargin < 1 )
    fprintf ( 1, '\n' );
    fprintf ( 1, 'GRIDS_DISPLAY:\n' );
    input_filename1 = input ( '  Enter the first filename:  ' );
  end
%
%  Second command line argument is file name2.
%
  if ( nargin < 2 )
    fprintf ( 1, '\n' );
    fprintf ( 1, 'GRIDS_DISPLAY:\n' );
    input_filename2 = input ( '  Enter the second filename:  ' );
  end
%
%  Third command line argument is the lower range.
%
  if ( nargin < 3 )
    fprintf ( 1, '\n' );
    fprintf ( 1, 'GRIDS_DISPLAY:\n' );
    rmin = input ( '  Enter the minimum value of the range:  ' );
  elseif ( ischar ( rmin ) )
    rmin = str2num ( rmin );
  end
%
%  Fourth command line argument is the upper range.
%
  if ( nargin < 4 )
    fprintf ( 1, '\n' );
    fprintf ( 1, 'GRIDS_DISPLAY:\n' );
    rmax = input ( '  Enter the maximum value of the range:  ' );
  elseif ( ischar ( rmax ) )
    rmax = str2num ( rmax );
  end
%
%  Load the data.
%  
  points1 = load ( input_filename1 )';

  dim_num1 = size ( points1, 1 );
  point_num1 = size ( points1, 2 );

  points2 = load ( input_filename2 )';

  dim_num2 = size ( points2, 1 );
  point_num2 = size ( points2, 2 );

  if ( dim_num1 ~= dim_num2 ) 
    fprintf ( 1, '\n' );
    fprintf ( 1, 'GRIDS_DISPLAY - Fatal error!\n' );
    fprintf ( 1, '  The spatial dimensions differ.\n' );
    error ( 'GRIDS_DISPLAY - Fatal error' );
  else
    dim_num = dim_num1;
  end
%
%  1D Plot
%
  if ( dim_num == 1 )
%
%  Clear the graphics frame;
%
    clf ( );
%
%  We have to name the axes in order to control the grid.
%
    axes_handle = axes;
%
%  Plot the points.
%
    points1(2,1:point_num1) = - 0.25 * ones(1,point_num1);
    h1 = scatter ( points1(1,:), points1(2,:), 80, 'r', 'filled' );
    hold ( 'on' );
    points2(2,1:point_num2) = 0.25 * ones(1,point_num2);
    h2 = scatter ( points2(1,:), points2(2,:), 80, 'b', 'filled' );
%
%  Request grid lines.
%
    grid ( 'on' );
%
%  Specify the location of the grid lines, and suppress labeling.
%
    rdelt = ( rmax - rmin ) / 10.0;
    ticks = rmin : rdelt : rmax;
  
    set ( axes_handle, 'xtick', ticks );
%   set ( axes_handle, 'xticklabel', [] );
    set ( axes_handle, 'ytick', ticks );
    set ( axes_handle, 'yticklabel', [] );
%
%  Make the plotting region slightly bigger than the data.
%
    range = rmax - rmin;
    rmin1 = rmin - 0.05 * range;
    rmax1 = rmax + 0.05 * range;
    axis ( [ rmin1, rmax1, rmin1, rmax1 ] );

    xlabel ( '<--- X axis --->', 'Fontsize', 16 );
%
%  Title
%
    title_string = strcat ( input_filename1, ' + ' );
    title_string = strcat ( title_string, input_filename2 );
    title_string = s_escape_tex ( title_string );
    title ( title_string, 'Fontsize', 24 );
%
%  Make an output file name.
%
    output_filename = 'grids_1d.png';
    print ( '-dpng', output_filename );
    fprintf ( 1, '\n' );
    fprintf ( 1, '  Graphics saved as "%s"\n', output_filename );

    hold ( 'off' );
%
%  2D Plot
%
  elseif ( dim_num == 2 )
%
%  Clear the graphics frame;
%
    clf ( );
%
%  We have to name the axes in order to control the grid.
%
    axes_handle = axes;
%
%  Plot the points.
%  TEMPORARY MOD FOR CLAYTON, GRAY + BLACK
%
    hold ( 'on' );

    if ( false )
      h1 = plot ( points2(1,:), points2(2,:), 'k.', 'Markersize', 15 );
      h2 = plot ( points1(1,:), points1(2,:), 'k.', 'Markersize', 25 );
    else
      h1 = plot ( points2(1,:), points2(2,:), 'ro', 'Markersize', 15 );
      h2 = plot ( points1(1,:), points1(2,:), 'b.', 'Markersize', 25 );
    end
%
%  Force the plotting region to be square, not rectangular.
%
    axis ( 'square' );
%
%  Request grid lines.
%
    grid ( 'on' );
%
%  Specify the location of the grid lines, and suppress labeling.
%
    rdelt = ( rmax - rmin ) / 10.0;
    ticks = rmin : rdelt : rmax;
  
    set ( axes_handle, 'xtick', ticks );
    set ( axes_handle, 'xticklabel', [] );
    set ( axes_handle, 'ytick', ticks );
    set ( axes_handle, 'yticklabel', [] );
%
%  Make the plotting region slightly bigger than the data.
%
    range = rmax - rmin;
    rmin1 = rmin - 0.05 * range;
    rmax1 = rmax + 0.05 * range;
    axis ( [ rmin1, rmax1, rmin1, rmax1 ] );

    if ( false )
      xlabel ( '--X axis--', 'Fontsize', 16 );
      ylabel ( '--Y axis--', 'Fontsize', 16 );
    end
%
%  Title
%
    if ( false )
      title_string = strcat ( input_filename1, ' + ' );
      title_string = strcat ( title_string, input_filename2 );
      title_string = s_escape_tex ( title_string );
      title ( title_string, 'Fontsize', 24 );
    end
%
%  Make an output file name.
%
    output_filename = 'grids_2d.png';
    print ( '-dpng', output_filename );
    fprintf ( 1, '\n' );
    fprintf ( 1, '  Graphics saved as "%s"\n', output_filename );

    hold ( 'off' );
%
%  3D plot
%
  elseif ( dim_num == 3 )
%
%  Clear the graphics frame;
%
    clf ( );
%
%  We have to name the axes in order to control the grid.
%
    axes_handle = axes;
%
%  Plot the points.
%
    hold ( 'on' );

    handle1 = scatter3 ( points1(1,:), points1(2,:), points1(3,:), ...
      'r', 'filled' );
    handle2 = scatter3 ( points2(1,:), points2(2,:), points2(3,:), ...
      'b', 'filled' );
%
%  Force the plotting region to be square, not rectangular.
%
    axis ( 'square' );
%
%  Request grid lines.
%
    grid ( 'on' );
%
%  Specify the location of the grid lines, and suppress labeling.
%
    rdelt = ( rmax - rmin ) / 10.0;
    ticks = rmin : rdelt : rmax;
  
    set ( axes_handle, 'xtick', ticks );
    set ( axes_handle, 'xticklabel', [] );
    set ( axes_handle, 'ytick', ticks );
    set ( axes_handle, 'yticklabel', [] );
    set ( axes_handle, 'ztick', ticks );
    set ( axes_handle, 'zticklabel', [] );
%
%  Make the plotting region slightly bigger than the data.
%
    range = rmax - rmin;
    rmin1 = rmin - 0.05 * range;
    rmax1 = rmax + 0.05 * range;
    axis ( [ rmin1, rmax1, rmin1, rmax1 ] );

    xlabel ( '--X axis--', 'Fontsize', 16 );
    ylabel ( '--Y axis--', 'Fontsize', 16 );
    zlabel ( '--Z axis--', 'Fontsize', 16 );
%
%  Title
%
    title_string = strcat ( input_filename1, ' + ' );
    title_string = strcat ( title_string, input_filename2 );
    title_string = s_escape_tex ( title_string );
    title ( title_string, 'Fontsize', 24 );
    
    view ( -10, 15 );
%
%  Make an output file name.
%
    output_filename = 'grids_3d.png';
    print ( '-dpng', output_filename );
    fprintf ( 1, '\n' );
    fprintf ( 1, '  Graphics saved as "%s"\n', output_filename );

    hold ( 'off' );

  end
%
%  Terminate.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'grids_display():\n' );
  fprintf ( 1, '  Normal end of execution.\n' );

  return
end
function s2 = s_escape_tex ( s1 )

%*****************************************************************************80
%
%% s_escape_tex() de-escapes TeX escape sequences.
%
%  Discussion:
%
%    In particular, every occurrence of the characters '\', '_',
%    '^', '{' and '}' will be replaced by '\\', '\_', '\^',
%    '\{' and '\}'.  A TeX interpreter, on seeing these character
%    strings, is then likely to return the original characters.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    19 January 2007
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    string S1, the string to be de-escaped.
%
%  Output:
%
%    string S2, a copy of the string, modified to avoid TeX escapes.
%
  s1_length = length ( s1 );

  s1_pos = 0;
  s2_pos = 0;
  s2 = [];

  while ( s1_pos < s1_length )

    s1_pos = s1_pos + 1;

    if ( s1(s1_pos) == '\' | ...
         s1(s1_pos) == '_' | ...
         s1(s1_pos) == '^' | ...
         s1(s1_pos) == '{' | ...
         s1(s1_pos) == '}' )
      s2_pos = s2_pos + 1;
      s2 = strcat ( s2, '\' );
    end

    s2_pos = s2_pos + 1;
    s2 = strcat ( s2, s1(s1_pos) );

  end

  return
end
function timestamp ( )

%*****************************************************************************80
%
%% timestamp() prints the current YMDHMS date as a timestamp.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    14 February 2003
%
%  Author:
%
%    John Burkardt
%
  t = now;
  c = datevec ( t );
  s = datestr ( c, 0 );
  fprintf ( 1, '%s\n', s );

  return
end
 
