No commit message

No commit message

da/cosmo/statevector.py

@@ -288,14 +288,15 @@ class StateVector(object):
         # Create members 1:nmembers and add to ensemble_members list
 
         for member in range(1, self.nmembers):
+            rands = np.random.uniform(low=-1., high=1., size=self.nparams)
  #           rands = np.random.randn(self.nparams)
 #            rands = np.clip(rands,-1,1)
-            if (member==1):
-                rands = np.array([1.836451,1.91572,1.719371,1.947495,1.788597,1.109507,0.9401712,2,0.003201536,1.2276,0])-1.
-            elif (member==2):
-                rands = np.array([2,0.04773442,1.809479,1.136452,2,0.8753765,1.488298,2,1.908293,0.02784704,0])-1.
-            elif (member==3):
-                rands = np.array([0.9136019,1.386204,2,0.1250395,1.815127,2,0.3002012,1.666945,0,0,0])-1.
+#            if (member==1):
+ #               rands = np.array([1.836451,1.91572,1.719371,1.947495,1.788597,1.109507,0.9401712,2,0.003201536,1.2276,0])-1.
+  #          elif (member==2):
+   #             rands = np.array([2,0.04773442,1.809479,1.136452,2,0.8753765,1.488298,2,1.908293,0.02784704,0])-1.
+    #        elif (member==3):
+     #           rands = np.array([0.9136019,1.386204,2,0.1250395,1.815127,2,0.3002012,1.666945,0,0,0])-1.
 
             newmember = EnsembleMember(member)
             newmember.param_values = np.dot(C, rands) + newmean