Changeset 2986

Timestamp:

May 26, 2006, 3:53:42 PM (19 years ago)

Author:

ole

Message:

More stochastic plotting

File:

• development/stochastic_study/plot_spread.py (modified) (2 diffs)

development/stochastic_study/plot_spread.py

@@ -50 +50 @@
 # Simple plot of timeseries for different realisations 
 for j in range(number_of_realisations):
-    #print j, data[:,j]
+    print j, data[:,j]
     plot(time, data[:,j], 'k-')
     
     xlabel('time(s)')
     ylabel('stage (m)')
-    title('Realisation %d of %d' %(j, number_of_realisations-1))    
+    title('Realisation %d of %d' %(j, number_of_realisations-1))    #
+#raw_input('Next')
 
+
+hold(False)
+# Plot spread of stage values for each timestep
+for i in range(200,project.number_of_timesteps):
+    # Plot histogram
+
+    w = data[i,:]
+    plot(w, 'k.')
+    xlabel('realisations')
+    ylabel('stage (m)')
+    title('Timestep %d of %d (t=%.2f)'\
+          %(i, project.number_of_timesteps-1, time[i]))        
     raw_input('Next')
 
 
-
+hold(False)
 # Plot histogram of stage values for each timestep
-for i in range(project.number_of_timesteps):
+for i in range(200,project.number_of_timesteps):
     # Plot histogram
 
@@ -78 +91 @@
     title('Timestep %d of %d (t=%.2f)'\
           %(i, project.number_of_timesteps-1, time[i]))        
-    raw_input('Next')
+    #raw_input('Next')