femwind/regular_mesh.m

   1 function X=regular_mesh(n,h,v)
   2 % X=regular_mesh(n,h)
   3 % Creates a regular mesh that is stretched in the vertical direction but
   4 % uniform horizationally.
   5 % Use plot_mesh_3d(X) to see the result.
   6 % in:
   7 %     n           size 3 number of cells in each direction
   8 %     h           size 3 step size in each direction
   9 %     v           size 1 vertical stretch factor
  10 % out:
  11 %     X           cell array with x y z grid coordinates
  12 %
  13 if ~exist('v','var')
  14     v=1;
  15 end
  16 zz = zeros(1,n(3)+1);
  17 for i=1:n(3)
  18     zz(i+1) = zz(i) + h(3) * v^(i-1);
  19 end
  20 [x,y,z] = ndgrid(h(1)*[0:n(1)],h(2)*[0:n(2)],zz);
  21 X = {x,y,z};
  22 fprintf('regular_mesh nodes %i %i %i step %g %g %g expansion %g size %g %g %g\n',...
  23     n,h,v,n(1)*h(1),n(2)*h(2),zz(end))
  24 end