% setplot3.m % called in plotclaw3 before plotting to set various parameters PlotType = 1; % 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 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 = linspace(-0.005,0.005,20); % 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 = 0.0; ySliceCoords = 0.0; zSliceCoords = []; % 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.0; y0 = 0.0; z0 = 1.0; PlotStyle = setplotstyle('ro');