Index: /issm/trunk-jpl/src/c/shared/Random/randomgenerator.h
===================================================================
--- /issm/trunk-jpl/src/c/shared/Random/randomgenerator.h	(revision 26590)
+++ /issm/trunk-jpl/src/c/shared/Random/randomgenerator.h	(revision 26591)
@@ -20,6 +20,4 @@
     double ubound;
 
-    int drnd() { return( _seed = ( a * _seed + c ) % m ); }
-
   public:
 
@@ -30,5 +28,8 @@
     void seed( unsigned int s ) { _seed = s; }
     unsigned int get_seed() { return _seed; }
-    double generator() { return (ubound-lbound)*(double) drnd()/ m + lbound; }
+    double generator(){
+      _seed = ( a * _seed + c ) % m;
+      return (ubound-lbound)*(double) _seed/ m + lbound;
+    }
 
 };
Index: /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_sampling.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_sampling.cpp	(revision 26590)
+++ /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_sampling.cpp	(revision 26591)
@@ -18,5 +18,5 @@
 
 	/*Define seed*/
-	rnd_normal_distribution distribution;
+	//rnd_normal_distribution distribution;
 	if(seed<0){
 		std::random_device rd;
@@ -29,5 +29,5 @@
 		seed = seed + 783728*my_rank; // change default seed for parallel simulations (by considering an arbitrary shif based on the rank number)
 	}
-	distribution.seed(seed);
+	//distribution.seed(seed);
 
 	int        *local_indices = NULL;
@@ -38,5 +38,5 @@
 	ppf->GetLocalSize(&M);
 	for(int i=0;i<M;i++){
-		rdnumber = distribution.generator();
+		rdnumber = 1.0;//distribution.generator();
 		ppf->SetValue(local_indices[i],rdnumber,INS_VAL);
 	}