fig_strain.m 1.55 KB
clear all;
close all;




kice=load('kice.dat');
dossier=['wavenumber'];
path=['/home/bauj0001/projets/WIM/output/',dossier,'/',dossier,'.nc'];
hcat=ncread(path,'thickness categories');
om=ncread(path,'omega');


fig=figure;

fig.Units='inches';
fig.Position=[7.0312 7.9479 5.3125 2.0417];
fig.OuterPosition=[7.0312 7.9479 5.3125 3.0417];

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



subplot(1,2,1)
W=9.81*kice(:,20)/om(20);
strain=hcat*0.5.*kice(:,20).^2.*W;
plot(hcat,strain,'color',[0.1 0.1 0.1],'linewidth',1.5)
wl(1)=2*pi/(om(20)^2/9.81);
hold on

W=9.81*kice(:,50)/om(50);
strain=hcat*0.5.*kice(:,50).^2.*W;
plot(hcat,strain,'color',[0.5 0.5 0.5],'linewidth',1.5)
wl(2)=2*pi/(om(50)^2/9.81);
hold on

W=9.81*kice(:,100)/om(100);
strain=hcat*0.5.*kice(:,100).^2.*W;
plot(hcat,strain,'color',[0.8 0.8 0.8],'linewidth',1.5)
wl(3)=2*pi/(om(100)^2/9.81);
xlabel('Ice thickness [m]','interpreter','latex')
ylabel('E','interpreter','latex')
legende=['$\lambda_{W}=$ 225m';'$\lambda_{W}=$ 85m ';'$\lambda_{W}=$ 30m '];
l=legend(legende);
set(l,'interpreter','latex')
set(l,'box','off')


subplot(1,2,2)
wl2=2*pi./kice(:,20);
plot(hcat,wl2,'color',[0.1 0.1 0.1],'linewidth',1.5)
hold on

wl2=2*pi./kice(:,50);
plot(hcat,wl2,'color',[0.5 0.5 0.5],'linewidth',1.5)
hold on

wl2=2*pi./kice(:,100);
plot(hcat,wl2,'color',[0.8 0.8 0.8],'linewidth',1.5)
hold on

xlabel('Ice thickness [m]','interpreter','latex')
y2=ylabel('$\lambda_{ice}$ [m]','interpreter','latex');