1 function [vorts] = vtxTrajectory(start, step, fin, initial,threeD,
dx,
dimSize,
dt)
2 %vtxTrajectory trajectory
plots of 2D vortex
data.
3 % vtxTrajectory takes the output of the postprocessed vortex positions
4 % from the matlab script
as loaded
by loadVtx
5 % Start step and fin are the initial, stepsize and
final datasets
6 % initial indicates
if you wish to highlight the
8 %
dx is the stepsize along one dimension and assumes the same over
x,
y 10 %
dt is the
time increment
for the spacetime diagram.
11 % Note: Since the processed
data assumes an initial starting
position, the
16 %Load vortices and convert
struct to
cell format
17 vorts = struct2cell(loadVtx(start,step,fin,0));
22 f0=csvread(['vort_ord_',int2str(start),'.csv'],1,0);
23 p=scatter((f0(:,1)-
dimSize/2)*
dx,(f0(:,2)-
dimSize/2)*
dx,40,[0. 0. 0.],'filled','
MarkerEdgeColor',[0. 0. 0.],'MarkerFaceColor',[0.4,0.4,0.4]); hold
on;
26 plot(([
vorts{1,
ii,:}]-
dimSize/2)*
dx,([
vorts{2,
ii,:}]-
dimSize/2)*
dx,'color',rand(3,1),'
LineWidth',2);
30 set(gca,
'TickLabelInterpreter',
'latex');
31 set(gca,
'FontName',
'Latin Modern Roman',
'FontSize',36);
axis square;
32 xlabel(
' $x$ (m)',
'Interpreter',
'latex');
ylabel(
' $y$ (m)',
'Interpreter',
'latex');
33 print(
'-dpng',
'-r300',strcat(
'Trajectory.png',
''));
35 if threeD==1% Generate 2+1 spacetime
plot of vortices
40 set(gca,
'TickLabelInterpreter',
'latex');
41 set(gca,
'FontName',
'Latin Modern Roman',
'FontSize',22);
42 xlabel(
'$x$ (m)',
'Interpreter',
'latex');
ylabel(
'$y$ (m)',
'Interpreter',
'latex');zlabel(
'$t$ (s)',
'Interpreter',
'latex');
43 print(
'-dpng',
'-r300',strcat(
'Trajectory3d.png',
''));
end if(length(DT.vertexAttachments{ii})==5) if plotit plot(x(ii)
% Imaginary time unitary evolution operators for position
stepsize and final datasets % initial indicates if you wish to highlight the % threeD also plots aspacetime(XYT) diagram in 3D. % dx is the stepsize along one dimension and assumes the same over x
ylabel(' $y$(m)', 'Interpreter', 'latex') % Defect marker size MarkerSize
stepsize and final datasets % initial indicates if you wish to highlight the % threeD also plots y % dimSize is the length of the and assumes the % data at is named while postprocessed data is vort_ord_xxx % Also
stepsize and final datasets % initial indicates if you wish to highlight the % threeD also plots y % dimSize is the length of the dimensions
and extended by(Dr?) Lee James O 'Riordan. % Ek
and print the results to file Essentially wraps % voronoi2dCellColour and saves the resulting plots
stepsize and final datasets % initial indicates if you wish to highlight the % threeD also plots y % dimSize is the length of the and assumes x
stepsize and final datasets % initial indicates if you wish to highlight the % threeD also plots y % dimSize is the length of the and assumes the % data at is named vort_arr_0
and assumes too high and very sparse Start at
and structures them % into Matlab structs for Vortex position and index % start step and fin are the initial step and final % dataset to load for the vortex postions % lostCheck turns on the ability to test for vortex losses and ensure % that these vortices do not contribute to the dataset % The Data should be in the format of [X Y Sign UID IsLost] % vorts
% Indexing needs to % be modified if you wish to use the ordered data sets % Calculate the Voronoi diagram of the resulting data
% axis([1e4 1e7 5e-18 1e-10])
and colour the % cells with the value of the orientational correlations defined at each % or with the area of the cell(col=0). % x
% If at least one of the indices is
Latin Modern latex print('-dpng','-r300', ['./Comp_CBAR_', int2str(iii),'.png'])
xlabel(' $x$(m)', 'Interpreter', 'latex')
end % idx is for idx is for type for ii