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 % get every other ring
        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