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
| d = dir('*.vtu');
for n = 1:numel(d)
V = xml_parseany(fileread(d(n).name));
Position = str2num(V.UnstructuredGrid{1}.Piece{1}.Points{1}.DataArray{1}.CONTENT);
connectivity = str2num(V.UnstructuredGrid{1}.Piece{1}.Cells{1}.DataArray{1}.CONTENT)+1;
U = str2num(V.UnstructuredGrid{1}.Piece{1}.PointData{1}.DataArray{1}.CONTENT);
[vertices{n},a,b] = unique(Position,'rows');
faces{n} = b(connectivity);
fvcdata{n} = U(a);
end
r = [min(min([fvcdata{:}])) max(max([fvcdata{:}]))];
figure
p = patch('vertices',[nan nan],'faces',nan, ...
'facevertexcdata',nan,'cdatamapping','scaled', ...
'facecolor','interp');
set(gca,'clim',r)
colorbar
for n=1:numel(vertices)
set(p,'vertices',vertices{n},'faces',faces{n}, ...
'facevertexcdata',fvcdata{n});
drawnow
title(sprintf('Pas : %d',n))
pause(1)
end |
