_RL wsqr,oldflux,newflux,entrainrate
oldflux=StartingFlux
newflux=oldflux+e2*radius*Wd(k2)*Dd(k2)*
entrainrate = (newflux - oldflux)/newflux
oldflux = newflux
real pdensity,wsqr,oldflux,newflux,entrainrate,adtemp
newflux=oldflux+e2*radius*Wd(k2)*Dd(k2)*0.50*