Changeset 1838
- Timestamp:
- Sep 16, 2005, 11:49:45 AM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
inundation/validation/Completed/LWRU2/extract_timeseries.py
r1835 r1838 4 4 import sys 5 5 from os import sep 6 from caching import cache 7 6 8 sys.path.append('..'+sep+'..'+sep) 7 9 … … 59 61 #Read model output 60 62 #filename = 'output.sww' 61 filename = 'lwru2.sww' 62 f = file_function(filename, 63 quantities = 'stage', 64 interpolation_points = gauges, 65 verbose = True) 63 #filename = 'lwru2.sww' 64 filename = 'lwru2_.05s.sww' 65 66 #f = file_function(filename, 67 # quantities = 'stage', 68 # interpolation_points = gauges, 69 # verbose = True) 70 71 f = cache(file_function,filename, 72 {'quantities': 'stage', 73 'interpolation_points': gauges, 74 'verbose': True}, 75 verbose = True) 66 76 67 77 … … 99 109 ion() 100 110 hold(False) 101 plot(reference_time, gauge_values[k], ' -k',102 reference_time, model, ' -b')111 plot(reference_time, gauge_values[k], 'r.-', 112 reference_time, model, 'k-') 103 113 title('Gauge %s' %name) 114 xlabel('time(s)') 115 ylabel('stage (m)') 116 legend(('Observed', 'Modelled'), shadow=True, loc='upper left') 104 117 savefig(name, dpi = 300) 105 118
Note: See TracChangeset
for help on using the changeset viewer.