_RL Ad(Nr),Wd(Nr),Dd(Nr),Md(Nr)
Wd(k)= - .5*(wVel(K)+wVel(K+1))
wsqr=Wd(k)*Wd(k)
StartingFlux=radius*radius*Wd(k)*Dd(k)
newflux=oldflux+e2*radius*Wd(k2)*Dd(k2)*
Wd(k2+1)=sqrt(wsqr)
radius=sqrt(newflux/(Wd(k2)*Dd(k2)))
Md(k)=Wd(k)*Ad(k)
& 'S/R OPPS_CALC: Md, Wd, and k are ',Md(k2),Wd(k2),k2
IF ( Wd(K2) .NE. 0. _d 0 ) dt = min(dt,drF(k2)/Wd(k2))
REAL Ad(km),Sd(km),Td(km),Wd(km),Dd(km),Md(km)
Wd(k) = 0.03
newflux=oldflux+e2*radius*Wd(k2)*Dd(k2)*0.50*
write(0,*)"Md, Wd, and k are ",Md(k2),Wd(k2),k2