%  setplot3.m
%  called in plotclaw3 before plotting to set various parameters

setviews  % set viewpoints so that view(xSlice), for example, can be used.

PlotType = 3;                % type of plot to produce:
			     % 1 = pcolor on slices (with optional contours)
			     % 2 = contour lines in 3d on transparent slices
			     % 3 = Schlieren plot on slices
			     % 4 = scatter plot of q vs. r
			     % 5 = isosurface plot (at given levels)

mq = 1;                      % which component of q to plot
UserVariable = 0;            % set to 1 to specify a user-defined variable
UserVariableFile = ' ';      % name of m-file mapping data to q
MappedGrid = 0;              % set to 1 if mapc2p.m exists for nonuniform grid
MaxFrames = 1000;            % max number of frames to loop over
MaxLevels = 6;               % max number of AMR levels

PlotData =  [1 1 1 0 0 0];       % Data on refinement level k is plotted only 
			         % if k'th component is nonzero
PlotGrid =  [0 0 0 0 0 0];       % Plot grid lines on each level?
PlotGridEdges =  [1 0 0 0 0 0];  % Plot edges of patches of each grid at
                                 % this level on slices?
PlotCubeEdges = [0 0 0 0 0 0];   % Plot edges of cube of refinement patch at
                                 % this level?


% ContourValues is a vector of contour lines that can be used with
% PlotType = 1,2,3.   Empty ==> no contour lines drawn
  ContourValues = [];     

% The next three parameters are vectors of x,y,z coordinates of 2d slices
% to be displayed for PlotType = 1,2,3.  Empty ==> no slices in that direction.
  xSliceCoords = [];     
  ySliceCoords = 0.5;
  zSliceCoords = 0.;

% For PlotType = 4 (Scatter plot)
% plot q(r) vs. r = sqrt((x-x0)^2 + (y-y0)^2 + (z-z0)^2);
% Use symbol PlotStyle(k) at refinement level k.
  x0 = 0.5;
  y0 = 0.5;
  z0 = 0.5;
  PlotStyle = setplotstyle('o','x','.','s','v','^');

% For PlotType = 5 (Isosurface plots)
% Note:  Lengths of SurfTransparency and SurfColors must greater than or
%        equal to the length of SurfValues..

  %IsosurfValues    =  [2.4 0.5 0.9];     % Plot surfaces at q = surfValue(i).
  IsosurfValues    =  [2.4 0.5];     % Plot surfaces at q = surfValue(i).

  IsosurfAlphas    =  [0.2   1  ];     % Transparency of each surface 
  %IsosurfAlphas    =  [0.2   1  0.5];     % Transparency of each surface 
                                          % (0=clear; 1=opaque)
                                          % NOTE: Your system must be able to 
                                          % use the OpenGL Renderer.

  IsosurfColors    = strvcat('b','r');      % Colors for each surface.
  %IsosurfColors    = strvcat('b','r','c');      % Colors for each surface.