Changeset 3507

Timestamp:

Aug 17, 2006, 6:20:53 PM (18 years ago)

Author:

ole

Message:

Work post suresh kumar visit

File:

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

development/stochastic_study/plot_spread.py

@@ -40 +40 @@
 #hold(False)
 
-
-for j in range(1):
+N=250
+for j in range(N):
     #print j, filenames[j]
     plot(time, data[:,j], 'k-')
@@ -54 +54 @@
     #raw_input('Next')
 
-
-raw_input('Stats')    
+title('%d realisations' %N)
+savefig('realisations', dpi = 300)
+#raw_input('Stats')    
  
 
@@ -61 +62 @@
 hold(False)
 # Plot spread of stage values for each timestep
-for i in range(0,project.number_of_timesteps):
+for i in range(0, project.number_of_timesteps, 10):
     # Plot histogram#
 
@@ -68 +69 @@
     xlabel('realisations')
     ylabel('stage (m)')
-    title('Study %s: spread at timestep %d of %d (t=%.2f)'\
-          %(study, i, project.number_of_timesteps-1, time[i]))        
-raw_input('Next')
+    #title('Study %s: spread at timestep %d of %d (t=%.2f)'\
+    #      %(study, i, project.number_of_timesteps-1, time[i]))
+    title('Spread at timestep %d of %d (t=%.2f)'\
+          %(i, project.number_of_timesteps-1, time[i]))
+
+    savefig('spread_%d' %i, dpi = 300)    
+    #raw_input('Next')
 
 show()