Changeset 4555 for anuga_validation/automated_validation_tests/okushiri_tank_validation/compare_timeseries_with_measures.py

Timestamp:

Jun 21, 2007, 3:23:11 PM (18 years ago)

Author:

ole

Message:

Computed maximal runup height and location in automated okushiri validation and
added one simple assertion. More could be done here.

File:

• anuga_validation/automated_validation_tests/okushiri_tank_validation/compare_timeseries_with_measures.py (modified) (2 diffs)

anuga_validation/automated_validation_tests/okushiri_tank_validation/compare_timeseries_with_measures.py

@@ -141 +141 @@
     #                            rtol=rtol, atol=atol)
     if plotting is False:
-        assert allclose(reference_value, computed_value, rtol=rtol, atol=atol), msg
+        assert allclose(reference_value, computed_value,
+                        rtol=rtol, atol=atol), msg
     
 
@@ -215 +216 @@
 
 
+# Check max runup
+
+from anuga.shallow_water.data_manager import get_maximum_inundation_elevation
+from anuga.shallow_water.data_manager import get_maximum_inundation_location
+from anuga.utilities.polygon import is_inside_polygon
+
+q = get_maximum_inundation_elevation(sww_filename)
+loc = get_maximum_inundation_location(sww_filename)
+
+
+print 'Max runup elevation: ', q
+print 'Max runup location:  ', loc
+
+
+from create_okushiri import gulleys
+assert is_inside_polygon(loc, gulleys)
+
+# FIXME more asserts here
+
+
+
+#msg = 'We got %f, should have been %f' %(q, q_max)
+#assert allclose(q, q_max, rtol=1.0/N), msg
+##print 'loc', loc, q
+#assert allclose(-loc[0]/2, q) # From topography formula 
+