-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathpost_process.m
42 lines (30 loc) · 919 Bytes
/
post_process.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
close all
times = [1 1000 2000 4000];
md = '0_1'
speed_history = sqrt(ux_history.^2 + uy_history.^2);
for t = times
% density rho
imagesc(reshape(rho_history(t,2:65,2:129),[64,128]))
colormap parula
axis image
colorbar horizontal
saveas(gcf, strcat('rho-md',md,'-t',int2str(t),'.eps'), 'epsc')
% ux
imagesc(reshape(ux_history(t,2:65,2:129),[64,128]))
colormap parula
axis image
colorbar horizontal
saveas(gcf, strcat('ux-md',md,'-t',int2str(t),'.eps'), 'epsc')
% uy
imagesc(reshape(uy_history(t,2:65,2:129),[64,128]))
colormap parula
axis image
colorbar horizontal
saveas(gcf, strcat('uy-md',md,'-t',int2str(t),'.eps'), 'epsc')
% speed
imagesc(reshape(speed_history(t,2:65,2:129),[64,128]))
colormap parula
axis image
colorbar horizontal
saveas(gcf, strcat('speed-md',md,'-t',int2str(t),'.eps'), 'epsc')
end