Changeset 5401 for anuga_work/development/boxingday08/run_boxingday_urs.py
- Timestamp:
- Jun 13, 2008, 10:29:14 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/development/boxingday08/run_boxingday_urs.py
r5393 r5401 40 40 # Application specific imports 41 41 import project # Definition of file names and polygons 42 import create_boundary43 42 44 43 #------------------------------------------------------------------------------ … … 53 52 54 53 timestamp = strftime('%Y%m%d_%H%M%S',localtime()) 55 basename = scenario + '_urs'54 basename = scenario[:4] + '_urs' 56 55 57 56 … … 70 69 # resolutions (maximal area of per triangle) for each polygon 71 70 #------------------------------------------------------------------------------ 72 extent_res = 1000000 0.071 extent_res = 1000000.0 73 72 contour20m_res = 50000.0 74 73 island_res = 5000.0 75 74 bay_res = 2000.0 76 patong_res = 100.075 patong_res = 400.0 77 76 78 77 # make polygon that contains land that does not affect result. … … 80 79 interior_regions = [[project.patong, patong_res], 81 80 [project.bay, bay_res], 82 [project.contour20m, contour20m_res] ,83 [project.island_north, island_res],84 [project.island_south, island_res],85 [project.island_south2, island_res]]81 [project.contour20m, contour20m_res]]#, 82 #[project.island_north, island_res], 83 #[project.island_south, island_res], 84 #[project.island_south2, island_res]] 86 85 87 86 #for coarse run to test gauges 88 extent_res = 10000000.089 interior_regions=None87 #extent_res = 10000000.0 88 #interior_regions=None 90 89 91 90 from Numeric import arange,allclose 92 print 'Beware not all boundary segments that are tagged as ocean in polyline are tagged here. missing are 79,80,81,82,83,84 which are set to other ocean' 93 boundary_tags={'ocean': arange(1,79).tolist(), 'otherocean': [0,89,84,83,82,81,80,79], 'land': [86, 87], 'both': [85,88]} 94 95 96 #trigs_min = number_mesh_triangles(interior_regions, create_boundary.bounding_polygon,extent_res) 91 boundary_tags={'ocean': arange(0,41).tolist(), 'otherocean': [41,44], 'land': [42,43]} 92 93 94 #trigs_min = number_mesh_triangles(interior_regions, project.bounding_polygon,extent_res) 97 95 #print 'Minimum number of traingles ', trigs_min 98 96 … … 104 102 105 103 _ = cache(create_mesh_from_regions, 106 create_boundary.bounding_polygon,104 project.bounding_polygon, 107 105 {'boundary_tags': boundary_tags, 108 106 'maximum_triangle_area': extent_res, … … 110 108 'interior_regions': interior_regions}, 111 109 verbose = True, 112 #dependencies = [' create_boundary.py']110 #dependencies = ['project.py'] 113 111 dependencies = ['run_boxingday_urs.py'] 114 112 #, evaluate=True … … 131 129 print domain.statistics() 132 130 133 domain.set_name(basename+ 'friction'+timestamp)131 domain.set_name(basename+timestamp) 134 132 #domain.set_name(basename) 135 133 domain.set_datadir(scenario) … … 208 206 # start when boundary stage first becomes non-zero. 209 207 210 for t in domain.evolve(yieldstep = 5.0, finaltime = 18000.0,208 for t in domain.evolve(yieldstep = 20.0, finaltime = 18000.0, 211 209 skip_initial_step = False): 212 210 if allclose(t,10800.0):
Note: See TracChangeset
for help on using the changeset viewer.