enderton/Diagnostics/DiagPlotMakePlot.m

% DiagPlotMakePlot is called by DiagPlot and cannot be used seperately.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                                 Make plot                               %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if isequal(cmp,'Dif')
    plotdata = data{inrow}{1} - data{inrow}{2};
else
    plotdata = data{inrow}{incol};
end


% Make gridded plot.  This should be fixed up to allow for non-constant
% height intervals, imagesc makes everything equally spaced in the vertical
% (and in the horizontal for that matter).
if isequal(pst,'Grd')
    if isCS
        merccube(xax{inrow}{incol},yax{inrow}{incol},plotdata);
    else
        imagesc(xax{inrow}{incol},yax{inrow}{incol},plotdata);
    end
    
% Make interpolated plot.
elseif isequal(pst,'Int')
    if isCS
        merccube(xax{inrow}{incol},yax{inrow}{incol},plotdata);
        shading interp;
    else
        pcolor(xax{inrow}{incol},yax{inrow}{incol},plotdata);
        shading interp;
    end
    
% Make contour (non-filled) plot.  Let MATLAB determine the contour
% intervals on a differencing plot.
elseif isequal(pst,'Con')
    if isequal(cmp,'Dif')
        [cs,h] = contour(xax{inrow}{incol},yax{inrow}{incol},plotdata);
    else
        [cs,h]=contour(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint);
        set(h,'linewidth',linewidth)
    end
    if UseConLabel
        clabel(cs,h,'fontsize',fs_clabel,'rotation',0);
    end
    
% Make contour (filled) plot.  Let MATLAB determine the contour
% intervals on a differencing plot.
elseif isequal(pst,'Cnf')
    if isequal(cmp,'Dif')
        [cs,h] = contourf(xax{inrow}{incol},yax{inrow}{incol},plotdata);
    else
        [cs,h] = contourf(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint);
    end
    if UseCnfLabel
        clabel(cs,h,'fontsize',fs_clabel,'rotation',0);
    end
    
% Make line plot.  If this is an overlaying comparison, cycle through the
% appropriate data of which to overlay.
elseif isequal(pst,'Lin')
    if ismember(cmp,{'OvC','OvE','OvF'})
        for intrl = 1:ntrl
            plot(xax{inrow}{intrl},data{inrow}{intrl},linecolors{intrl},'linewidth',linewidth);
        end
    else
        plot(xax{inrow}{incol},plotdata,linecolors{1},'linewidth',linewidth);
    end
    
% Unknown plot style encountered.
else
    error(['Oh dear!  Undefined PlotStyle:  ',pst]);
end
1	% DiagPlotMakePlot is called by DiagPlot and cannot be used seperately.
2
3
4	%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5	% Make plot %
6	%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
7	if isequal(cmp,'Dif')
8	plotdata = data{inrow}{1} - data{inrow}{2};
9	else
10	plotdata = data{inrow}{incol};
11	end
12
13
14	% Make gridded plot. This should be fixed up to allow for non-constant
15	% height intervals, imagesc makes everything equally spaced in the vertical
16	% (and in the horizontal for that matter).
17	if isequal(pst,'Grd')
18	if isCS
19	merccube(xax{inrow}{incol},yax{inrow}{incol},plotdata);
20	else
21	imagesc(xax{inrow}{incol},yax{inrow}{incol},plotdata);
22	end
23
24	% Make interpolated plot.
25	elseif isequal(pst,'Int')
26	if isCS
27	merccube(xax{inrow}{incol},yax{inrow}{incol},plotdata);
28	shading interp;
29	else
30	pcolor(xax{inrow}{incol},yax{inrow}{incol},plotdata);
31	shading interp;
32	end
33
34	% Make contour (non-filled) plot. Let MATLAB determine the contour
35	% intervals on a differencing plot.
36	elseif isequal(pst,'Con')
37	if isequal(cmp,'Dif')
38	[cs,h] = contour(xax{inrow}{incol},yax{inrow}{incol},plotdata);
39	else
40	[cs,h]=contour(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint);
41	set(h,'linewidth',linewidth)
42	end
43	if UseConLabel
44	clabel(cs,h,'fontsize',fs_clabel,'rotation',0);
45	end
46
47	% Make contour (filled) plot. Let MATLAB determine the contour
48	% intervals on a differencing plot.
49	elseif isequal(pst,'Cnf')
50	if isequal(cmp,'Dif')
51	[cs,h] = contourf(xax{inrow}{incol},yax{inrow}{incol},plotdata);
52	else
53	[cs,h] = contourf(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint);
54	end
55	if UseCnfLabel
56	clabel(cs,h,'fontsize',fs_clabel,'rotation',0);
57	end
58
59	% Make line plot. If this is an overlaying comparison, cycle through the
60	% appropriate data of which to overlay.
61	elseif isequal(pst,'Lin')
62	if ismember(cmp,{'OvC','OvE','OvF'})
63	for intrl = 1:ntrl
64	plot(xax{inrow}{intrl},data{inrow}{intrl},linecolors{intrl},'linewidth',linewidth);
65	end
66	else
67	plot(xax{inrow}{incol},plotdata,linecolors{1},'linewidth',linewidth);
68	end
69
70	% Unknown plot style encountered.
71	else
72	error(['Oh dear! Undefined PlotStyle: ',pst]);
73	end