fig3_article.m 2.72 KB
clear all;
close all;


dossier=['fig1_1'];
path=['/home/bauj0001/projets/WIM/output/',dossier,'/',dossier,'.nc'];


Fsize=ncread(path,'floe size');
strain1=load('/home/bauj0001/projets/WIM/output/strain1.dat');
beta1=load('/home/bauj0001/projets/WIM/output/beta1.dat');
Q1=load('/home/bauj0001/projets/WIM/output/Q1.dat');

strain2=load('/home/bauj0001/projets/WIM/output/strain2.dat');
beta2=load('/home/bauj0001/projets/WIM/output/beta2.dat');
Q2=load('/home/bauj0001/projets/WIM/output/Q2.dat');

strain3=load('/home/bauj0001/projets/WIM/output/strain3.dat');
beta3=load('/home/bauj0001/projets/WIM/output/beta3.dat');
Q3=load('/home/bauj0001/projets/WIM/output/Q3.dat');

Lmin=pi*0.5*((5e6*2^3)/(3*10*(1-0.3^2)))^0.25;
Lmin=Lmin*2;
wavelength=Fsize*2;


fig=figure;

fig.Units='inches';
fig.Position=[7.4167 2.0833 2.9583 5.5312];
fig.OuterPosition=[7.4167 2.0833 2.9583 6.5312];

fig.PaperUnits = 'inches';
fig.PaperPosition = [0 0 2.9583 5.5312];
fig.PaperPositionMode = 'manual';
fig.PaperOrientation='landscape';
fig.PaperSize=[2.9583 5.5312];




subplot(3,1,1)
plot(wavelength(find(wavelength>Lmin)),strain1(find(wavelength>Lmin)),'color',[0.1 0.1 0.1],'linewidth',1.5)
hold on
plot(wavelength(find(wavelength>Lmin)),strain2(find(wavelength>Lmin)),'color',[0.5 0.5 0.5],'linewidth',1.5)
plot(wavelength(find(wavelength>Lmin)),strain3(find(wavelength>Lmin)),'color',[0.8 0.8 0.8],'linewidth',1.5)
a=ylabel('$E_s$');
b=xlabel('$\lambda_{\mathrm{ice}}$ [m]');
set(a,'interpreter','latex')
set(b,'interpreter','latex')
c=legend('Hs=0.1 m','Hs=0.5 m ','Hs=1 m');
set(c,'box','off')
hold on
plot([Lmin Lmin],ylim,'color',[0.85 0.325 0.098])
e=plot(wavelength(find(wavelength<Lmin)),strain1(find(wavelength<Lmin)),'color',[0.1 0.1 0.1],'linewidth',1.5);
set(e,'linestyle',':')
f=plot(wavelength(find(wavelength<Lmin)),strain2(find(wavelength<Lmin)),'color',[0.5 0.5 0.5],'linewidth',1.5);
set(f,'linestyle',':')
g=plot(wavelength(find(wavelength<Lmin)),strain3(find(wavelength<Lmin)),'color',[0.8 0.8 0.8],'linewidth',1.5);
set(g,'linestyle',':')
j=text(55,2e-3,'$\lambda_{MIN}$');
set(j,'interpreter','latex','rotation',90)


subplot(3,1,2)
plot(Fsize,Q1,'color',[0.1 0.1 0.1],'linewidth',1.5)
hold on
plot(Fsize,Q2,'color',[0.5 0.5 0.5],'linewidth',1.5)
plot(Fsize,Q3,'color',[0.8 0.8 0.8],'linewidth',1.5)
a=ylabel('$Q$');
b=xlabel('Floe size $l$ [m]');
set(a,'interpreter','latex')
set(b,'interpreter','latex')

subplot(3,1,3)
plot(Fsize,beta1(end,:),'color',[0.1 0.1 0.1],'linewidth',1.5)
hold on
plot(Fsize,beta2(end,:),'color',[0.5 0.5 0.5],'linewidth',1.5)
plot(Fsize,beta3(end,:),'color',[0.8 0.8 0.8],'linewidth',1.5)
a=ylabel('$\beta(l=400,l'')$');
b=xlabel('Floe size $l$ [m]');
set(a,'interpreter','latex')
set(b,'interpreter','latex')

saveas(fig,'test','pdf')