Index: /issm/trunk-jpl/src/c/shared/Random/randomgenerator.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/Random/randomgenerator.cpp	(revision 26596)
+++ /issm/trunk-jpl/src/c/shared/Random/randomgenerator.cpp	(revision 26597)
@@ -2,20 +2,2 @@
  * \brief random number generating functions
  */
-
-#include <iostream>
-#include "./randomgenerator.h"
-
-#undef M_PI
-#define M_PI 3.141592653589793238462643
-
-//rnd_uniform_distribution::rnd_uniform_distribution(){/*{{{*/
-
-//		a   = 1103515245; 	// BSD Formula
-//		c  = 12345;					// BSD Formula
-//		m = 2147483648;			// BSD Formula
-//		_seed = 0;
-//		lbound = 0.0;
-//		ubound = 1.0;
-//		return;
-//}
-/*}}}*/
Index: /issm/trunk-jpl/src/c/shared/Random/randomgenerator.h
===================================================================
--- /issm/trunk-jpl/src/c/shared/Random/randomgenerator.h	(revision 26596)
+++ /issm/trunk-jpl/src/c/shared/Random/randomgenerator.h	(revision 26597)
@@ -2,72 +2,2 @@
  * \brief prototypes for randomgenerator.h
  */
-
-//#ifndef _RANDOMGENERATOR_H_
-//#define _RANDOMGENERATOR_H_
-
-#undef M_PI
-#define M_PI 3.141592653589793238462643
-
-class rnd_uniform_distribution
-{
-
-  private:
-    int a;
-    int c;
-    unsigned int m;
-    unsigned _seed;
-    double lbound;
-    double ubound;
-
-  public:
-
-    /*constructors, destructors: */
-    rnd_uniform_distribution() : _seed( 0 ), a( 1103515245 ), c( 12345 ), m( 2147483648 ), lbound(0.0), ubound(1.0) {}
-    rnd_uniform_distribution(double a_1,double a_2) : _seed( 0 ), a( 1103515245 ), c( 12345 ), m( 2147483648 ), lbound(a_1), ubound(a_2) {}
-    ~rnd_uniform_distribution(){}
-
-    void seed( unsigned int s ) { _seed = s; }
-    unsigned int get_seed() { return _seed; }
-    double generator(){
-      _seed = ( a * _seed + c ) % m;
-      return (ubound-lbound)*(double) _seed/ m + lbound;
-    }
-
-};
-
-class rnd_normal_distribution
-{
-
-  private:
-    unsigned _seed;
-    double mean;
-    double sdev;
-
-  public:
-
-    /*constructors, destructors: */
-    rnd_normal_distribution() : _seed( 0 ), mean( 0), sdev(1.0) {}
-    rnd_normal_distribution(double m,double s) : _seed( 0 ), mean( m ), sdev(s) {}
-    ~rnd_normal_distribution(){}
-
-    void seed( unsigned int s ) { _seed = s; }
-    double generator()
-    {
-      rnd_uniform_distribution unifdistri;
-      unifdistri.seed(_seed);
-
-      double u1 = unifdistri.generator();
-      double u2 = unifdistri.generator();
-
-      double R = sqrt(-2*log(u1));
-      double theta = 2*M_PI*u2;
-
-      seed(unifdistri.get_seed());
-
-      return mean + sdev * (R*cos(theta));
-
-    }
-
-};
-
-//#endif //ifndef _RANDOMGENERATOR_H_
Index: /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_sampling.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_sampling.cpp	(revision 26596)
+++ /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_sampling.cpp	(revision 26597)
@@ -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);
 	}