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
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
|
Path='.\'
Case = 'zoom'
Angle = 'test'
[A.ind, A.Time]= textread(['Pdt_simu_20min.dat'],'%f %10s');
[P.loc]= textread(['loc_Profils_2m.dat'],'%f');
%% 900 = nombre de pas de temps
%%% s = nombre de profils
%%
for i=1:length (A.Time)
%% Load des pas de temps à ploter
load([Path, 'RES_surf_', Case, '.mat'], ['Vksi_' char(char(A.Time(i)))]);
load([Path, 'RES_surf_', Case, '.mat'], ['Veta_' char(char(A.Time(i)))]);
load([Path, 'RES_surf_', Case, '.mat'], ['Watlev_' char(char(A.Time(i)))]);
eval(['Vksi_pdt=Vksi_' char(char(A.Time(i)))]);
eval(['Veta_pdt=Veta_' char(char(A.Time(i)))]);
eval(['Watlev_pdt=Watlev_' char(char(A.Time(i)))]);
Cur_pdt = sqrt(Vksi_pdt.^2 + Veta_pdt.^2);
Mask_pdt = Cur_pdt*0;
Watlev_Mask = Watlev_pdt + Mask_pdt;
%% Figures 2D de watlev et courants
%%% ici en OFF
fig_watlev=figure;
set(fig_watlev,'visible','off');
pcolor(double(Watlev_Mask)); shading flat
set(gca, 'DataAspectRatio', [1 1 1])
colorbar
caxis([-2 4])
%% Série temporelle extraite sur profils line1D
% Localisation des profils en config ZOOM
for p=1:length (P.loc)
% Watlev_prof1 = Watlev_Mask((P.loc(p)),:);
Watlev_prof1 = Watlev_Mask(:,(P.loc(p)));
% figure(1);
% plot(Watlev_Mask((p),:),'-b'),hold on;
isnan(Watlev_prof1);
IndWatlevNull1 = find(isnan(Watlev_prof1) == 1);
Runup_prof1(p) = Watlev_prof1(IndWatlevNull1(end)+1);
ResProf_pdt(i,p) = Runup_prof1(p);
end
%%% clear all
clear Vksi_pdt
clear Veta_pdt
clear Watlev_pdt
clear Mask_pdt
clear Cur_pdt
clear Watlev_Mask
clear (['Vksi_' char(char(A.Time(i)))])
clear (['Veta_' char(char(A.Time(i)))])
clear (['Watlev_' char(char(A.Time(i)))])
end |
