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
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	for k=1:ny
13	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