reflct.f

Go to the documentation of this file.

      SUBROUTINE reflct(I,RI2)
*
*
*       Boundary reflection.
*       --------------------
*
      include 'common6.h'
*
*
      rdot = x(1,i)*xdot(1,i) + x(2,i)*xdot(2,i) + x(3,i)*xdot(3,i)
*
*       Avoid multiple reflections for inward moving particles.
      IF (rdot.LT.0.0) go to 50
*
*       Obtain radial velocity with respect to the inertial centre.
      ri = sqrt(ri2)
      rdot = rdot/ri
*
*       Form radial velocity for mirror reflection.
      rdot = -2.0*rdot
*
*       Set reflected velocity at the boundary and initialize X0DOT.
      DO 10 k = 1,3
          xdot(k,i) = xdot(k,i) + rdot*x(k,i)/ri
          x0dot(k,i) = xdot(k,i)
   10 CONTINUE
*
*       Predict current coordinates & velocities for all neighbours.
      nnb = list(1,i)
      CALL xvpred(i,nnb)
*
*       Initialize force polynomials & time-steps.
      CALL fpoly1(i,i,0)
      CALL fpoly2(i,i,0)
*
      ri2 = -1.0
*       Return of negative argument indicates new polynomials & steps.
*
   50 RETURN
*
      END