Changeset 5401
- Timestamp:
- Jun 13, 2008, 10:29:14 AM (15 years ago)
- Location:
- anuga_work/development/boxingday08
- Files:
-
- 3 edited
-
run_boxingday_most.py (modified) (7 diffs)
-
run_boxingday_polyline.py (modified) (9 diffs)
-
run_boxingday_urs.py (modified) (8 diffs)
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/development/boxingday08/run_boxingday_most.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 + 'source'54 basename = scenario[:4] + '_most' 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.0 75 patong_res = 400.0 76 77 77 78 78 79 # make polygon that contains land that does not affect result. 79 80 interior_regions = [[project.patong, patong_res], 80 81 [project.bay, bay_res], 81 [project.contour20m, contour20m_res] ,82 [project.island_north, island_res],83 [project.island_south, island_res],84 [project.island_south2, island_res]]82 [project.contour20m, contour20m_res]]#, 83 #[project.island_north, island_res], 84 #[project.island_south, island_res], 85 #[project.island_south2, island_res]] 85 86 86 87 #for coarse run to test gauges 87 extent_res = 10000000.0 88 interior_regions=None88 #extent_res = 1000000.0 #causes inaccurate stage 89 #interior_regions=None 89 90 90 91 from Numeric import arange,allclose 91 boundary_tags={'ocean': arange( 1,85).tolist(), 'otherocean': [0,89], 'land': [86, 87], 'both': [85,88]}92 93 94 #trigs_min = number_mesh_triangles(interior_regions, create_boundary.bounding_polygon,extent_res)92 boundary_tags={'ocean': arange(0,41).tolist(), 'otherocean': [41,44], 'land': [42,43]} 93 94 95 #trigs_min = number_mesh_triangles(interior_regions, project.bounding_polygon,extent_res) 95 96 #print 'Minimum number of traingles ', trigs_min 96 97 … … 102 103 103 104 _ = cache(create_mesh_from_regions, 104 create_boundary.bounding_polygon,105 project.bounding_polygon, 105 106 {'boundary_tags': boundary_tags, 106 107 'maximum_triangle_area': extent_res, 107 108 'filename': meshname, 108 109 'interior_regions': interior_regions}, 109 verbose = True ,110 dependencies = ['run_boxingday_most.py']110 verbose = True#, 111 #dependencies = ['run_boxingday_most.py'] 111 112 #, evaluate=True 112 113 ) … … 128 129 print domain.statistics() 129 130 130 domain.set_name(basename+ 'friction'+timestamp)131 domain.set_name(basename+timestamp) 131 132 domain.set_datadir(scenario) 132 133 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum']) … … 197 198 print 'Available boundary tags', domain.get_boundary_tags() 198 199 Bf = File_boundary(project.boundary_most_out+'.sww', 199 domain, time_thinning=1, verbose = True)200 domain, time_thinning=1, use_cache=True,verbose = True) 200 201 Br = Reflective_boundary(domain) 201 202 Bd = Dirichlet_boundary([tide,0.0,0.0]) … … 219 220 # start when boundary stage first becomes non-zero. 220 221 221 for t in domain.evolve(yieldstep = 5.0, finaltime = 18000.0,222 for t in domain.evolve(yieldstep = 20.0, finaltime = 18000.0, 222 223 skip_initial_step = False): 223 224 if allclose(t,10800.0): -
anuga_work/development/boxingday08/run_boxingday_polyline.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 … … 54 53 55 54 timestamp = strftime('%Y%m%d_%H%M%S',localtime()) 56 basename = scenario + '_polyline'55 basename = scenario[:4] + '_polyline' 57 56 58 57 … … 71 70 # resolutions (maximal area of per triangle) for each polygon 72 71 #------------------------------------------------------------------------------ 73 extent_res = 1000000 0.072 extent_res = 1000000.0 74 73 contour20m_res = 50000.0 75 74 island_res = 5000.0 76 75 bay_res = 2000.0 77 patong_res = 100.076 patong_res = 400.0 78 77 79 78 … … 82 81 interior_regions = [[project.patong, patong_res], 83 82 [project.bay, bay_res], 84 [project.contour20m, contour20m_res] ,85 [project.island_north, island_res],86 [project.island_south, island_res],87 [project.island_south2, island_res]]83 [project.contour20m, contour20m_res]]#, 84 #[project.island_north, island_res], 85 #[project.island_south, island_res], 86 #[project.island_south2, island_res]] 88 87 89 88 90 89 #for coarse run to test gauges 91 extent_res = 10000000.092 interior_regions=None90 #extent_res = 10000000.0 91 #interior_regions=None 93 92 94 93 from Numeric import arange,allclose 95 boundary_tags={'ocean': arange( 1,85).tolist(), 'otherocean': [0,89], 'land': [86, 87], 'both': [85,88]}96 97 #trigs_min = number_mesh_triangles(interior_regions, create_boundary.bounding_polygon,extent_res)94 boundary_tags={'ocean': arange(0,41).tolist(), 'otherocean': [41,44], 'land': [42,43]} 95 96 #trigs_min = number_mesh_triangles(interior_regions, project.bounding_polygon,extent_res) 98 97 #print 'Minimum number of traingles ', trigs_min 99 98 … … 105 104 106 105 _ = cache(create_mesh_from_regions, 107 create_boundary.bounding_polygon,106 project.bounding_polygon, 108 107 {'boundary_tags': boundary_tags, 109 108 'maximum_triangle_area': extent_res, … … 111 110 'interior_regions': interior_regions}, 112 111 verbose = True, 113 dependencies = [' create_boundary.py']112 dependencies = ['project.py'] 114 113 #, evaluate=True 115 114 ) … … 130 129 print domain.statistics() 131 130 132 domain.set_name(basename+ 'friction'+timestamp)131 domain.set_name(basename+timestamp) 133 132 #domain.set_name(basename) 134 133 domain.set_datadir(scenario) … … 168 167 #------------------------------------------------------------------------------ 169 168 boundary_urs_in='data/boxing' 170 boundary_urs_out= create_boundary.base_name169 boundary_urs_out=project.base_name 171 170 172 171 print 'Available boundary tags', domain.get_boundary_tags() 173 172 Bf = File_boundary(boundary_urs_out+'.sts', 174 173 domain, time_thinning=1, 175 use_cache= False,verbose = True,boundary_polygon=create_boundary.bounding_polygon)174 use_cache=True,verbose = True,boundary_polygon=project.bounding_polygon) 176 175 177 176 Br = Reflective_boundary(domain) … … 196 195 # start when boundary stage first becomes non-zero. 197 196 198 for t in domain.evolve(yieldstep = 5.0, finaltime = 18000.0,197 for t in domain.evolve(yieldstep = 20.0, finaltime = 18000.0, 199 198 skip_initial_step = False): 200 199 if allclose(t,10800.0): -
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.
