0071         <span class="keyword">if</span> Mod.nx(k+1)&gt;Mod.nx(k)+1, Mod.nx(k+1)=Mod.nx(k)+1; <span class="keyword">end</span>
0072         <span class="keyword">if</span> Mod.ny(k+1)&gt;Mod.ny(k)+1, Mod.ny(k+1)=Mod.ny(k)+1; <span class="keyword">end</span>
0073     <span class="keyword">end</span>
0074 <span class="keyword">case</span> 3, <span class="comment">% quadratic</span>
0075     Mod.nx=round(diff(Mod.z)/Mod.dx);
0076     Mod.ny=round(diff(Mod.z)/Mod.dx);
0077 <span class="keyword">otherwise</s