function contour_plot ( )

%*****************************************************************************80
%
%% volcano_fill_contour() makes a filled contour plot of volcano height.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    03 May 2019
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'contour_plot():\n' );
  fprintf ( 1, '  Make a filled color contour plot of volcano height.\n' );
%
%  Get the spatial data.
%
  [ ~, ~, L, ~, ~, ~, ~ ] = mhd_parameters ( );

  xvec = linspace ( 0.0, L, 21 );
  yvec = linspace ( -1.0, 1.0, 11 );
  
  [ X, Y ] = meshgrid ( xvec, yvec );
%
%  Evaluate the solution.
%
  [ U, V, P, B1, B2 ] = mhd_solution ( X, Y );
%
%  Plot the pressure.
%
  contourf ( X, Y, P );

  title ( 'pressure contours' )

  filename = 'pressure_contour.png';
  print ( '-dpng', filename )
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Graphics saved as "%s"\n', filename );

  return
end

