quicwind/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 % in:
   6 %     n           number of cells in each direction
   7 %     h           step size in each direction
   8 %     v           vertical stretch factor
   9 % out:
  10 %     Mmat        matrix computing fluxes in normal directions on mesh
  11 %
  12 if ~exist('v','var')
  13     v=1;
  14 end
  15 zz = zeros(1,n(3));
  16 zz(1) = 0;
  17 for i=2:n(3)+1
  18     zz(i) = zz(i-1) + h(3) * v^(i-2);
  19 end
  20 [x,y,z] = ndgrid(h(1)*[0:n(1)],h(2)*[0:n(2)],zz);
  21 X = {x,y,z};