From 83d75279ae773491eecd48ca864b7deed16d634d Mon Sep 17 00:00:00 2001
From: sjplimp <sjplimp@f3b2605a-c512-4ea7-a41b-209d697bcdaa>
Date: Tue, 7 Jan 2014 21:03:02 +0000
Subject: [PATCH] git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@11162
 f3b2605a-c512-4ea7-a41b-209d697bcdaa

---
 src/GRANULAR/pair_gran_hooke_history.cpp | 10 +++++++---
 1 file changed, 7 insertions(+), 3 deletions(-)

diff --git a/src/GRANULAR/pair_gran_hooke_history.cpp b/src/GRANULAR/pair_gran_hooke_history.cpp
index ccebe88f0e..83e232b11e 100644
--- a/src/GRANULAR/pair_gran_hooke_history.cpp
+++ b/src/GRANULAR/pair_gran_hooke_history.cpp
@@ -474,14 +474,18 @@ void PairGranHookeHistory::init_style()
   int itype;
   for (i = 1; i <= atom->ntypes; i++) {
     onerad_dynamic[i] = onerad_frozen[i] = 0.0;
-    if (ipour >= 0) 
+    if (ipour >= 0) {
+      itype = i;
       onerad_dynamic[i] = 
         *((double *) modify->fix[ipour]->extract("radius",itype));
-    if (idep >= 0) 
+    }
+    if (idep >= 0) {
+      itype = i;
       onerad_dynamic[i] = 
         *((double *) modify->fix[idep]->extract("radius",itype));
+    }
   }
-    
+
   double *radius = atom->radius;
   int *mask = atom->mask;
   int *type = atom->type;