Commit 963a18cb authored by Jérémy Baudry's avatar Jérémy Baudry
Browse files

version problems

parent fe7bdc5e
...@@ -36,6 +36,9 @@ plot1=subplot(2,2,3:4); ...@@ -36,6 +36,9 @@ plot1=subplot(2,2,3:4);
set(hAx,'color','none') set(hAx,'color','none')
set(hAx(1),'ycolor','b') set(hAx(1),'ycolor','b')
set(hAx(2),'ycolor','g') set(hAx(2),'ycolor','g')
hAx(1).XAxis.Color = 'w';
set(hLine1,'linewidth',2,'color','b') set(hLine1,'linewidth',2,'color','b')
set(hLine2,'linewidth',2,'color','g') set(hLine2,'linewidth',2,'color','g')
h2=gca; h2=gca;
...@@ -83,6 +86,10 @@ set(curs,'UpdateFcn', @plotother) ...@@ -83,6 +86,10 @@ set(curs,'UpdateFcn', @plotother)
evalin('base', 'set(linex,''visible'',''off'')') evalin('base', 'set(linex,''visible'',''off'')')
end end
get(h1,'ylim')
linex=plot(h1,[pos(1) pos(1)],get(h1,'ylim'),'linewidth',1.5,'color',[0.2 0.2 0.6],'linestyle','--');
assignin('base','linex',linex)
FSD1=reshape(FSD(posx,posx,:,:),length(Fsize),length(hcat)); FSD1=reshape(FSD(posx,posx,:,:),length(Fsize),length(hcat));
IDT1=IDT(posx,:); IDT1=IDT(posx,:);
for i=1:length(hcat) for i=1:length(hcat)
...@@ -93,29 +100,32 @@ set(curs,'UpdateFcn', @plotother) ...@@ -93,29 +100,32 @@ set(curs,'UpdateFcn', @plotother)
subplot(2,2,1) subplot(2,2,1)
title('FSTD')
colormap(cmap) colormap(cmap)
pcolor(hcat,Fsize,nfloe); pcolor(hcat,Fsize,nfloe);
shading interp shading flat
c=colorbar; c=colorbar;
set(c,'color','w')
xlabel('Floe thickness [m]','interpreter','latex') xlabel('Floe thickness [m]','interpreter','latex')
ylabel('Floe size','interpreter','latex') ylabel('Floe size','interpreter','latex')
ylabel(c,'normalized number of floes','interpreter','latex') ylabel(c,'normalized number of floes','interpreter','latex')
ax = gca;
ax.XAxis.Color = 'w';
ax.YAxis.Color = 'w';
subplot(2,2,2)
subplot(2,2,2)
title('FSD')
FSD2(posx,:)=sum(FSD1,2);
nfloe2=((sum(FSD1,2)*dx^2)./Fsize.^2)/sum((sum(FSD1,2)*dx^2)./Fsize.^2); nfloe2=((sum(FSD1,2)*dx^2)./Fsize.^2)/sum((sum(FSD1,2)*dx^2)./Fsize.^2);
plot(Fsize,nfloe2,'linewidth',2,'color','m') plot(Fsize,nfloe2,'linewidth',2,'color','m')
set(gca,'color','none') set(gca,'color','none')
ax = gca;
ax.XAxis.Color = 'w';
ax.YAxis.Color = 'w';
xlabel('Floe size [m]','interpreter','latex') xlabel('Floe size [m]','interpreter','latex')
ylabel('normalized number of floes','interpreter','latex') ylabel('normalized number of floes','interpreter','latex')
grid on grid on
get(h1,'ylim')
linex=plot(h1,[pos(1) pos(1)],get(h1,'ylim'),'linewidth',1.5,'color',[0.2 0.2 0.6],'linestyle','--');
assignin('base','linex',linex)
end end
...@@ -125,7 +135,7 @@ end ...@@ -125,7 +135,7 @@ end
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment