58 lines
1.4 KiB
Matlab
58 lines
1.4 KiB
Matlab
% Assume output from LevelEllipsoid has been read in
|
|
doinside=1;
|
|
red=[179,27,27]/255;
|
|
white=[1,1,1];
|
|
black=[0,0,0];
|
|
blue=[0,19,56]/100;
|
|
green=[9,45,27]/100;
|
|
gray=[0.9,0.9,0.9];
|
|
thick=1;
|
|
pltsize=[5.9 4.12];
|
|
set(gcf,'Units','pixels');
|
|
set(gcf,'Position',50+150*[0 0 pltsize]);
|
|
set(gcf,'Units','normalized');
|
|
hold off;
|
|
if ~doinside
|
|
fill([p{1}(:,1);0],[p{1}(:,2);0], gray, 'EdgeColor', 'none');
|
|
hold on
|
|
end
|
|
nq=size(q,2);
|
|
for i=1:nq;
|
|
plot(q{i}(:,1),q{i}(:,2), 'Color', green);
|
|
if i == 1, hold on; end
|
|
end
|
|
np=size(p,2);
|
|
for i=1:np
|
|
color = blue;
|
|
if i == 1
|
|
plot(p{i}(:,1),p{i}(:,2), 'Color', red, 'LineWidth', thick);
|
|
else
|
|
plot(p{i}(:,1),p{i}(:,2), 'Color', blue);
|
|
end
|
|
end
|
|
if doinside
|
|
nq=size(qa,2);
|
|
for i=1:nq;
|
|
plot(qa{i}(:,1),qa{i}(:,2), 'Color', green);
|
|
end
|
|
np=size(pa,2);
|
|
for i=1:np
|
|
plot(pa{i}(:,1),pa{i}(:,2), 'Color', blue);
|
|
end
|
|
plot([0, 0.6], [0, 0], 'Color', black, 'LineWidth', thick);
|
|
end
|
|
hold off
|
|
xlabel('R'); ylabel('Z');
|
|
axis equal;
|
|
axis([0,xmax,0,ymax]);
|
|
set(gcf,'PaperOrientation', 'landscape');
|
|
set(gcf,'PaperSize',pltsize);
|
|
set(gcf,'PaperPosition',[0 0 pltsize]);
|
|
set(gca, 'XTick',[0:3]);
|
|
set(gca, 'YTick',[0:2]);
|
|
set(gca, 'LooseInset',[0.07 0.09 0.03 0.02]);
|
|
ylabelh=get(gca,'ylabel');
|
|
set(ylabelh,'rotation',0);
|
|
set(ylabelh,'Position', get(ylabelh, 'Position') + [-0.1 0.2 0]);
|
|
print('-dsvg', ['normal-gravity-potential-', num2str(doinside), '.svg']);
|