Changeset 4877

Timestamp:

Dec 6, 2007, 6:26:08 PM (17 years ago)

Author:

nick

Message:

update conical island validation

File:

• anuga_validation/conical_island/run_circular.py (modified) (10 diffs)

anuga_validation/conical_island/run_circular.py

@@ -11 +11 @@
 from anuga.shallow_water import Dirichlet_boundary, Time_boundary
 from anuga.shallow_water.data_manager import get_maximum_inundation_data
-from anuga.abstract_2d_finite_volumes.util import file_function
+from anuga.abstract_2d_finite_volumes.util import file_function, gauges_sww2csv
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
 from anuga.utilities.polygon import read_polygon, plot_polygons
@@ -100 +100 @@
 r2 = 2.5
 d=0.75
-res=0.0005
+res=5.
 
 poly = []
@@ -118 +118 @@
 
 poly.append(poly1)
-internal_poly.append([poly1,.1])
+internal_poly.append([poly1,.1*res])
 
 poly2 = []
@@ -128 +128 @@
     
 poly.append(poly2)
-internal_poly.append([poly2,.01])
+internal_poly.append([poly2,.01*res])
 
 poly3 = []
@@ -138 +138 @@
     
 poly.append(poly3)
-internal_poly.append([poly3,1])
+internal_poly.append([poly3,1*res])
 
 poly4 = []
@@ -148 +148 @@
     
 poly.append(poly4)
-internal_poly.append([poly2,.01])
+internal_poly.append([poly2,.01*res])
 
 plot_polygons(poly,'poly.png')
@@ -174 +174 @@
     fileName = "test.csv"
     file = open(fileName,"w")
-    file.write("x,y,elevation \n")
+    file.write("easting,northing,elevation \n")
     x = 12.96
     y = 13.80
@@ -257 +257 @@
     domain.write_time()
    
-for t in domain.evolve(yieldstep = 1, finaltime = 75
+for t in domain.evolve(yieldstep = 1, finaltime = 50
                        ,skip_initial_step = True): 
     domain.write_time()
@@ -274 +274 @@
 #d=2
 #poly = {}
+
 for i,angle in enumerate(angles):
 #    angle = ((angle-180)/180.0)*pi
@@ -288 +289 @@
     
     print 'maximum_inundation_data',angle, run_up, x_y 
-    
-
-
-
-
-
+
+
+
+
+
+
+