osse/utils/hmaker.m


nx=120;
ny=31;
nz=29;
nf=2;
vectorsize=nx*ny*nz*nf;
dim=1;
fid = fopen('hfile.txt','w')

for i=1:nf
   for j=1:nz
      for k=1:ny 
        for l=1:nx
            if (any(j==[5 10 15 20 25]) && ... % 5 levels
                any(i==[1 2]) && ... % u, v
                mod(k,2)==0)
%          if mod(dim,359)==0
                 fprintf(fid,'%i\n',dim);
            end
            dim=dim+1;
         end
      end
   end
end


1	afe	1.1
2			nx=120;
3			ny=31;
4			nz=29;
5			nf=2;
6			vectorsize=nxnynz*nf;
7			dim=1;
8			fid = fopen('hfile.txt','w')
9
10			for i=1:nf
11			for j=1:nz
12	afe	1.2	for k=1:ny
13	afe	1.1	for l=1:nx
14			if (any(j==[5 10 15 20 25]) && ... % 5 levels
15			any(i==[1 2]) && ... % u, v
16			mod(k,2)==0)
17			% if mod(dim,359)==0
18			fprintf(fid,'%i\n',dim);
19			end
20			dim=dim+1;
21			end
22			end
23			end
24			end
25
26