 Mesh.cellattr=aa(:);
0020 Mesh.cell=zeros(numel(Mod.M)*size(uninode,1),4);
0021 <span class="keyword">for</span> k=1:nz-1,
0022     <span class="keyword">for</span> j=1:ny-1,
0023         <span class="keyword">for</span> i=1:nx-1,
0024             cellxy=[i i+1 i+nx+1 i+nx ]+ny*(j-1);
0025             cellp=[nxy*(k-1)+cellxy nxy*k+cellxy];
0026             nc=i-1+(nx-1)*(j-1)+(nx-1)*(ny-1)*(k-1);
0027             Mesh.cell(nc*nt+(1:nt),:)=cellp(uninode);
0028         <span class="keyword">end</span>
0029  