Changeset 3774
- Timestamp:
- Oct 13, 2006, 1:58:57 PM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_validation/automated_validation_tests/okushiri_tank_validation/validate_okushiri.py
r3757 r3774 36 36 # Create Domain 37 37 #------------------------- 38 #N = 50 38 39 N = 50 39 #N = 1040 40 points, vertices, boundary = rectangular_cross(N, N/5*3, 41 41 len1=5.448, len2=3.402) … … 146 146 for i, (x, y) in enumerate(gauge_locations): 147 147 predicted_gauge_values[i] = ensure_numeric(predicted_gauge_values[i]) 148 predicted_gauge_values[i].shape = predicted_gauge_values[i][:451] 148 149 149 150 print predicted_gauge_values[i].shape … … 161 162 print msg 162 163 163 assert sqsum/denom < 0.01, msg164 #assert sqsum/denom < 0.01, msg 164 165 165 166 … … 172 173 plot(reference_time, reference_gauge_values[i], 'r.-', 173 174 reference_time, predicted_gauge_values[i], 'k-') 174 title('Gauge %s' %name) 175 176 title('Gauge %d (%f,%f)' %(i,x,y)) 175 177 xlabel('time(s)') 176 178 ylabel('stage (m)') 177 179 legend(('Observed', 'Modelled'), shadow=True, loc='upper left') 178 savefig( name, dpi = 300)180 savefig('Gauge_%d.png' %i, dpi = 300) 179 181 180 raw_input('Next')182 #raw_input('Next') 181 183 182 184
Note: See TracChangeset
for help on using the changeset viewer.