quicwind/block_example.m

   1 function block_example
   2 disp('A simple block example: input wind speed is zero inside, uniform outside')
   3 n=[500,500,50];
   4 h=[2,2,1];
   5 w=[1,1,1];
   6 u0=grad3z(zeros(n),h);  % staggered mesh wind arrays
   7 u0{3}(:,:,1)=0;       % zero vertical speed on the ground
   8 u0{1}(:,:,:)=1;       % unit homogeneous speed in direction 1
   9 u0{1}(200:300,200:300,1:20)=0; % except zero in a block
  10 u0{2}(:,:,:)=0;
  11 u = mass_cons_int(u0,h,w,'check');
  12 plot_wind(u,h,5)
  13 disp('horizontal velocity component about the middle of the leading edge of the block')
  14 squeeze(u{1}(197:203,250,25:-1:1))'
  15 disp('note small inside and speedup above the top of the block')
  16 disp('vertical velocity component about the middle of the leading edge of the block')
  17 squeeze(u{3}(197:203,250,25:-1:1))'
  18 disp('note positive and zero at ground level')
  19 disp('vertical velocity component about the middle of the trailing edge of the block')
  20 squeeze(u{3}(297:303,250,25:-1:1))'
  21 disp('note negative and zero at ground level')
  22 end