Changeset 4800 for anuga_validation/conical_island/run_circular.py
 Timestamp:
 Nov 7, 2007, 4:09:45 PM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

anuga_validation/conical_island/run_circular.py
r4799 r4800 124 124 #convert Degs to Rads 125 125 angle = ((angle180)/180.0)*pi 126 p = get_xy(x,y,2. ,angle)126 p = get_xy(x,y,2.4,angle) 127 127 poly2.append([p[0],p[1]]) 128 128 … … 134 134 #convert Degs to Rads 135 135 angle = ((angle180)/180.0)*pi 136 p = get_xy(x,y,1. 1,angle)136 p = get_xy(x,y,1.6,angle) 137 137 poly3.append([p[0],p[1]]) 138 138 139 139 poly.append(poly3) 140 140 internal_poly.append([poly3,1]) 141 142 poly4 = [] 143 for i,angle in enumerate(angles1): 144 #convert Degs to Rads 145 angle = ((angle180)/180.0)*pi 146 p = get_xy(x,y,1.1,angle) 147 poly4.append([p[0],p[1]]) 148 149 poly.append(poly4) 150 internal_poly.append([poly2,.01]) 141 151 142 152 plot_polygons(poly,'poly.png')
Note: See TracChangeset
for help on using the changeset viewer.