diff --git a/pj.mod b/pj.mod
index 4002e11..a79dbdc 100644
--- a/pj.mod
+++ b/pj.mod
@@ -77,8 +77,13 @@ FUNCTION invl(mean (ms)) (ms) {
 	}
 }
 VERBATIM
+#ifndef NRN_VERSION_GTEQ_8_2_0
 double nrn_random_pick(void* r);
 void* nrn_random_arg(int argpos);
+#define RANDCAST
+#else
+#define RANDCAST (Rand*)
+#endif
 ENDVERBATIM
 
 FUNCTION erand() {
@@ -89,7 +94,7 @@ VERBATIM
 		: each instance. However, the corresponding hoc Random
 		: distribution MUST be set to Random.negexp(1)
 		*/
-		_lerand = nrn_random_pick(_p_donotuse);
+		_lerand = nrn_random_pick(RANDCAST _p_donotuse);
 	}else{
 		/* only can be used in main thread */
 		if (_nt != nrn_threads) {
diff --git a/readme.txt b/readme.txt
index 48a1a69..01d2744 100644
--- a/readme.txt
+++ b/readme.txt
@@ -1,6 +1,6 @@
 NEURON files from the paper:
 
-Using Strahler’s analysis to reduce up to 200-fold the run time 
+Using Strahler’s analysis to reduce up to 200-fold the run time 
 of realistic neuron models, by A. Marasco, A.Limongiello, and M.Migliore
 (2013) Sci. Rep. 3, 2934; DOI:10.1038/srep02934 (in press).
 
@@ -16,5 +16,8 @@ Reduced Models" button, then select e4cb3a1.hoc in the morphologies
 folder. The simulation will reproduce the traces in Fig.2D (left) of
 the paper (where the last 1000 ms are shown). 
 
+Note: pj.mod updated to be valid for 9.0 and legacy NEURON versions.
+
 Questions on the code should be addressed to 
-Alessandro Limongiello: ax.magnum@gmail.com
\ No newline at end of file
+Alessandro Limongiello: ax.magnum@gmail.com
+