femwind/hzc2icl.m

   1 function [icl1,icl2]=hzc2icl(hzc,n)
   2 % [icl1,icl2]=hzc2icl(hzc,n)
   3 % In:
   4 %   hzc     vector size 2, coarsening factors
   5 %   n       vector size 3, fine mesh size
   6 % Out:
   7 %   icl1    coarse indices in the horizontal direction 1
   8 %   icl2    coarse indices in the horizontal direction 2
   9 icl1=1:hzc(1):n(1);
  10 if icl1(end) ~= n(1)
  11     icl1(end+1)=n(1); % extend, end must be coarse
  12 end
  13 icl2=1:hzc(2):n(2);
  14 if icl2(end) ~= n(2)
  15     icl2(end+1)=n(2);
  16 end
  17 end