From b0a5f4a2805404a18e775d657904ba6f49b04d82 Mon Sep 17 00:00:00 2001
From: sjplimp <sjplimp@f3b2605a-c512-4ea7-a41b-209d697bcdaa>
Date: Tue, 18 Nov 2014 15:41:05 +0000
Subject: [PATCH] git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@12717
 f3b2605a-c512-4ea7-a41b-209d697bcdaa

---
 src/velocity.cpp | 19 ++++++++++---------
 1 file changed, 10 insertions(+), 9 deletions(-)

diff --git a/src/velocity.cpp b/src/velocity.cpp
index 7e6ef94d54..b264af0fcb 100644
--- a/src/velocity.cpp
+++ b/src/velocity.cpp
@@ -229,13 +229,14 @@ void Velocity::create(double t_desired, int seed)
 
     for (i = 1; i <= natoms; i++) {
       if (dist_flag == 0) {
-        vx = random->uniform();
-        vy = random->uniform();
-        vz = random->uniform();
+        vx = random->uniform() - 0.5;
+        vy = random->uniform() - 0.5;
+        vz = random->uniform() - 0.5;
       } else {
         vx = random->gaussian();
         vy = random->gaussian();
         vz = random->gaussian();
+        sum += vx;
       }
       m = atom->map(i);
       if (m >= 0 && m < nlocal) {
@@ -264,9 +265,9 @@ void Velocity::create(double t_desired, int seed)
     for (i = 0; i < nlocal; i++) {
       if (mask[i] & groupbit) {
         if (dist_flag == 0) {
-          vx = random->uniform();
-          vy = random->uniform();
-          vz = random->uniform();
+          vx = random->uniform() - 0.5;
+          vy = random->uniform() - 0.5;
+          vz = random->uniform() - 0.5;
         } else {
           vx = random->gaussian();
           vy = random->gaussian();
@@ -289,9 +290,9 @@ void Velocity::create(double t_desired, int seed)
       if (mask[i] & groupbit) {
         random->reset(seed,x[i]);
         if (dist_flag == 0) {
-          vx = random->uniform();
-          vy = random->uniform();
-          vz = random->uniform();
+          vx = random->uniform() - 0.5;
+          vy = random->uniform() - 0.5;
+          vz = random->uniform() - 0.5;
         } else {
           vx = random->gaussian();
           vy = random->gaussian();