在编程计算颗粒碰撞时,某一时间步长结束时,得出两个颗粒的速度,如何将这一速度在下一时间步长进行调用呢?

 upnew(n1)=up(n1)*(1-(1+erc)*rratiom)
 upnew(n2)=(m(n1)*up(n1)+m(n2)*up(n2)-
     &    m(n1)*upnew(n1))/m(n2)
 ratiov=up(n1)/vp(n2)
 slid=3.5*uf*(1+erc)
 if(ratiov.lt.slid) then
   vpnew(n1)=vp(n1)*(1-2/(7*rratiom))
 else
   vpnew(n1)=vp(n1)*(1-uf*(1+erc)*ratiov*rratiom)
 endif
 vpnew(n2)=(m(n1)*vp(n1)+m(n2)*vp(n2)-
     &  m(n1)*vpnew(n1))/m(n2)
 wpnew(n1)=vpnew(n1)*atan(IP*180/pi)
 wpnew(n2)=vpnew(n2)*atan(IP*180/pi)