4 |
|
|
5 |
|
!! @author Charles F. Vardeman II |
6 |
|
!! @author Matthew Meineke |
7 |
< |
!! @version $Id: do_Forces.F90,v 1.7 2003-04-01 16:50:14 chuckv Exp $, $Date: 2003-04-01 16:50:14 $, $Name: not supported by cvs2svn $, $Revision: 1.7 $ |
7 |
> |
!! @version $Id: do_Forces.F90,v 1.8 2003-04-04 19:57:01 mmeineke Exp $, $Date: 2003-04-04 19:57:01 $, $Name: not supported by cvs2svn $, $Revision: 1.8 $ |
8 |
|
|
9 |
|
module do_Forces |
10 |
|
use force_globals |
246 |
|
|
247 |
|
!! save current configuration, construct neighbor list, |
248 |
|
!! and calculate forces |
249 |
< |
call saveNeighborList(q) |
249 |
> |
call saveNeighborList(nlocal, q) |
250 |
|
|
251 |
|
neighborListSize = size(list) |
252 |
|
nlist = 0 |
313 |
|
|
314 |
|
! save current configuration, contruct neighbor list, |
315 |
|
! and calculate forces |
316 |
< |
call saveNeighborList(q) |
316 |
> |
call saveNeighborList(natoms, q) |
317 |
|
|
318 |
|
neighborListSize = size(list) |
319 |
|
|
543 |
|
|
544 |
|
if ( is_DP_i .and. is_DP_j ) then |
545 |
|
|
546 |
< |
call do_dipole_pair(i, j, d, r, rijsq, pot, u_l, f, t, & |
546 |
> |
call do_dipole_pair(natoms,i, j, d, r, rijsq, pot, u_l, f, t, & |
547 |
|
do_pot, do_stress) |
548 |
|
if (FF_uses_RF .and. SimUsesRF()) then |
549 |
|
call accumulate_rf(i, j, r, u_l) |