function dw=der2(z)
x=0.1;
u0=5;
k=0.4;
zoo=2e-5;
y=1;
a=.8;
Lz=10;
Lx=300;
if x<0.5
      us=a*sin(2*pi*x);
      dus=a*2*pi*cos(2*pi*x)/Lx;
      
else;
      us=0;
      dus=0;
      
end;
zo=zoo*(1-.9*sin(2*pi*x));
dzo=-zoo*.9*2*pi*cos(2*pi*x)/Lx;
dlg=dzo/zo;
up=us+(u0-us).*log(z./zo)/log(10/zoo);
upp=u0/log(10/zoo)*log(z/zoo);
u=upp.*(1-exp(-z/Lz))+up./exp(z/Lz);
dw=-dus.*(1-log(z/zo)/log(10/zoo))-(u0-us)/log(10/zoo).*dlg;
dw=dw/exp(z/Lz);

   
   
   