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); |