1 |
% This is a matlab script that generates the input data |
% This is a matlab script that generates the input data |
2 |
|
|
3 |
|
% $Header$ |
4 |
|
% $Name$ |
5 |
|
|
6 |
% Dimensions of grid |
% Dimensions of grid |
7 |
nx=100; |
nx=100; |
8 |
ny=100; |
ny=100; |
51 |
qs=Q/9; qs( find(r2 > Rc*Rc) )=0.; |
qs=Q/9; qs( find(r2 > Rc*Rc) )=0.; |
52 |
Qc=Qc+qs; |
Qc=Qc+qs; |
53 |
end ; end |
end ; end |
54 |
fid=fopen('Qnet.bin','w','b'); fwrite(fid,Qc,'real*8'); fclose(fid); |
%fid=fopen('Qnet_p64.bin','w','b'); fwrite(fid,Qc,'real*8'); fclose(fid); |
55 |
|
fid=fopen('Qnet_p32.bin','w','b'); fwrite(fid,Qc,'real*4'); fclose(fid); |
56 |
|
|
57 |
var=2*pi*[0:1000]/1000; xl=Rc*cos(var); yl=Rc*sin(var); |
var=2*pi*[0:1000]/1000; xl=Rc*cos(var); yl=Rc*sin(var); |
58 |
figure(1);clf; |
figure(1);clf; |