DO LND_i=2,LND_iend
LND_x(LND_i)=LND_x(LND_i-1)+LND_dx
LND_tmp = log(LND_x(LND_i))-LND_mu
LND_pdf(LND_i)= 1.
& / (LND_x(LND_i)*LND_sigma*sqrt(2*3.1416))
DO LND_i=1,LND_iend
IF ( LND_x(LND_i).GT.Hlimit(k) ) k=k+1
& + LND_pdf(LND_i)*LND_dx
& + LND_pdf(LND_i)*LND_x(LND_i)*LND_dx
INTEGER LND_i, LND_iend