% draw clength planting paper

close all;

datadir = '../data/';
Dstr = '1.5';
Ddble = str2double(Dstr);
dataname = strcat('rhoxiD', Dstr, '_NN.dat');

psize = 18; % font size
msize = 9; % marker size

data = dlmread(strcat(datadir,dataname), '', 1, 0);
xafit = data(data(:,1)>30, 1);
yafit = log(data(data(:,1)>30, 3));
xafit = [xafit, ones(numel(xafit), 1)];
[b, bint] = regress(yafit, xafit);
m1 = b(1);
b1 = b(2);
errb = (bint(2,2)-bint(2,1))/2;
xregs = linspace(1,200);
yregs = xregs*m1+b1;

plot(data(:,1), data(:,3));
hold on;
plot(xregs, exp(yregs), 'r--');
plot(29.83, 5.35, '*b', 'markersize', msize, 'linewidth', 1.5);
set(gca, 'yscale', 'log');
xlim([0, 40]);
xticks(0:10:40);
xlabel('$\beta F d$', 'interpreter', 'latex');
ylabel('$\xi_y$', 'interpreter', 'latex');
set(gca, 'fontname', 'times new roman', 'fontsize', psize);
%% insetA
insetA = axes('Position', [0.2, 0.62, 0.3, 0.28]);

plot(data(:,1), data(:,3));
hold on;
plot(xregs, exp(yregs), 'r--');

set(insetA, 'yscale', 'log');
xticks([0, 100, 200]);
yticks([1e0, 1e5, 1e10]);
box off;

set(gca, 'fontname', 'times new roman', 'fontsize', 16);

%% insetB
insetB = axes('Position', [0.595, 0.25, 0.3, 0.28]);

plot(data(:,2), data(:,4));
hold on;
plot(1.0526, 0.43, '*b', 'markersize', 9, 'linewidth', 1.25);
xlim([0.8, 1.14]);
ylim([0.1, 0.5]);
yticks([0.1, 0.3, 0.5]);
xlabel('$\rho$', 'interpreter', 'latex');
ylabel('$\xi_{\delta x}$', 'interpreter', 'latex');

set(gca, 'fontname', 'times new roman', 'fontsize', 16);
annotation('textbox',[0.2, 0.82, 0.1, 0.1], 'String', '(a)', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);
annotation('textbox',[0.59, 0.44, 0.1, 0.1], 'String', '(b)', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);

print('../figures/clength.eps', '-depsc');