1 |
molod |
1.1 |
% load COS & SIN of rotation angle: |
2 |
|
|
fid=fopen('proj_cs32_2uEvN.bin','r','b'); |
3 |
|
|
nx = 192; |
4 |
|
|
ny = 32; |
5 |
|
|
uvEN=fread(fid,nx*ny*2,'real*8'); fclose(fid); |
6 |
|
|
uvEN=reshape(uvEN,[prod(size(uvEN))/2 2]); |
7 |
|
|
AngleCS=uvEN(:,1); |
8 |
|
|
AngleSN=uvEN(:,2); |
9 |
|
|
CS = reshape(AngleCS,[nx ny]); |
10 |
|
|
SN = reshape(AngleSN,[nx ny]); |
11 |
|
|
for fi = 1:6 |
12 |
|
|
is = (fi - 1)*32; |
13 |
|
|
nTUu(:,:,fi) = CS(is+[1:32],1:32); |
14 |
|
|
nTUv(:,:,fi) = -SN(is+[1:32],1:32); |
15 |
|
|
nTVu(:,:,fi) = SN(is+[1:32],1:32); |
16 |
|
|
nTVv(:,:,fi) = CS(is+[1:32],1:32); |
17 |
|
|
end |
18 |
|
|
|
19 |
|
|
load ../set_001/TUV.mat |
20 |
|
|
|
21 |
|
|
oTUu = TUu; |
22 |
|
|
oTUv = TUv; |
23 |
|
|
oTVu = TVu; |
24 |
|
|
oTVv = TVv; |
25 |
|
|
|
26 |
|
|
clear TUu TUv TVu TVv |
27 |
|
|
|
28 |
|
|
figure(1) |
29 |
|
|
subplot(1,1,1) |
30 |
|
|
for fi = 1:6 |
31 |
|
|
|
32 |
|
|
iif = (fi-1)*4; |
33 |
|
|
subplot(6,4,iif+1), surf(oTUu(:,:,fi)), view(2), shading flat |
34 |
|
|
subplot(6,4,iif+2), surf(oTUv(:,:,fi)), view(2), shading flat |
35 |
|
|
subplot(6,4,iif+3), surf(oTVu(:,:,fi)), view(2), shading flat |
36 |
|
|
subplot(6,4,iif+4), surf(oTVv(:,:,fi)), view(2), shading flat |
37 |
|
|
|
38 |
|
|
end |
39 |
|
|
|
40 |
|
|
figure(2) |
41 |
|
|
subplot(1,1,1) |
42 |
|
|
for fi = 1:6 |
43 |
|
|
|
44 |
|
|
iif = (fi-1)*4; |
45 |
|
|
subplot(6,4,iif+1), surf(nTUu(:,:,fi)), view(2), shading flat |
46 |
|
|
subplot(6,4,iif+2), surf(nTUv(:,:,fi)), view(2), shading flat |
47 |
|
|
subplot(6,4,iif+3), surf(nTVu(:,:,fi)), view(2), shading flat |
48 |
|
|
subplot(6,4,iif+4), surf(nTVv(:,:,fi)), view(2), shading flat |
49 |
|
|
|
50 |
|
|
end |
51 |
|
|
|
52 |
|
|
TUu = nTUu; |
53 |
|
|
TUv = nTUv; |
54 |
|
|
TVu = nTVu; |
55 |
|
|
TVv = nTVv; |
56 |
|
|
|
57 |
|
|
save TUV_from_proj_cs32_2uEvN TUu TUv TVu TVv |
58 |
|
|
|
59 |
|
|
|
60 |
|
|
|