% all units cgs
u=50;       % cm/s
d50=.0250;  %cm  (250 micron)
zr=100;
z0=1;  % cm
c_d=z0tocd(z0,zr);
ustar=u*sqrt(c_d);
ustar_c=1.38;
q_b=mpm(d50,ustar_c,ustar);

