4 |
|
|
5 |
|
!! @author Charles F. Vardeman II |
6 |
|
!! @author Matthew Meineke |
7 |
< |
!! @version $Id: do_Forces.F90,v 1.49 2004-04-15 16:18:26 tim Exp $, $Date: 2004-04-15 16:18:26 $, $Name: not supported by cvs2svn $, $Revision: 1.49 $ |
7 |
> |
!! @version $Id: do_Forces.F90,v 1.50 2004-04-22 21:33:55 tim Exp $, $Date: 2004-04-22 21:33:55 $, $Name: not supported by cvs2svn $, $Revision: 1.50 $ |
8 |
|
|
9 |
|
module do_Forces |
10 |
|
use force_globals |
421 |
|
end if |
422 |
|
call zero_work_arrays() |
423 |
|
|
424 |
+ |
|
425 |
|
do_pot = do_pot_c |
426 |
|
do_stress = do_stress_c |
427 |
|
|
871 |
|
#endif |
872 |
|
|
873 |
|
|
873 |
– |
|
874 |
|
end subroutine do_force_loop |
875 |
|
|
876 |
|
subroutine do_pair(i, j, rijsq, d, do_pot, do_stress, u_l, A, f, t, pot) |
888 |
|
real ( kind = dp ), intent(inout) :: d(3) |
889 |
|
integer :: me_i, me_j |
890 |
|
|
891 |
+ |
|
892 |
|
r = sqrt(rijsq) |
893 |
|
|
894 |
|
#ifdef IS_MPI |