_RL wsqr,oldflux,newflux,entrainrate
newflux=oldflux+e2*radius*Wd(k2)*Dd(k2)*
PlumeEntrainment(k2+1) = newflux/StartingFlux
IF(newflux.LE.0.0) then
entrainrate = (newflux - oldflux)/newflux
oldflux = newflux
radius=sqrt(newflux/(Wd(k2)*Dd(k2)))
real pdensity,wsqr,oldflux,newflux,entrainrate,adtemp
newflux=oldflux+e2*radius*Wd(k2)*Dd(k2)*0.50*
IF(newflux.LT.0.0) then