sensitivity_ec.m 827 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=['ec1';'ec2';'ec3';'ec4';'ec5'];
dossier=cellstr(dossier);
Pc=[3e-5 3.5e-5 4e-5 4.5e-5 5e-5 5e-4];

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

   subplot(1,2,1)
   Dave(find(Dave==0 | Dave>=396))=NaN;
   plot(x,Dave)
   hold on
   
   subplot(1,2,2)
   Lmiz=length(find(Dave>0 & Dave<396))*(x(2)-x(1));
   plot(Pc(i),Lmiz,'ko','markerFacecolor','k')
   hold on
   
end

figure