From 365d1190a2caaf1004f1f2ee09a35528afdd85d5 Mon Sep 17 00:00:00 2001
From: sjplimp <sjplimp@f3b2605a-c512-4ea7-a41b-209d697bcdaa>
Date: Tue, 4 Aug 2015 19:26:51 +0000
Subject: [PATCH] git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@13830
 f3b2605a-c512-4ea7-a41b-209d697bcdaa

---
 src/REPLICA/prd.cpp | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/src/REPLICA/prd.cpp b/src/REPLICA/prd.cpp
index 4f03981295..a99c201480 100644
--- a/src/REPLICA/prd.cpp
+++ b/src/REPLICA/prd.cpp
@@ -333,11 +333,11 @@ void PRD::command(int narg, char **arg)
     // decrement clock by random time at which 1 or more events occurred
 
     int frac_t_event = t_event;
-    for (int i = 0; i <= fix_event->ncoincident; i++) {
+    for (int i = 0; i < fix_event->ncoincident; i++) {
       int frac_rand = static_cast<int> (random_clock->uniform() * t_event);
       frac_t_event = MIN(frac_t_event,frac_rand);
     }
-    int decrement = frac_t_event*universe->nworlds;
+    int decrement = (t_event - frac_t_event)*universe->nworlds;
     clock -= decrement;
 
     // share event across replicas
@@ -508,6 +508,7 @@ void PRD::dephase()
       if (compute_event->compute_scalar() > 0.0) {
         fix_event->restore_state_dephase();
         update->ntimestep -= t_dephase;
+        log_event();
       } else {
         fix_event->restore_state_quench();
         done = 1;