-
Notifications
You must be signed in to change notification settings - Fork 0
/
plotContour.m
73 lines (70 loc) · 2.8 KB
/
plotContour.m
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
%
% [X,Z] = meshgrid([0.019:0.001:0.171],[0.019:0.001:0.171]);
% name = {'U', 'V', 'W1', 'W2', 'TKE', 'ReStrUV', 'ReStrUW', 'ReStrVW'};
% field = [4,5,6,7,11,12,13,14];
%
% figure1 = figure;
%
% annotation(figure1,'textbox',...
% [0.385119047619048 0.49379932356257 0.25 0.0375263441236606],...
% 'String',{'Downstream direction is towards the right.'},...
% 'Interpreter','latex',...
% 'HorizontalAlignment','center',...
% 'FontSize',14,...
% 'FitBoxToText','off');
%
% for i = 1 : 8
% Data=griddata(Loc_GenStat_ReStr_good(:,1), Loc_GenStat_ReStr_good(:,3), Loc_GenStat_ReStr_good(:,field(i)), X,Z,'cubic');
% subplot(2,4,i);
% me = mesh(X,Z,Data);
% hold on;
% dumfac = max(max(Data));
% dummy = ones(1,54).*dumfac+1;
% scatter3(Loc_GenStat_ReStr_good(:,1),Loc_GenStat_ReStr_good(:,3),dummy', '.k');
% axis tight;
% view(0,90);
% hold on;
% colorbar;
% title(sprintf('%s', name{i}));
% %axes([0.019 0.152 0.019 0.1714]);
% annotation(figure1,'rectangle',...
% [0.12987012987013 0.580357142857143 0.00285528756957348 0.0431547619047608],...
% 'FaceColor',[0 0 0]);
% annotation(figure1,'rectangle',...
% [0.243042671614102 0.580357142857142 0.00285528756957348 0.0431547619047608],...
% 'FaceColor',[0 0 0]);
% end
[X,Z] = meshgrid([0:0.001:0.342],[0.011:0.001:0.171]);
name = {'U', 'V', 'W1', 'W2', 'TKE', 'ReStrUV', 'ReStrUW', 'ReStrVW'};
field = [4,5,6,7,11,12,13,14];
for i = 1 : 8
figure(i)
Data=griddata(Loc_GenStats_ReStr(:,1), Loc_GenStats_ReStr(:,3), Loc_GenStats_ReStr(:,field(i)), X,Z,'cubic');
me = surf(X,Z,Data);
hold on;
dumfac = max(max(Data));
dummy = ones(1,181).*dumfac+1;
scatter3(Loc_GenStats_ReStr(:,1),Loc_GenStats_ReStr(:,3),dummy', '.k');
axis tight;
view(0,90);
hold on;
colorbar;
me.EdgeColor = 'none';
title(sprintf('%s', name{i}));
annotation(figure(i),'rectangle',...
[0.155357142857143 0.111904761904762 0.0724285714285713 0.228571428571429]);
annotation(figure(i),'rectangle',...
[0.459928571428571 0.111904761904762 0.0668571428571429 0.226190476190477]);
annotation(figure(i),'rectangle',...
[0.767071428571428 0.111904761904762 0.0311428571428575 0.226190476190477]);
annotation(figure(i),'rectangle',...
[0.788499999999999 0.114285714285715 0.00792857142857206 0.133333333333335],...
'FaceColor',[0 0 0]);
annotation(figure(i),'rectangle',...
[0.486714285714285 0.10952380952381 0.0150714285714287 0.133333333333335],...
'FaceColor',[0 0 0]);
annotation(figure(i),'rectangle',...
[0.181357142857143 0.111904761904762 0.0150714285714286 0.133333333333334],...
'FaceColor',[0 0 0]);
saveas(figure(i), sprintf('%s.bmp', name{i}));
end