--- trunk/OOPSE/libmdtools/do_Forces.F90 2004/01/05 22:12:11 895 +++ trunk/OOPSE/libmdtools/do_Forces.F90 2004/01/05 22:18:52 897 @@ -4,7 +4,7 @@ !! @author Charles F. Vardeman II !! @author Matthew Meineke -!! @version $Id: do_Forces.F90,v 1.42 2004-01-05 22:12:11 chuckv Exp $, $Date: 2004-01-05 22:12:11 $, $Name: not supported by cvs2svn $, $Revision: 1.42 $ +!! @version $Id: do_Forces.F90,v 1.43 2004-01-05 22:18:52 chuckv Exp $, $Date: 2004-01-05 22:18:52 $, $Name: not supported by cvs2svn $, $Revision: 1.43 $ module do_Forces use force_globals @@ -228,11 +228,12 @@ contains real( kind = DP ) :: rijsq real(kind=dp),dimension(3) :: d real(kind=dp) :: rfpot, mu_i, virial - integer :: me_i + integer :: me_i, me_j logical :: is_dp_i integer :: neighborListSize integer :: listerror, error integer :: localError + integer :: propPack_i, propPack_j real(kind=dp) :: listSkin = 1.0 @@ -262,14 +263,14 @@ contains #ifdef IS_MPI if (.not.allocated(propertyMapI)) then - allocate(propertyMapI(5,getNrow()) + allocate(propertyMapI(5,nrow)) endif do i = 1, nrow me_i = atid_row(i) #else if (.not.allocated(propertyMapI)) then - allocate(propertyMapI(5,getNlocal()) + allocate(propertyMapI(5,nlocal)) endif do i = 1, natoms @@ -297,14 +298,14 @@ contains #ifdef IS_MPI if (.not.allocated(propertyMapJ)) then - allocate(propertyMapJ(5,getNcol()) + allocate(propertyMapJ(5,ncol)) endif do j = 1, ncol me_j = atid_col(j) #else if (.not.allocated(propertyMapJ)) then - allocate(propertyMapJ(5,getNlocal()) + allocate(propertyMapJ(5,nlocal)) endif do j = 1, natoms