Changeset 3035 for development/stochastic_study/plot_spread.py

Timestamp:

May 31, 2006, 5:03:59 PM (18 years ago)

Author:

ole

Message:

Working with Suresh at ACFR

File:

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

development/stochastic_study/plot_spread.py

@@ -15 +15 @@
 # Initialise
 # Read in all realisations-timeseries
-study = 'cyclone1'
-study = 'cyclone3'
-#study = 'nautilus3'
+
+#study = 'nautilus1' #~70 realisations, blocks of 100, sequential, stddev= 0.0006
+#study = 'nautilus3' #~100 realisations, blocks of 10, sequential, stddev= 0.0013 
+study = 'cyclone1' #~4000 realisations, blocks of 100, parallel, stddev= 0.0006
+#study = 'cyclone3' #~120 realisations, blocks of 10, parallel, stddev= 0.0013
 
 time, data, filenames = read_realisations(study, 
                                           #max_realisations = 200,
-                                          gauge_number=2,
-                                          use_cache=False)
+                                          gauge_number=0,
+                                          sorting='randomised',
+                                          #sorting='numerical',
+                                          verbose=True,
+                                          use_cache=True)
 number_of_realisations = data.shape[1]
+number_of_timesteps = data.shape[0]
 print 'Read %d realisations' %number_of_realisations
+
 
 # Plot
@@ -34 +41 @@
 
 
-for j in range(20):
+for j in range(1):
     #print j, filenames[j]
     plot(time, data[:,j], 'k-')
@@ -45 +52 @@
           %(study, shortname, j, number_of_realisations-1))
     
-    #title('Study %s: timeseries for %s (realisation %d of %d)'\
-    #      %(study, filenames[j], j, number_of_realisations-1))    
-    raw_input('Next')
+    #raw_input('Next')
+
+
+raw_input('Stats')    
  
 
@@ -53 +61 @@
 hold(False)
 # Plot spread of stage values for each timestep
-for i in range(300,320): #project.number_of_timesteps):
+for i in range(0,project.number_of_timesteps):
     # Plot histogram#
 
@@ -62 +70 @@
     title('Study %s: spread at timestep %d of %d (t=%.2f)'\
           %(study, i, project.number_of_timesteps-1, time[i]))        
-    raw_input('Next')
+raw_input('Next')
 
+show()
+import sys; sys.exit()
 
 hold(False)