1 |
edhill |
1.1 |
% |
2 |
|
|
% Ed Hill |
3 |
|
|
% |
4 |
|
|
% Quick script to double the lengths/areas around the face edges. |
5 |
|
|
|
6 |
|
|
fin = 'a.face001.bin'; |
7 |
|
|
fout = 'b.face001.bin'; |
8 |
|
|
|
9 |
|
|
sNx = 102; |
10 |
|
|
sNy = 41; |
11 |
|
|
|
12 |
|
|
vnall = { 'XC','YC','dxF','dyF','rA','XG','YG','dxV', ... |
13 |
|
|
'dyU','rAz','dxC','dyC','rAw','rAs','dxG','dyG' }; |
14 |
|
|
|
15 |
|
|
gvars = { 'XC','YC','DXF','DYF','RA','XG','YG','DXV', ... |
16 |
|
|
'DYU','RAZ','DXC','DYC','RAW','RAS','DXG','DYG' }; |
17 |
|
|
|
18 |
|
|
nslab = (sNx + 1)*(sNy + 1); |
19 |
|
|
|
20 |
|
|
grid = {}; |
21 |
|
|
for ivar = 1:length(gvars) |
22 |
|
|
|
23 |
|
|
disp([ ' reading : ' gvars{ivar} ]); |
24 |
|
|
|
25 |
|
|
gid = fopen(fin, 'r', 'ieee-be'); |
26 |
|
|
fseek(gid,(ivar-1)*nslab*8,'bof'); |
27 |
|
|
tmp = reshape(fread(gid,nslab,'real*8',0,'ieee-be'),[sNx+1,sNy+1]); |
28 |
|
|
fclose(gid); |
29 |
|
|
grid.(gvars{ivar}) = tmp; |
30 |
|
|
|
31 |
|
|
end |
32 |
|
|
|
33 |
|
|
vn_fonly = { 'XC' 'YC' 'dxF' 'dyF' 'rA' 'XG' 'YG' 'dxG' 'dyG' }; |
34 |
|
|
vn_fedge = { 'rAw' 'rAs' 'rAz' 'dxV' 'dyU' 'dxC' 'dyC' }; |
35 |
|
|
|
36 |
|
|
gnew = grid; |
37 |
|
|
|
38 |
|
|
% RAW |
39 |
|
|
gnew.RAW(1 ,:) = 2*grid.RAW(1 ,:); % Bot |
40 |
|
|
gnew.RAW(end,:) = 2*grid.RAW(end,:); % Top |
41 |
|
|
subplot(1,2,1), surf(grid.RAW) |
42 |
|
|
title('old'), shading interp, cax = caxis; |
43 |
|
|
subplot(1,2,2), surf(gnew.RAW) |
44 |
|
|
title('new'), shading interp, caxis(cax); |
45 |
|
|
% RAS |
46 |
|
|
gnew.RAS(:,1 ) = 2*grid.RAS(:,1 ); % L |
47 |
|
|
gnew.RAS(:,end) = 2*grid.RAS(:,end); % R |
48 |
|
|
subplot(1,2,1), surf(grid.RAS) |
49 |
|
|
title('old'), shading interp, cax = caxis; |
50 |
|
|
subplot(1,2,2), surf(gnew.RAS) |
51 |
|
|
title('new'), shading interp, caxis(cax); |
52 |
|
|
% RAZ |
53 |
|
|
gnew.RAZ(1 ,:) = 2*gnew.RAZ(1 ,:); % Bot |
54 |
|
|
gnew.RAZ(end,:) = 2*gnew.RAZ(end,:); % Top |
55 |
|
|
gnew.RAZ(:,1 ) = 2*gnew.RAZ(:,1 ); % L |
56 |
|
|
gnew.RAZ(:,end) = 2*gnew.RAZ(:,end); % R |
57 |
|
|
subplot(1,2,1), surf(grid.RAZ) |
58 |
|
|
title('old'), shading interp, cax = caxis; |
59 |
|
|
subplot(1,2,2), surf(gnew.RAZ) |
60 |
|
|
title('new'), shading interp, caxis(cax); |
61 |
|
|
% DXV |
62 |
|
|
gnew.DXV(1 ,:) = 2*grid.DXV(1 ,:); % Bot |
63 |
|
|
gnew.DXV(end,:) = 2*grid.DXV(end,:); % Top |
64 |
|
|
subplot(1,2,1), surf(grid.DXV) |
65 |
|
|
title('old'), shading interp, cax = caxis; |
66 |
|
|
subplot(1,2,2), surf(gnew.DXV) |
67 |
|
|
title('new'), shading interp, caxis(cax); |
68 |
|
|
% DYU |
69 |
|
|
gnew.DYU(:,1 ) = 2*grid.DYU(:,1 ); % L |
70 |
|
|
gnew.DYU(:,end) = 2*grid.DYU(:,end); % R |
71 |
|
|
subplot(1,2,1), surf(grid.DYU) |
72 |
|
|
title('old'), shading interp, cax = caxis; |
73 |
|
|
subplot(1,2,2), surf(gnew.DYU) |
74 |
|
|
title('new'), shading interp, caxis(cax); |
75 |
|
|
% DXC |
76 |
|
|
gnew.DXC(1 ,:) = 2*grid.DXC(1 ,:); % Bot |
77 |
|
|
gnew.DXC(end,:) = 2*grid.DXC(end,:); % Top |
78 |
|
|
subplot(1,2,1), surf(grid.DXC) |
79 |
|
|
title('old'), shading interp, cax = caxis; |
80 |
|
|
subplot(1,2,2), surf(gnew.DXC) |
81 |
|
|
title('new'), shading interp, caxis(cax); |
82 |
|
|
% DYC |
83 |
|
|
gnew.DYC(:,1 ) = 2*grid.DYC(:,1 ); % L |
84 |
|
|
gnew.DYC(:,end) = 2*grid.DYC(:,end); % R |
85 |
|
|
subplot(1,2,1), surf(grid.DYC) |
86 |
|
|
title('old'), shading interp, cax = caxis; |
87 |
|
|
subplot(1,2,2), surf(gnew.DYC) |
88 |
|
|
title('new'), shading interp, caxis(cax); |
89 |
|
|
|
90 |
|
|
fid = fopen(fout, 'w', 'ieee-be'); |
91 |
|
|
for i = 1:length(gvars) |
92 |
|
|
fwrite(fid,gnew.(gvars{i}),'real*8',0,'ieee-be'); |
93 |
|
|
end |
94 |
|
|
fclose(fid); |
95 |
|
|
|