vis3d/wind_uv_at_h.m

   1 function p=wind_uv_at_h(p,heights)
   2 % p=add2struct_wind_at_h(p,levels)
   3 % read
   4 % in:
   5 % p           structure from wrfatm2struct
   6 % heights     1D array, heights above terrain to compute the wind
   7 % out: additional fields in p computed, value(:,:,i) is value at heights(i)
   8 % p.uch,p.vch         wind at cell centers (theta points, not staggered)
   9 % p.uah,p.vah         wind at u-points, v-points (staggered)
  10
  11 if any(p.z0(:)<=0),
  12     error('roughness height z0 must be positive')
  13 end
  14
  15 % log interpolation of the wind at center points to height
  16 p.uch=log_interp_vert(p.uc,p.height,p.z0,heights);
  17 p.vch=log_interp_vert(p.vc,p.height,p.z0,heights);
  18
  19 p.heights=heights;
  20
  21 end