HDF5_parser/thesis_parser.m at main · cw646/HDF5_parser · GitHub

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
72
73
74
clear all;
start = 1.0;
count = 13800.0;
%read = h5read('301/particles0010.h5','/table',start,count)
store = [];
store1 = [];
shear1 = [];
visco1= [];
strain1= [];
test4 = [];
octa = [];
sq = [];
dx= [];


% for i=1:999
%
% FileName  = sprintf('301/particles%04d.h5', i);
%
% read = h5read(FileName,'/table',start,count);
% store(i) = read.tau_yz;
% shear(i) = read.svars_4;
%
% end


 j = 1;

FileName  = sprintf('results/folder/particles000000.h5')

read = h5read(FileName,'/table',start,count)

max_x = max(read.coord_x);
min_x = min(read.coord_x);
store1(j) = max_x - min_x;
store2(j) = mean(read.svars_1)
shear1(j) = mean(read.svars_6)

j = j+1;

for i = 1000:1000:210000

    i
FileName  = sprintf('results/folder/particles%06d.h5', i);

read = h5read(FileName,'/table',start,count);

max_x = max(read.coord_x);
min_x = min(read.coord_x);
store1(j) = max_x - min_x
store2(j) = mean(read.svars_1)
shear1(j) = mean(read.svars_6)
dx(j) = store1(j) - store1(j-1);
dxmm(j) = dx(j)*1000;
j = j+1;


end

%[zz] = meshgrid(store1)


 figure
 plot(store1)
 figure
 plot(dx)
 %surf(shear1,strain1,zz)
 %scatter(shear1,store1)
 %plot(visco1, store1, 'x-')


%scatter3(shear1,store1, strain5)
%  plot(strain3,'o-')
%  plot(strain4,'x-')