Bonjour,

Je traite des données sorties d'un modèle hydrodynamique. Le problème ait que je reprends le script de mon prédécesseur en n'ayant pas exactement les mêmes sorties que lui...

Pour moi il traite à partir d'une structure, des matrices. Voici son script :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
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
Je ne peux pas appliquer son script car il a des données sous forme de structure et moi sous forme de matrices. J'aurais donc aimé, le modifier un peu pour reproduire la même chose sachant que moi j'ai déjà toutes les matrices avec le pas de temps mais que je n'ai pas ce pas de temps A.Time(i) pour le mettre dans les calculs après..

en fait je pense qu'il me faut savoir comment "créer" ce pas de temps A.Time qui est sous forme 00000_000 - 00009_999 etc