Changeset 3010
- Timestamp:
- May 29, 2006, 2:53:46 PM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
development/stochastic_study/plot_spread.py
r3006 r3010 3 3 4 4 # Standard modules 5 import sys, os 5 import sys, os, string 6 6 7 7 # Related major packages … … 15 15 # Initialise 16 16 # Read in all realisations-timeseries 17 time, data, filenames = read_realisations('cyclone1', 17 #study = 'cyclone1' 18 study = 'nautilus3' 19 20 time, data, filenames = read_realisations(study, 18 21 max_realisations = 200, 19 22 use_cache=True) 20 #time, data, filenames = read_realisations('test', exclude='_0_',21 # max_realisations = 200,22 # use_cache=True)23 #time, data, filenames = read_realisations('nautilus1', exclude='_0_')24 23 number_of_realisations = data.shape[1] 25 24 … … 32 31 #hold(False) 33 32 34 # Simple plot of timeseries for different realisations 33 35 34 for j in range(20): 36 print j, filenames[j]35 #print j, filenames[j] 37 36 plot(time, data[:,j], 'k-') 38 37 39 38 xlabel('time(s)') 40 39 ylabel('stage (m)') 41 title('Timeseries for realisation %d of %d' %(j, number_of_realisations-1)) 40 41 shortname = string.join(os.path.split(filenames[j])[-1].split('_')[-2:], '_') 42 title('Study %s: timeseries for %s (realisation %d of %d)'\ 43 %(study, shortname, j, number_of_realisations-1)) 44 45 #title('Study %s: timeseries for %s (realisation %d of %d)'\ 46 # %(study, filenames[j], j, number_of_realisations-1)) 42 47 raw_input('Next') 43 48 … … 53 58 xlabel('realisations') 54 59 ylabel('stage (m)') 55 title('S pread at timestep %d of %d (t=%.2f)'\56 %( i, project.number_of_timesteps-1, time[i]))60 title('Study %s: spread at timestep %d of %d (t=%.2f)'\ 61 %(study, i, project.number_of_timesteps-1, time[i])) 57 62 raw_input('Next') 58 63 … … 74 79 xlabel('stage (m)') 75 80 ylabel('frequency') 76 title(' Histogram at timestep %d of %d (t=%.2f)'\77 %( i, project.number_of_timesteps-1, time[i]))81 title('Study %s: histogram at timestep %d of %d (t=%.2f)'\ 82 %(study, i, project.number_of_timesteps-1, time[i])) 78 83 raw_input('Next') 79 84
Note: See TracChangeset
for help on using the changeset viewer.