matlab program for wind profile by adil rasool

>> radius= 10

>> z=3

>> alphaa=0

>> lamda=7

>> r= 0.1*radius

>> tr=(1/z)*(16*pi)*(r)*(sind((1/3)*(atand(radius/(lamda*r)))))^2

>> trr=tr/10

>> x=[0*trr, 0.2*trr, 0.4*trr, 0.6*trr, 0.8*trr, 1.0*trr]

>> yup=[0*trr, 0.1*trr, 0.11*trr, 0.09*trr, 0.05*trr, 0*trr]

>> ydown=[0*trr, 0.03*trr, 0.02*trr, 0.013*trr, 0.008*trr, 0*trr]

>> h1=line(x,yup)

>> hold

>> h2=line(x,ydown)

>> dir=[0 0 1]

>> alphaconnected=((2/3)*(atand(radius/lamda*r)))-alphaa

>> rotate(h1,dir,alphaconnected)

>> rotate(h2,dir,alphaconnected)