fig8_article.m 1.32 KB
clear all;
close all;


for i=1:100
    
    
    dossier=['pc_ec',num2str(i-1)];
   path=['/home/bauj0001/projets/WIM/output/',dossier,'/',dossier,'.nc'];
   Dave=ncread(path,'Dave');
   
   x=ncread(path,'x_axis');

   Lmiz(i)=length(find(Dave>0 & Dave<396))*(x(2)-x(1));
  
 i
   
end

Lmiz2=reshape(Lmiz,10,10);
Ec=[3e-5 4e-5 5e-5 6e-5 7e-5 8e-5 9e-5 1e-4 2e-4 3e-4];

legende=['$P_c=0.10$';'$P_c=0.15$';'$P_c=0.20$';'$P_c=0.25$';'$P_c=0.30$'...
    ;'$P_c=0.35$';'$P_c=0.40$';'$P_c=0.45$';'$P_c=0.50$';'$P_c=0.55$'];
legende=cellstr(legende);

cmap=jet(10);
M=['*','o','+','*','o','+','*','o','+','*'];
    


fig=figure;

fig.Units='inches';
fig.Position=[7.0729 6.7396 4.2708 2.600];
fig.OuterPosition=[7.0729 6.7396 4.2708 3.600];

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


for i=1:10
 
    
    plot(Ec,Lmiz2(i,:),'color',cmap(i,:),'Marker',M(i));
    hold on
   
end  
  xlabel('$\varepsilon_c$','interpreter','latex')
  ylabel('$L_{MIZ}$ [km]','interpreter','latex')
  a=gca;
%   a.Position(2)=a.Position(2)+0.2;
%   a.Position(3)=a.Position(3)-0.2;
%   a.Position(4)=a.Position(4)-0.2;
  
  
  l=legend(char(legende));
  set(l,'interpreter','latex','Box','off','location','eastoutside');