Changeset 2975

Timestamp:

May 25, 2006, 5:21:14 PM (18 years ago)

Author:

ole

Message:

Stochastic work

Location:

development/stochastic_study

Files:

• plot_spread.py (modified) (4 diffs)
• project.py (modified) (1 diff)

development/stochastic_study/plot_spread.py

@@ -14 +14 @@
 
 # Initialise
-gauge_name = project.gauge_names[0]
+gauge = '%s.txt' %project.gauge_names[0]
 number_of_realisations = project.number_of_realisations
 
@@ -25 +25 @@
 j = 0 # Count realisations
 for filename in os.listdir('.'):
-    if filename.startswith(project.basename) and filename.endswith('%s.txt' %gauge_name):
+    if filename.startswith(project.basename) and filename.endswith(gauge):
         if j < data.shape[1]:        
             print 'Reading filename %s (column %d)' %(filename, j)
@@ -51 +51 @@
 for j in range(number_of_realisations):
     #print j, data[:,j]
-    plot(data[:,j], 'k-')
+    plot(time, data[:,j], 'k-')
     
-    raw_input('Next')
     xlabel('time(s)')
     ylabel('stage (m)')
-    title('Realisation %d of %d' %(j, number_of_realisations))    
-    
+    title('Realisation %d of %d' %(j, number_of_realisations-1))    
+
+    raw_input('Next')
+
 
 
@@ -68 +69 @@
     bins = create_bins(w, project.number_of_bins)
     print 'bins', bins
-    hist = histogram(w, bins)
-    print 'hist', hist
+    hist = histogram(w, bins, relative=True)
+    #print 'hist', hist
     
     #plot(w, 'k.')
-    plot(hist, 'k.')
+    plot(bins, hist, 'k.')
     xlabel('stage (m)')
-    ylabel('count')
+    ylabel('frequency')
     title('Timestep %d of %d (t=%.2f)'\
-          %(i, project.number_of_timesteps, time[i]))        
+          %(i, project.number_of_timesteps-1, time[i]))        
     raw_input('Next')
 

development/stochastic_study/project.py

@@ -19 +19 @@
 
 # Stats (Suresh ?)
-number_of_realisations = 1
+number_of_realisations = 9
 #std_dev = 0.0026  #Range is 26.035 cm
 std_dev = 0.0013  #Range is 26.035 cm