Commit 43a95ea8 authored by Jérémy Baudry's avatar Jérémy Baudry

'clean'

parent 5fb59e80
File deleted
compute-2-8.local
processing 2% processing 3% processing 4% processing 5% processing 6% processing 7% processing 8% processing 9% processing 10% processing 11% processing 12% processing 13% processing 14% processing 16% processing 16% processing 17% processing 18% processing 19% processing 20% processing 21% processing 22% processing 23% processing 24% processing 25% processing 26% processing 28% processing 28% processing 29% processing 31% processing 31% processing 32% processing 33% processing 34% processing 35% processing 36% processing 37% processing 38% processing 39% processing 40% processing 41% processing 42% processing 43% processing 44% processing 45% processing 46% processing 47% processing 48% processing 49% processing 50% processing 51% processing 52% processing 53% processing 55% processing 55% processing 56% processing 57% processing 58% processing 59% processing 61% processing 61% processing 62% processing 63% processing 64% processing 65% processing 66% processing 67% processing 68% processing 69% processing 70% processing 71% processing 72% processing 73% processing 74% processing 75% processing 76% processing 77% processing 78% processing 79% processing 80% processing 81% processing 82% processing 83% processing 84% processing 85% processing 86% processing 87% processing 88% processing 89% processing 90% processing 91% processing 92% processing 93% processing 94% processing 95% processing 96% processing 97% processing 98% processing 99% processing 100%
simulation completed!
compute-2-7.local
processing 2% processing 3% processing 4% processing 5% processing 6% processing 7% processing 8% processing 9% processing 10% processing 11% processing 12% processing 13% processing 14% processing 16% processing 16% processing 17% processing 18% processing 19% processing 20% processing 21% processing 22% processing 23% processing 24% processing 25% processing 26% processing 28% processing 28% processing 29% processing 31% processing 31% processing 32% processing 33% processing 34% processing 35% processing 36% processing 37% processing 38% processing 39% processing 40% processing 41% processing 42% processing 43% processing 44% processing 45% processing 46% processing 47% processing 48% processing 49% processing 50% processing 51% processing 52% processing 53% processing 55% processing 55% processing 56% processing 57% processing 58% processing 59% processing 61% processing 61% processing 62% processing 63% processing 64% processing 65% processing 66% processing 67% processing 68% processing 69% processing 70% processing 71% processing 72% processing 73% processing 74% processing 75% processing 76% processing 77% processing 78% processing 79% processing 80% processing 81% processing 82% processing 83% processing 84% processing 85% processing 86% processing 87% processing 88% processing 89% processing 90% processing 91% processing 92% processing 93% processing 94% processing 95% processing 96% processing 97% processing 98% processing 99% processing 100%
simulation completed!
compute-2-8.local
processing 2% processing 3% processing 4% processing 5% processing 6% processing 7% processing 8% processing 9% processing 10% processing 11% processing 12% processing 13% processing 14% processing 16% processing 16% processing 17% processing 18% processing 19% processing 20% processing 21% processing 22% processing 23% processing 24% processing 25% processing 26% processing 28% processing 28% processing 29% processing 31% processing 31% processing 32% processing 33% processing 34% processing 35% processing 36% processing 37% processing 38% processing 39% processing 40% processing 41% processing 42% processing 43% processing 44% processing 45% processing 46% processing 47% processing 48% processing 49% processing 50% processing 51% processing 52% processing 53% processing 55% processing 55% processing 56% processing 57% processing 58% processing 59% processing 61% processing 61% processing 62% processing 63% processing 64% processing 65% processing 66% processing 67% processing 68% processing 69% processing 70% processing 71% processing 72% processing 73% processing 74% processing 75% processing 76% processing 77% processing 78% processing 79% processing 80% processing 81% processing 82% processing 83% processing 84% processing 85% processing 86% processing 87% processing 88% processing 89% processing 90% processing 91% processing 92% processing 93% processing 94% processing 95% processing 96% processing 97% processing 98% processing 99% processing 100%
simulation completed!
compute-2-7.local
processing 2% processing 3% processing 4% processing 5% processing 6% processing 7% processing 8% processing 9% processing 10% processing 11% processing 12% processing 13% processing 14% processing 16% processing 16% processing 17% processing 18% processing 19% processing 20% processing 21% processing 22% processing 23% processing 24% processing 25% processing 26% processing 28% processing 28% processing 29% processing 31% processing 31% processing 32% processing 33% processing 34% processing 35% processing 36% processing 37% processing 38% processing 39% processing 40% processing 41% processing 42% processing 43% processing 44% processing 45% processing 46% processing 47% processing 48% processing 49% processing 50% processing 51% processing 52% processing 53% processing 55% processing 55% processing 56% processing 57% processing 58% processing 59% processing 61% processing 61% processing 62% processing 63% processing 64% processing 65% processing 66% processing 67% processing 68% processing 69% processing 70% processing 71% processing 72% processing 73% processing 74% processing 75% processing 76% processing 77% processing 78% processing 79% processing 80% processing 81% processing 82% processing 83% processing 84% processing 85% processing 86% processing 87% processing 88% processing 89% processing 90% processing 91% processing 92% processing 93% processing 94% processing 95% processing 96% processing 97% processing 98% processing 99% processing 100%
simulation completed!
compute-2-8.local
processing 2% processing 3% processing 4% processing 5% processing 6% processing 7% processing 8% processing 9% processing 10% processing 11% processing 12% processing 13% processing 14% processing 16% processing 16% processing 17% processing 18% processing 19% processing 20% processing 21% processing 22% processing 23% processing 24% processing 25% processing 26% processing 28% processing 28% processing 29% processing 31% processing 31% processing 32% processing 33% processing 34% processing 35% processing 36% processing 37% processing 38% processing 39% processing 40% processing 41% processing 42% processing 43% processing 44% processing 45% processing 46% processing 47% processing 48% processing 49% processing 50% processing 51% processing 52% processing 53% processing 55% processing 55% processing 56% processing 57% processing 58% processing 59% processing 61% processing 61% processing 62% processing 63% processing 64% processing 65% processing 66% processing 67% processing 68% processing 69% processing 70% processing 71% processing 72% processing 73% processing 74% processing 75% processing 76% processing 77% processing 78% processing 79% processing 80% processing 81% processing 82% processing 83% processing 84% processing 85% processing 86% processing 87% processing 88% processing 89% processing 90% processing 91% processing 92% processing 93% processing 94% processing 95% processing 96% processing 97% processing 98% processing 99% processing 100%
simulation completed!
compute-2-7.local
processing 2% processing 3% processing 4% processing 5% processing 6% processing 7% processing 8% processing 9% processing 10% processing 11% processing 12% processing 13% processing 14% processing 16% processing 16% processing 17% processing 18% processing 19% processing 20% processing 21% processing 22% processing 23% processing 24% processing 25% processing 26% processing 28% processing 28% processing 29% processing 31% processing 31% processing 32% processing 33% processing 34% processing 35% processing 36% processing 37% processing 38% processing 39% processing 40% processing 41% processing 42% processing 43% processing 44% processing 45% processing 46% processing 47% processing 48% processing 49% processing 50% processing 51% processing 52% processing 53% processing 55% processing 55% processing 56% processing 57% processing 58% processing 59% processing 61% processing 61% processing 62% processing 63% processing 64% processing 65% processing 66% processing 67% processing 68% processing 69% processing 70% processing 71% processing 72% processing 73% processing 74% processing 75% processing 76% processing 77% processing 78% processing 79% processing 80% processing 81% processing 82% processing 83% processing 84% processing 85% processing 86% processing 87% processing 88% processing 89% processing 90% processing 91% processing 92% processing 93% processing 94% processing 95% processing 96% processing 97% processing 98% processing 99% processing 100%
simulation completed!
compute-2-8.local
processing 2% processing 3% processing 4% processing 5% processing 6% processing 7% processing 8% processing 9% processing 10% processing 11% processing 12% processing 13% processing 14% processing 16% processing 16% processing 17% processing 18% processing 19% processing 20% processing 21% processing 22% processing 23% processing 24% processing 25% processing 26% processing 28% processing 28% processing 29% processing 31% processing 31% processing 32% processing 33% processing 34% processing 35% processing 36% processing 37% processing 38% processing 39% processing 40% processing 41% processing 42% processing 43% processing 44% processing 45% processing 46% processing 47% processing 48% processing 49% processing 50% processing 51% processing 52% processing 53% processing 55% processing 55% processing 56% processing 57% processing 58% processing 59% processing 61% processing 61% processing 62% processing 63% processing 64% processing 65% processing 66% processing 67% processing 68% processing 69% processing 70% processing 71% processing 72% processing 73% processing 74% processing 75% processing 76% processing 77% processing 78% processing 79% processing 80% processing 81% processing 82% processing 83% processing 84% processing 85% processing 86% processing 87% processing 88% processing 89% processing 90% processing 91% processing 92% processing 93% processing 94% processing 95% processing 96% processing 97% processing 98% processing 99% processing 100%
simulation completed!
compute-2-7.local
processing 2% processing 3% processing 4% processing 5% processing 6% processing 7% processing 8% processing 9% processing 10% processing 11% processing 12% processing 13% processing 14% processing 16% processing 16% processing 17% processing 18% processing 19% processing 20% processing 21% processing 22% processing 23% processing 24% processing 25% processing 26% processing 28% processing 28% processing 29% processing 31% processing 31% processing 32% processing 33% processing 34% processing 35% processing 36% processing 37% processing 38% processing 39% processing 40% processing 41% processing 42% processing 43% processing 44% processing 45% processing 46% processing 47% processing 48% processing 49% processing 50% processing 51% processing 52% processing 53% processing 55% processing 55% processing 56% processing 57% processing 58% processing 59% processing 61% processing 61% processing 62% processing 63% processing 64% processing 65% processing 66% processing 67% processing 68% processing 69% processing 70% processing 71% processing 72% processing 73% processing 74% processing 75% processing 76% processing 77% processing 78% processing 79% processing 80% processing 81% processing 82% processing 83% processing 84% processing 85% processing 86% processing 87% processing 88% processing 89% processing 90% processing 91% processing 92% processing 93% processing 94% processing 95% processing 96% processing 97% processing 98% processing 99% processing 100%
simulation completed!
compute-2-8.local
processing 2% processing 3% processing 4% processing 5% processing 6% processing 7% processing 8% processing 9% processing 10% processing 11% processing 12% processing 13% processing 14% processing 16% processing 16% processing 17% processing 18% processing 19% processing 20% processing 21% processing 22% processing 23% processing 24% processing 25% processing 26% processing 28% processing 28% processing 29% processing 31% processing 31% processing 32% processing 33% processing 34% processing 35% processing 36% processing 37% processing 38% processing 39% processing 40% processing 41% processing 42% processing 43% processing 44% processing 45% processing 46% processing 47% processing 48% processing 49% processing 50% processing 51% processing 52% processing 53% processing 55% processing 55% processing 56% processing 57% processing 58% processing 59% processing 61% processing 61% processing 62% processing 63% processing 64% processing 65% processing 66% processing 67% processing 68% processing 69% processing 70% processing 71% processing 72% processing 73% processing 74% processing 75% processing 76% processing 77% processing 78% processing 79% processing 80% processing 81% processing 82% processing 83% processing 84% processing 85% processing 86% processing 87% processing 88% processing 89% processing 90% processing 91% processing 92% processing 93% processing 94% processing 95% processing 96% processing 97% processing 98% processing 99% processing 100%
simulation completed!
compute-2-7.local
processing 2% processing 3% processing 4% processing 5% processing 6% processing 7% processing 8% processing 9% processing 10% processing 11% processing 12% processing 13% processing 14% processing 16% processing 16% processing 17% processing 18% processing 19% processing 20% processing 21% processing 22% processing 23% processing 24% processing 25% processing 26% processing 28% processing 28% processing 29% processing 31% processing 31% processing 32% processing 33% processing 34% processing 35% processing 36% processing 37% processing 38% processing 39% processing 40% processing 41% processing 42% processing 43% processing 44% processing 45% processing 46% processing 47% processing 48% processing 49% processing 50% processing 51% processing 52% processing 53% processing 55% processing 55% processing 56% processing 57% processing 58% processing 59% processing 61% processing 61% processing 62% processing 63% processing 64% processing 65% processing 66% processing 67% processing 68% processing 69% processing 70% processing 71% processing 72% processing 73% processing 74% processing 75% processing 76% processing 77% processing 78% processing 79% processing 80% processing 81% processing 82% processing 83% processing 84% processing 85% processing 86% processing 87% processing 88% processing 89% processing 90% processing 91% processing 92% processing 93% processing 94% processing 95% processing 96% processing 97% processing 98% processing 99% processing 100%
simulation completed!
&waves_parameters
Tm =6
Hs =1
disp =0
/
&model_parameter
nbin =10
dx =5000
Cfl =0.1
/
&spectrum_parameters
nfreq =60
Tmin =2.5
Tmax =20
alpha_s =0.2044
beta_s =1.2500
gamma_s =3.3
/
&ice_parameters
cice =0.5
hice =1
D0 =500
gam =1.5
Dmin =20
/
clear all;
close all;
dossier=['temp1'];
path=['/home/bauj0001/projets/WIM/output/',dossier,'/',dossier,'.nc'];
x=ncread(path,'x_axis');
t=ncread(path,'time');
om=ncread(path,'omega');
spectre=ncread(path,'Spectrum');
Dave=ncread(path,'Dave');
Dmax=ncread(path,'Dmax');
FSD=ncread(path,'Floe size distribution');
Fsize=ncread(path,'floe size');
thick=ncread(path,'Ice thickness');
conc=ncread(path,'Ice concentration');
Hs=ncread(path,'significant height');
IDT=ncread(path,'Ice Thickness Distribution');
hcat=ncread(path,'thickness categories');
fig=figure(1);
u = fig.Units;
fig.Units = 'inches';
set(fig,'outerposition',[3 3 12 10])
cm=colormap('jet');
cm1(1,:)=[1 1 1];
for i=2:10
cm1(i,:)=[1-0.08*i 1-0.06*i 1];
end
cmap(1:10,:)=cm1;
cmap(11:59,:)=cm(16:end,:);
w = waitforbuttonpress;
%figure
%filename = 'joint_distribution_anim.gif';
for j=1:length(t)
FSD1=reshape(FSD(j,j,:,:),length(Fsize),length(hcat));
IDT1=IDT(j,:);
for i=1:length(hcat)
FSD1(:,i)=FSD1(:,i)*IDT1(i);
nfloe(:,i)=((FSD1(:,i)*(5000)^2)./Fsize.^2) ;
end
nfloe=nfloe/sum(sum(nfloe));
FSD2(j,:)=sum(FSD1,2);
subplot(2,2,1)
hh=surf(hcat,Fsize,nfloe);
shading interp
%
zoom(1.5)
caxis([0 0.06])
colormap(cmap)
axis([0.1 10 0 400 0 0.1])
view(108,60)
xlabel('Ice thickness [m]')
ylabel('Floe size [m]')
zlabel('area fraction')
subplot(2,2,2)
nfloe2=((FSD2(j,:)*5000^2)./Fsize.^2')/sum((FSD2(j,:)*5000^2)./Fsize.^2');
plot((Fsize),nfloe2)
axis([0 400 0 1])
xlabel('Floe size [m]')
ylabel('normalized number of floe')
subplot(2,2,3:4)
plot(x,thick)
hold on
plot(x(j),thick(j),'ro','markerfacecolor','r')
hold off
xlabel('distance [km]')
ylabel('Average ice thickness [m]')
%pause(0.02)
w = waitforbuttonpress;
% drawnow
% frame = getframe(1);
% im = frame2im(frame);
% [A,map] = rgb2ind(im,256);
% if j == 1;
% imwrite(A,map,filename,'gif','LoopCount',Inf,'DelayTime',0);
% else
% imwrite(A,map,filename,'gif','WriteMode','append','DelayTime',0);
% end
end
clear all;
close all;
for i=1:100
dossier=['pc_ec',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));
i
end
Lmiz2=reshape(Lmiz,10,10);
Ec=[3e-5 4e-5 5e-5 6e-5 7e-5 8e-5 9e-5 1e-4 2e-4 3e-4];
legende=['$P_c=0.10$';'$P_c=0.15$';'$P_c=0.20$';'$P_c=0.25$';'$P_c=0.30$'...
;'$P_c=0.35$';'$P_c=0.40$';'$P_c=0.45$';'$P_c=0.50$';'$P_c=0.55$'];
legende=cellstr(legende);
cmap=jet(10);
M=['*','o','+','*','o','+','*','o','+','*'];
fig=figure;
fig.Units='inches';
fig.Position=[7.0729 6.7396 4.2708 2.600];
fig.OuterPosition=[7.0729 6.7396 4.2708 3.600];
fig.PaperUnits = 'inches';
fig.PaperPosition = [0 0 4.2708 2.600];
fig.PaperPositionMode = 'manual';
fig.PaperOrientation='landscape';
fig.PaperSize=[4.2708 2.600];
for i=1:10
plot(Ec,Lmiz2(i,:),'color',cmap(i,:),'Marker',M(i));
hold on
end
xlabel('$\varepsilon_c$','interpreter','latex')
ylabel('$L_{MIZ}$ [km]','interpreter','latex')
a=gca;
% a.Position(2)=a.Position(2)+0.2;
% a.Position(3)=a.Position(3)-0.2;
% a.Position(4)=a.Position(4)-0.2;
l=legend(char(legende));
set(l,'interpreter','latex','Box','off','location','eastoutside');
%------------------------------------------------------------------------------------------------------
figure
cmap=jet(10);
for j=1:10
dossier=['ec',num2str(j)];
path=['/home/bauj0001/projets/WIM/output/',dossier,'/',dossier,'.nc'];
x=ncread(path,'x_axis');
t=ncread(path,'time');
om=ncread(path,'omega');
spectre=ncread(path,'Spectrum');
Dave=ncread(path,'Dave');
Dmax=ncread(path,'Dmax');
FSD=ncread(path,'Floe size distribution');
Fsize=ncread(path,'floe size');
thick=ncread(path,'Ice thickness');
conc=ncread(path,'Ice concentration');
Hs=ncread(path,'significant height');
IDT=ncread(path,'Ice Thickness Distribution');
hcat=ncread(path,'thickness categories');
FSD1=reshape(FSD(10,10,:,:),length(Fsize),length(hcat));
IDT1=IDT(10,:);
for i=1:length(hcat)
FSD1(:,i)=FSD1(:,i)*IDT1(i);
nfloe(:,i)=((FSD1(:,i)*(5000)^2)./Fsize.^2) ;
end
nfloe=nfloe/sum(sum(nfloe));
FSD2=sum(FSD1,2);
nfloe2=((FSD2*5000^2)./Fsize.^2)/sum((FSD2*5000^2)./Fsize.^2);
plot((Fsize),nfloe2,'color',cmap(j,:))
axis([0 400 0 1])
xlabel('Floe size [m]')
ylabel('normalized number of floe')
hold on
axis([0 250 0 0.2])
a=a+10
end
legende=['$\varepsilon_c=3\times 10^{-5}$';'$\varepsilon_c=4\times 10^{-5}$';'$\varepsilon_c=5\times 10^{-5}$'...
;'$\varepsilon_c=6\times 10^{-5}$';'$\varepsilon_c=7\times 10^{-5}$';'$\varepsilon_c=8\times 10^{-5}$'...
;'$\varepsilon_c=9\times 10^{-5}$';'$\varepsilon_c=1\times 10^{-4}$';'$\varepsilon_c=2\times 10^{-4}$'...
;'$\varepsilon_c=3\times 10^{-4}$'];
legende=cellstr(legende);
l=legend(char(legende));
set(l,'interpreter','latex','Box','off','location','eastoutside');
%
\ No newline at end of file
clear all;
close all;
dossier=['Pc9'];
path=['/home/bauj0001/projets/WIM/output/',dossier,'/',dossier,'.nc'];
x=ncread(path,'x_axis');
t=ncread(path,'time');
om=ncread(path,'omega');
spectre=ncread(path,'Spectrum');
Dave=ncread(path,'Dave');
Dmax=ncread(path,'Dmax');
FSD=ncread(path,'Floe size distribution');
Fsize=ncread(path,'floe size');
thick=ncread(path,'Ice thickness');
conc=ncread(path,'Ice concentration');
Hs=ncread(path,'significant height');
IDT=ncread(path,'Ice Thickness Distribution');
hcat=ncread(path,'thickness categories');
FSD1=reshape(FSD(:,10,:,:),size(FSD,1),size(FSD,3),size(FSD,4));
IDT1=IDT(25,:)';
for i=1:length(hcat)
FSD1(:,:,i)=FSD1(:,:,i).*IDT1(i);
end
a=1;
for i=1:length(t)
FSD2(:,i)=sum(reshape(FSD1(i,:,:),size(FSD1,2),size(FSD1,3)),2);
FSD2(find(FSD2(:,i)<0.0001),i)=0;
if FSD2(end,i)<0.999
FSD3(:,a)=FSD2(:,i);
a=a+1;
end
end
a=1;
FSD4(:,a)=FSD3(:,1);
for i=2:size(FSD3,2)
if max(FSD3(:,i))==max(FSD3(:,i-1))
else
a=a+1;
FSD4(:,a)=FSD3(:,i);
end
end
cmap=jet(size(FSD4,2));
for i=1:size(FSD4,2)
for j=1:size(FSD4,1)
CFSD(j,i)=sum(FSD4(j:end,i))
end
subplot(1,2,1)
plot(Fsize,FSD4(:,i),'color',cmap(i,:))
hold on
xlabel('Floe size [m]','interpreter','latex')
ylabel('Area fraction','interpreter','latex')
subplot(1,2,2)
plot(Fsize,CFSD(:,i),'color',cmap(i,:))
hold on
xlabel('Floe size [m]','interpreter','latex')
ylabel(' cumulative area fraction','interpreter','latex')
end
c=colorbar;
colormap(jet)
ylabel(c,'Time [min]')
This diff is collapsed.
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