fig5_article.m 666 Bytes
clear all;
close all;



fig=figure;

fig.Units='inches';
fig.Position=[7.0729 5.6146 5.8333 4.3750];
fig.OuterPosition=[7.0729 5.6146 5.8333 5.3750];

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



for i=1:100
    
    
    dossier=['thickhs',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));
  
 
  
   
end

Lmiz=reshape(Lmiz,10,10);
for i=1:10
    plot(hs,Lmiz(:,i))
    hold on
end