taining node,cell and cellattr</span>
0007 
0008 Mesh=[];Mesh.dim=3;
0009 [X,Y,Z]=ndgrid(Mod.x,Mod.y,Mod.z);
0010 Mesh.node=[X(:) Y(:) Z(:)];
0011 Mesh.cell=[];
0012 nx=length(Mod.x);ny=length(Mod.y);nz=length(Mod.z);
0013 nxy=nx*ny;
0014 <span class="comment">% uninode=[1 2 7 3;1 3 7 4;1 7 8 4;6 7 1 2;6 5 1 7;5 8 1 7];</span>
0015 uninode=[1 2 4 5;2 3 4 7;2 4 5 7;6 5 7 2;5 8 7 4];
0016 <span class="comment">% uninode=uninode(:,[1 2 4 3]);</span>
0017 nt=size(uninode,1);
0018 aa=repmat(Mod.M(:)',nt,1);
0019