diff --git a/readme.html b/readme.html
index 4162162..6a9334f 100644
--- a/readme.html
+++ b/readme.html
@@ -17,4 +17,6 @@
 20120330 CaIntraCellDyn.mod had it's solve method updated from euler
 to derivimplicit.  For more information see:
 http://www.neuron.yale.edu/phpBB/viewtopic.php?f=28&t=592
+
+20241125 Updated hoc files to use Random123
 </pre></html>
diff --git a/wdr-complete-model-exportsyns.hoc b/wdr-complete-model-exportsyns.hoc
index 53983cb..b9868d6 100755
--- a/wdr-complete-model-exportsyns.hoc
+++ b/wdr-complete-model-exportsyns.hoc
@@ -179,6 +179,7 @@ fih = new FInitializeHandler("loadqueue()")
 proc loadqueue() { local ii, jj, event_time  localobj r
         
     r = new Random(123456789) // use a different seed if you need different random streams; "123456789" is the seed used in the paper figure
+    r.Random123()
     
     //load Ad-fiber spike times
     for ii=0,stim_times.size()-1 {
@@ -223,6 +224,7 @@ objref ampa_interneuron_conn[N_IC], nmda_interneuron_conn[N_IC], gabaa_interneur
 
 objref r
 r = new Random(123456789) // use a different seed if you need different random streams; the seed used for the figures in the paper was 123456789
+r.Random123()
 
 for i=0,N_IC-1 {
     
diff --git a/wdr-complete-model.hoc b/wdr-complete-model.hoc
index 2ad7ff2..9061f24 100755
--- a/wdr-complete-model.hoc
+++ b/wdr-complete-model.hoc
@@ -64,6 +64,7 @@ for i=0,n_stim_sets-1 {
 /*
 objref r, stim_times
 r = new Random(666)
+r.Random123()
 r.negexp(T2)
 stim_times  = new Vector(n_stim_sets, 0.0)
 last_time = start_time
@@ -184,6 +185,7 @@ fih = new FInitializeHandler("loadqueue()")
 proc loadqueue() { local ii, jj, event_time  localobj r
         
     r = new Random(123456789) // use a different seed if you need different random streams; "123456789" is the seed used in the paper figures
+    r.Random123()
     
     //load Ad-fiber spike times
     for ii=0,stim_times.size()-1 {
@@ -228,6 +230,7 @@ objref ampa_interneuron_conn[N_IC], nmda_interneuron_conn[N_IC], gabaa_interneur
 
 objref r
 r = new Random(123456789) // use a different seed if you need different random streams; the seed used for the figures in the paper was 123456789
+r.Random123()
 
 for i=0,N_IC-1 {