function [tau,reff]=mie2lognormal(N0,lambda)
r=logspace(-2.5,1.5,100);
rm1=0.05;
rm2=0.8;
sig1=2;
sig2=2;
N=N0(1)*lognormal(r,rm1,sig1)+N0(2)*lognormal(r,rm2,sig2);
reff=trapz(r,N.*r.^3)/trapz(r,N.*r.^2);
n=ref_sulfate(lambda);
tau=mie2size(n,lambda,N,r);

return
semilogx(r,N.*r.^2)

