pan>, <span class="comment">%logarithmically ascending</span>
0078     Mod.nx=round(refx.*((0:K-1)*0.13+1));
0079     Mod.ny=round(refx.*((0:K-1)*0.13+1));
0080 <span class="keyword">end</span>
0081 Mod.ncells=0;
0082 <span class="keyword">for</span> k=1:K,
0083     Mod.M{k}=ones(floor(nnx/Mod.nx(k)),floor(nny/Mod.ny(k)))*Mod.Bg(k);
0084     Mod.ncells=Mod.ncells+prod(size(Mod.M{k}));
0085 <span class="keyword">end</span>
0086 message(sprintf(<span class="string">'Creating para model: %d cells, rho=%.1f'</s