Changeset 4799

Timestamp:

Nov 7, 2007, 3:42:51 PM (17 years ago)

Author:

nick

Message:

working on validation

File:

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

anuga_validation/conical_island/run_circular.py

@@ -103 +103 @@
 
 poly = []
+poly1 = []
 internal_poly=[]
-for i,angle in enumerate(angles):
+for i,angle in enumerate(angles1):
+    #convert Degs to Rads
     angle = ((angle-180)/180.0)*pi
 #    p1 = get_xy(x,y,r1,angle-0.01745*d)
@@ -112 +114 @@
 #    poly[i] = [[p1[0],p1[1]],[p2[0],p2[1]],[p3[0],p3[1]],[p4[0],p4[1]]]
 
-    p = get_xy(x,y,4,angle)
-    poly.insert(0,[p[0],p[1]])
-internal_poly.insert(0,(poly,1))
-
-poly1 = []
-for i,angle in enumerate(angles):
+    p = get_xy(x,y,3.6,angle)
+    poly1.append([p[0],p[1]])
+
+poly.append(poly1)
+internal_poly.append([poly1,.1])
+
+poly2 = []
+for i,angle in enumerate(angles1):
+    #convert Degs to Rads
     angle = ((angle-180)/180.0)*pi
     p = get_xy(x,y,2.,angle)
-    poly1.insert(0,[p[0],p[1]])
-    
-internal_poly.insert(1,(poly1,.1))
-
-poly2 = []
-for i,angle in enumerate(angles):
+    poly2.append([p[0],p[1]])
+    
+poly.append(poly2)
+internal_poly.append([poly2,.01])
+
+poly3 = []
+for i,angle in enumerate(angles1):
+    #convert Degs to Rads
     angle = ((angle-180)/180.0)*pi
-    p = get_xy(x,y,1.,angle)
-    poly2.insert(0,[p[0],p[1]])
-    
-internal_poly.insert(2,(poly2,1))
-
-print internal_poly
-
-plot_polygons(internal_poly,'poly.png',)
+    p = get_xy(x,y,1.1,angle)
+    poly3.append([p[0],p[1]])
+    
+poly.append(poly3)
+internal_poly.append([poly3,1])
+
+plot_polygons(poly,'poly.png')
 
 poly_all= [[0,7],[0,25],[30,25],[30,7]]
-#poly1=[[9,10],[9,17.6],[16.6,17.6],[16.6,10]]
-#internal_poly.insert(0,(poly1,0.025))
 create_mesh_from_regions(poly_all,
                          boundary_tags={'wall': [0,2],'wave': [3], 'buffer':[1]},