Changeset 4877


Ignore:
Timestamp:
Dec 6, 2007, 6:26:08 PM (16 years ago)
Author:
nick
Message:

update conical island validation

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_validation/conical_island/run_circular.py

    r4800 r4877  
    1111from anuga.shallow_water import Dirichlet_boundary, Time_boundary
    1212from anuga.shallow_water.data_manager import get_maximum_inundation_data
    13 from anuga.abstract_2d_finite_volumes.util import file_function
     13from anuga.abstract_2d_finite_volumes.util import file_function, gauges_sww2csv
    1414from anuga.pmesh.mesh_interface import create_mesh_from_regions
    1515from anuga.utilities.polygon import read_polygon, plot_polygons
     
    100100r2 = 2.5
    101101d=0.75
    102 res=0.0005
     102res=5.
    103103
    104104poly = []
     
    118118
    119119poly.append(poly1)
    120 internal_poly.append([poly1,.1])
     120internal_poly.append([poly1,.1*res])
    121121
    122122poly2 = []
     
    128128   
    129129poly.append(poly2)
    130 internal_poly.append([poly2,.01])
     130internal_poly.append([poly2,.01*res])
    131131
    132132poly3 = []
     
    138138   
    139139poly.append(poly3)
    140 internal_poly.append([poly3,1])
     140internal_poly.append([poly3,1*res])
    141141
    142142poly4 = []
     
    148148   
    149149poly.append(poly4)
    150 internal_poly.append([poly2,.01])
     150internal_poly.append([poly2,.01*res])
    151151
    152152plot_polygons(poly,'poly.png')
     
    174174    fileName = "test.csv"
    175175    file = open(fileName,"w")
    176     file.write("x,y,elevation \n")
     176    file.write("easting,northing,elevation \n")
    177177    x = 12.96
    178178    y = 13.80
     
    257257    domain.write_time()
    258258   
    259 for t in domain.evolve(yieldstep = 1, finaltime = 75
     259for t in domain.evolve(yieldstep = 1, finaltime = 50
    260260                       ,skip_initial_step = True):
    261261    domain.write_time()
     
    274274#d=2
    275275#poly = {}
     276
    276277for i,angle in enumerate(angles):
    277278#    angle = ((angle-180)/180.0)*pi
     
    288289   
    289290    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.