--- trunk/mdtools/md_code/simulation_module.F90	2003/01/27 21:28:19	249
+++ trunk/mdtools/md_code/simulation_module.F90	2003/01/28 22:16:55	252
@@ -82,7 +82,7 @@ contains
     real(kind = dp ), intent(in), dimension(3) :: box
     real(kind = dp ), intent(in) :: rlist
     real(kind = dp ), intent(in) :: rcut
-
+    integer :: alloc_stat
     if( setsim ) return  ! simulation is already initialized
     setSim = .true.
 
@@ -92,7 +92,10 @@ contains
     thisSim%rcut         = rcut
     thisSim%rcutsq       = rcut * rcut
     thisSim%rcut6        = thisSim%rcutsq * thisSim%rcutsq * thisSim%rcutsq
-
+    
+    if (.not. allocated(q0)) then
+       allocate(q0(3,nLRParticles),stat=alloc_stat)
+    endif
   end subroutine setSimulation
 
   function getNparticles() result(nparticles)
@@ -137,7 +140,7 @@ contains
     
 #ifdef MPI 
     dispmx_tmp = 0.0E0_DP
-    do i = 1, nlocal
+    do i = 1, thisSim%nLRparticles
        dispmx_tmp = max( abs ( q(1,i) - q0(1,i) ), dispmx )
        dispmx_tmp = max( abs ( q(2,i) - q0(2,i) ), dispmx )
        dispmx_tmp = max( abs ( q(3,i) - q0(3,i) ), dispmx )
@@ -145,7 +148,8 @@ contains
     call mpi_allreduce(dispmx_tmp,dispmx,1,mpi_double_precision, & 
        mpi_max,mpi_comm_world,mpi_err)
 #else
-    do i = 1, natoms
+
+    do i = 1, thisSim%nLRparticles
        dispmx = max( abs ( q(1,i) - q0(1,i) ), dispmx )
        dispmx = max( abs ( q(2,i) - q0(2,i) ), dispmx )
        dispmx = max( abs ( q(3,i) - q0(3,i) ), dispmx )
@@ -225,7 +229,7 @@ contains
     
     thisrcut = thisSim%rcut
     if(present(rcut2)) rcut2 = thisSim%rcutsq
-    if(present(rcut2)) rcut6 = thisSim%rcut6
+    if(present(rcut6)) rcut6 = thisSim%rcut6
 
   end subroutine getRcut