fig5_article.m 666 Bytes
Newer Older
Jérémy Baudry's avatar
Jérémy Baudry committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42
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