4 |
|
|
5 |
|
!! @author Charles F. Vardeman II |
6 |
|
!! @author Matthew Meineke |
7 |
< |
!! @version $Id: do_Forces.F90,v 1.65 2004-05-27 19:51:18 tim Exp $, $Date: 2004-05-27 19:51:18 $, $Name: not supported by cvs2svn $, $Revision: 1.65 $ |
7 |
> |
!! @version $Id: do_Forces.F90,v 1.66 2004-05-28 15:21:37 gezelter Exp $, $Date: 2004-05-28 15:21:37 $, $Name: not supported by cvs2svn $, $Revision: 1.66 $ |
8 |
|
|
9 |
|
module do_Forces |
10 |
|
use force_globals |
428 |
|
nAtomsInCol = getNatomsInCol(plan_atom_col) |
429 |
|
nGroupsInRow = getNgroupsInRow(plan_group_row) |
430 |
|
nGroupsInCol = getNgroupsInCol(plan_group_col) |
431 |
+ |
write(*,*) nAtomsInRow, nAtomsInCol, nGroupsInRow, nGroupsInCol |
432 |
+ |
write(*,*) pot_local |
433 |
|
#else |
434 |
|
natoms = nlocal |
435 |
|
#endif |
479 |
|
endif |
480 |
|
|
481 |
|
do loop = loopStart, loopEnd |
482 |
+ |
write(*,*) 'doign loop ', loop |
483 |
|
|
484 |
|
! See if we need to update neighbor lists |
485 |
|
! (but only on the first time through): |
486 |
|
if (loop .eq. loopStart) then |
487 |
< |
call checkNeighborList(nGroups, q_group, listSkin, update_nlist) |
487 |
> |
#ifdef IS_MPI |
488 |
> |
call checkNeighborList(nGroupsInRow, q_group_row, listSkin, & |
489 |
> |
update_nlist) |
490 |
> |
#else |
491 |
> |
call checkNeighborList(nGroups, q_group, listSkin, & |
492 |
> |
update_nlist) |
493 |
> |
#endif |
494 |
|
endif |
495 |
|
|
496 |
|
if (update_nlist) then |
497 |
|
!! save current configuration and construct neighbor list |
498 |
|
#ifdef IS_MPI |
499 |
< |
call saveNeighborList(nGroupsInRow, q_group) |
499 |
> |
call saveNeighborList(nGroupsInRow, q_group_row) |
500 |
|
#else |
501 |
|
call saveNeighborList(nGroups, q_group) |
502 |
|
#endif |
511 |
|
iend = nGroups - 1 |
512 |
|
#endif |
513 |
|
outer: do i = istart, iend |
514 |
+ |
|
515 |
+ |
write(*,*) 'doing ', i |
516 |
|
|
517 |
|
if (update_nlist) point(i) = nlist + 1 |
518 |
|
|