sensitivity_Pc.m 881 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];


dossier=['test1';'test2';'test3';'test4'];
dossier=cellstr(dossier);
Pc=[0.15 0.2 0.3 0.5];

for i=1:4
    
   path=['/home/bauj0001/projets/WIM/output/',char(dossier(i)),'/',char(dossier(i)),'.nc'];
   Dave=ncread(path,'Dave');
   x=ncread(path,'x_axis');
   
   Hs=ncread(path,'significant height');
   subplot(1,2,1)
   plot(x,Dave)
   hold on
   
   subplot(1,2,2)
   Lmiz=length(find(Dave>0 & Dave<396.7623))*(x(2)-x(1));
   plot(Pc(i),Lmiz,'ko','markerFacecolor','k')
   hold on
   
end

figure
 for i=1:size(Hs,1)
     
 hs2(i)=Hs(i,i);
 end
 plot(x,hs2)