Changeset 4877
 Timestamp:
 Dec 6, 2007, 6:26:08 PM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

anuga_validation/conical_island/run_circular.py
r4800 r4877 11 11 from anuga.shallow_water import Dirichlet_boundary, Time_boundary 12 12 from anuga.shallow_water.data_manager import get_maximum_inundation_data 13 from anuga.abstract_2d_finite_volumes.util import file_function 13 from anuga.abstract_2d_finite_volumes.util import file_function, gauges_sww2csv 14 14 from anuga.pmesh.mesh_interface import create_mesh_from_regions 15 15 from anuga.utilities.polygon import read_polygon, plot_polygons … … 100 100 r2 = 2.5 101 101 d=0.75 102 res= 0.0005102 res=5. 103 103 104 104 poly = [] … … 118 118 119 119 poly.append(poly1) 120 internal_poly.append([poly1,.1 ])120 internal_poly.append([poly1,.1*res]) 121 121 122 122 poly2 = [] … … 128 128 129 129 poly.append(poly2) 130 internal_poly.append([poly2,.01 ])130 internal_poly.append([poly2,.01*res]) 131 131 132 132 poly3 = [] … … 138 138 139 139 poly.append(poly3) 140 internal_poly.append([poly3,1 ])140 internal_poly.append([poly3,1*res]) 141 141 142 142 poly4 = [] … … 148 148 149 149 poly.append(poly4) 150 internal_poly.append([poly2,.01 ])150 internal_poly.append([poly2,.01*res]) 151 151 152 152 plot_polygons(poly,'poly.png') … … 174 174 fileName = "test.csv" 175 175 file = open(fileName,"w") 176 file.write(" x,y,elevation \n")176 file.write("easting,northing,elevation \n") 177 177 x = 12.96 178 178 y = 13.80 … … 257 257 domain.write_time() 258 258 259 for t in domain.evolve(yieldstep = 1, finaltime = 75259 for t in domain.evolve(yieldstep = 1, finaltime = 50 260 260 ,skip_initial_step = True): 261 261 domain.write_time() … … 274 274 #d=2 275 275 #poly = {} 276 276 277 for i,angle in enumerate(angles): 277 278 # angle = ((angle180)/180.0)*pi … … 288 289 289 290 print 'maximum_inundation_data',angle, run_up, x_y 290 291 292 293 294 295 291 292 293 294 295 296 297
Note: See TracChangeset
for help on using the changeset viewer.