function volcano_surface ( )

%*****************************************************************************80
%
%% volcano_surface makes a 3d surface plot of height data near a volcano.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    19 April 2019
%
%  Author:
%
%    John Burkardt
%
  timestamp ( );
  fprintf ( 1, '\n' );
  fprintf ( 1, 'volcano_surface:\n' );
  fprintf ( 1, '  MATLAB/Octave version %s:\n', version ( ) );
  fprintf ( 1, '  Make a surface plot of volcano height.\n' );

  Z = csvread ( 'volcano_data.csv' );

  [ m, n ] = size ( Z );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Z values stored as %d by %d array\n', m, n );

  surf ( Z );
  title ( 'Surface plot of volcano height', 'Fontsize', 16 );

  filename = 'volcano_surface.png';
  print ( '-dpng', filename )
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Graphics saved as "%s"\n', filename );
%
%  Terminate.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'volcano_surface:\n' );
  fprintf ( 1, '  Normal end of execution.\n' );
  timestamp ( );

  return
end

