%% velocity data of single defect where another defect is at the apex

v = [-.36, -.56, -.5, .25, 1.3, 2.73];

chi = [.4, .5, .6, .66666666667, .7, .75];

%% velocity data for single defect at origin

v = 1.0e-05*[0.3070, 0.2996, 0.2562, 0.1627, 0.0028, -0.2287, -0.5066, -0.7282];

chi = [0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8];

%% data for defect trajectories

chi = 0.1;

% data for x = 0, y = 40:
x0_40 = [-0.5,  1.5,  3.5,  5.5,  6.5,  8.5,  9.5, 10.5, 10.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 10.5,  8.5,  6.5, -0.5, -0.5, -0.5];
y0_40 = [39.5, 36.5, 34.5, 32.5, 30.5, 29.5, 27.5, 25.5, 23.5, 21.5, 20.5, 18.5, 16.5, 14.5, 12.5, 10.5,  8.5,  5.5, -0.5, -0.5, -0.5];

% data for x = 20, y = 40:
x20_40 = [19.5, 21.5, 22.5, 23.5, 24.5, 25.5, 26.5, 26.5, 27.5, 28.5, 29.5, 30.5, 31.5, 32.5, 33.5, 34.5, 34.5, 35.5, 36.5, 37.5, 38.5];
y20_40 = [39.5, 37.5, 35.5, 34.5, 33.5, 32.5, 31.5, 30.5, 29.5, 28.5, 28.5, 27.5, 27.5, 26.5, 26.5, 26.5, 25.5, 25.5, 25.5, 25.5, 25.5];

% data for x = 20, y = 40, free boundary conditions:
xFree20_40 = [19.5, 21.5, 22.5, 23.5, 24.5, 26.5, 27.5, 29.5, 31.5];
yFree20_40 = [39.5, 37.5, 35.5, 34.5, 33.5, 31.5, 30.5, 28.5, 26.5];

% data for x = 40, y = 40:
x40_40 = [39.5, 40.5, 41.5, 43.5, 44.5, 45.5, 46.5, 47.5, 49.5, 50.5, 51.5, 53.5, 54.5, 56.5, 58.5, 59.5, 61.5, 63.5, 65.5, 67.5, 69.5];
y40_40 = [39.5, 38.5, 37.5, 36.5, 35.5, 35.5, 34.5, 34.5, 34.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 33.5, 34.5];

% data for x = -10, y = 0:
xm10_0 = [-10.5, -9.5, -8.5, -8.5, -7.5, -6.5, -5.5, -4.5, -3.5, -1.5, -1.5];
ym10_0 = [-0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5];

% data for x = -10, y = 20:
xm10_20 = [-10.5, -8.5, -6.5, -4.5, -3.5, -1.5, -0.5, -0.5, -0.5, -0.5, -0.5];
ym10_20 = [19.5, 17.5, 15.5, 13.5, 10.5, 7.5, 0.5, -0.5, -0.5, -0.5, -0.5];

% data for x = -10, y = 20, free boundary conditions:
xFreem10_20 = [-10.5, -8.5, -6.5, -4.5, -1.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5];
yFreem10_20 = [19.5, 17.5, 15.5, 13.5, 10.5, 7.5, -0.5, -0.5, -0.5, -0.5, -0.5];

% data for x = -10, y = 30:
xm10_30 = [-10.5,  -7.5,  -4.5,  -2.5,  -0.5,   0.5,   1.5,  -0.5,  -0.5, -0.5,  -0.5,  -0.5,  -0.5,  -0.5,  -0.5,  -0.5,  -0.5,  -0.5, -0.5,  -0.5,  -0.5];
ym10_30 = [29.5, 26.5, 23.5, 20.5, 17.5, 13.5,  7.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5];

% data for x = -10, y = 40:
xm10_40 = [-10.5, -7.5, -5.5, -2.5, -1.5, 0.5, 2.5, 3.5, 4.5, 5.5, 5.5, 6.5, 5.5, 4.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5];
ym10_40 = [39.5, 37.5, 34.5, 32.5, 30.5, 28.5, 26.5, 24.5, 22.5, 19.5, 17.5, 14.5, 11.5, 7.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5, -0.5];
