1 |
|
nx=32;ny=32;o=3;nt=6; |
2 |
|
N=[nx+2*o ny+2*o nt]; |
3 |
|
rtype='real*8'; |
4 |
|
endian='l'; |
5 |
|
|
6 |
|
fid=fopen('fort.10','r',endian); |
7 |
|
% Centered |
8 |
|
xc=f77read(fid,N,rtype); |
9 |
|
yc=f77read(fid,N,rtype); |
10 |
|
dxf=f77read(fid,N,rtype); |
11 |
|
dyf=f77read(fid,N,rtype); |
12 |
|
ra=f77read(fid,N,rtype); |
13 |
|
% Corners |
14 |
|
xg=f77read(fid,N,rtype); |
15 |
|
yg=f77read(fid,N,rtype); |
16 |
|
dxv=f77read(fid,N,rtype); |
17 |
|
dyu=f77read(fid,N,rtype); |
18 |
|
raz=f77read(fid,N,rtype); |
19 |
|
% Staggered |
20 |
|
dxc=f77read(fid,N,rtype); |
21 |
|
dyc=f77read(fid,N,rtype); |
22 |
|
raw=f77read(fid,N,rtype); |
23 |
|
ras=f77read(fid,N,rtype); |
24 |
|
dxg=f77read(fid,N,rtype); |
25 |
|
dyg=f77read(fid,N,rtype); |
26 |
|
aW=f77read(fid,N,rtype); |
27 |
|
aS=f77read(fid,N,rtype); |
28 |
|
pC=f77read(fid,N,rtype); |
29 |
|
pW=f77read(fid,N,rtype); |
30 |
|
pS=f77read(fid,N,rtype); |
31 |
|
fclose(fid); |